Experimental studies on the changes of rock resistivity image and anisotropy
-
摘要: 在两组人工样品自由表面以中心点为基准对称布设3条辐射状测线，对样品实施单轴应力加载和卸载后，利用电阻率层析成像方法构建了相应的视电阻率相对变化图像，并计算和绘制了表征裂隙产生和发展速率的视电阻率各向异性系数λ*以及表征裂隙产生和发展方位的各向异性主轴方位角α随应力和深度的变化曲线.结果表明:所有测线所对应的RRC图像均随着应力的变化呈现出相同的变化趋势，即在加载阶段，随着应力的增加，视电阻率相对变化图像中电阻率降低区域逐渐收缩，而电阻率升高区域逐渐扩张，在卸载阶段，随着应力的减小，电阻率降低区域继续收缩，电阻率升高区域继续扩张；样品中的高阻体对其所在部位及附近区域的电阻率增幅有较大影响，而对横越高阻体测线的视电阻率相对变化图像的趋势性变化无影响；对于原始电性为各向异性的样品，随着应力的增加，其各向异性程度降低；裂隙主要在岩样的浅部产生和发展，而在较深部位的裂隙产生和发展的速率相对较低.上述结果有助于解释和理解地震、火山活动和大型地质构造运动引起的视电阻率及其各向异性的变化特征，电阻率层析成像方法可能成为目前地震电阻率观测方法的有益补充.Abstract: Apparent resistivity data was acquired during the uniaxial compression on two sets of man-made samples. Then we constructed the relative resistivity change (RRC) images corresponding to three radial measuring lines intersecting with the center of a sample surface using electrical resistivity tomography, and plotted the curves of apparent resistivity anisotropy factor λ* and azimuthal angle of anisotropy axis α versus stress and depth. λ* and α represent the rate and direction of crack generation and development respectively. Our results indicate that all RRC images show the same change trend with the change of stress. With the increase of stress, the resistivity-decreased region (RDR) in the RRC images would shrink gradually, while the resistivity-increased region (RIR) would expand gradually. During the process of unload-ing, with the decrease of stress, the RIR continues to expand, and RDR conti-nues to shrink. The high-resistivity block embedded in a sample has a great influence on the resistivity-increased amplitude at its location and surroundings, but little effect on the trending change of resistivity image. For the samples with originally electrical anisotropy, λ* decreases with the increase of stress; Cracks appeared and developed mainly in the shallower part of a rock sample, while in the deeper part, the rate of crack generation and development is much lower, which can help to explain and understand the changes in resistivity and its anisotropy caused by earthquakes, volcanic activities and large-scale tectonic movements. This method could be a useful complement to the current seismic resistivity observation methods.
-
引言
地震是构造应力作用下岩体突然断裂和错动所致。在这个过程中，随着应力的不断加载，地壳介质的物性参数也不断变化。及至孕震晚期，震源区应力相对集中的地壳介质的电性结构会发生变化，从而引起大地电场的变化（毛桐恩等，1999），这种地电场变化已在多次中强地震前观测到，属于中短期异常（Varotsos，Alexopoulos，1984；黄清华，刘涛，2006；马钦忠，2008；田山等，2009；安张辉等，2017；席继楼，2019）。但由于电磁环境复杂，有效的地震电信号经常被淹没在各种干扰中，使得震前异常电信号提取面临较大困难。目前常规使用的地震电信号提取方法有长短极距比值法（田山等，2009）、极化方位计算法（毛桐恩等，1999）、常规波形分析法（钱复业，赵玉林，2005）、频谱分析法（范莹莹等，2010）等，利用上述方法处理1975年海城MS7.3地震、2006年文安MS5.1地震、2008年汶川MS8.0地震、2013年芦山MS7.0地震和2017年九寨沟MS7.0地震等多次中强地震的观测资料时，均提取到与地震有关的地电场异常变化，但这些方法在异常提取、机理解释和抗干扰性方面还存在一定的局限。因此在复杂电磁环境中提取相对较弱的地震电磁信号，亟需将数学、信号处理与地震电磁物理过程相结合，才能更好地降低或消除地电场数据中干扰成分对预测的影响（黄清华，2005）。
近年来，相关学者对大地电场的机理、特征、数值模拟以及大地电场的物理解析等方面进行了深入的研究（黄清华，刘涛，2006；叶青等，2007），特别是谭大诚等（2010，2011，2012，2013，2019）基于大地电场的潮汐机理，将源于空间电流系和潮汐作用的大地电场与岩体裂隙结构联系起来（陈全等，2019），建立了大地电场的岩体裂隙水（电荷）渗流（移动）模型，并逐步发展出地电场优势方位角分析方法，藉此探寻震中附近的岩体裂隙结构变化。谭大诚等（2019）和艾萨等（2020）利用地电场优势方位角分析方法，在近几年的部分中强地震前也提取到了显著的异常信息。尤其是在部分地震前会出现多个场地的优势方位角变化在时间上出现了准同步的异常现象，这也加强了异常信息的信度。此外，地电场优势方位角分析方法机理解释明确，并对一般空间电磁干扰、附近电磁环境干扰等具有一定程度的抗干扰性（辛建村等，2017；张波等，2020），故该方法具有较高的可靠性和较好的应用前景。
2021年5月22日青海省果洛州玛多发生MS7.4地震（34.59°N，98.34°E），震源深度为17 km，震中位于巴颜喀拉地块内部，是一次典型的左旋走滑地震事件，局部兼有正断分量，在地表形成了走向北西—南东、长约70 km的破裂面。该地震是2008年汶川MS8.0地震以来中国大陆发生的震级最大的一次地震，打破了青海地区56个月的M6平静以及10年的M7平静。大武台地电场场地距本次地震震中约172 km，为距震中最近的地球物理场观测台站，但大武台的原始观测资料并未呈显著的地电异常现象。为了探究玛多MS7.4地震前大武台地电场是否存在地震电信号，本文拟利用大地电场优势方位角方法来提取玛多地震前大武台地电场的异常变化，对异常特征进行分析研究，并结合该方法在大武台记录的以往震例中的应用，梳理出强震前地电场优势方位角的异常变化特征，归纳总结该方法在青海地区强震前的异常判定指标，旨在为未来大震前后地电场的观测分析和跟踪研判提供参考依据。
1. 观测台站及数据分析方法
1.1 大武地震台
大武地震台位于青海省果洛藏族自治州玛沁县大武镇，海拔3 800 m，距离镇区3 km。目前大武台有新旧两个地电观测场地，相距约36 km。旧的大武台地电场（以下简称为大武台（旧））地位于台站院内，新的大武台地电场（以下简称为大武台（新））位于甘德县青珍乡的一个草滩，两套地电场台站周围均无高大建筑物、水库、湖泊、河流和大型厂矿等重大环境干扰源，且分别位于东昆仑断裂带两侧约10 km和25 km （图1），而东昆仑活动断裂带是一条规模巨大、构造活动十分强烈的全新世活动断裂，自1900年以来沿该断裂带已发生了10余次M7强震，最大地震为1947年3月果洛达日M7.7地震。
![]()
图 1 玛多MS7.4地震震中及附近地电场台站分布Figure 1. Location of the epicenter of Maduo MS7.4 earthquake and distribution of the geoelectric stations大武台地电场（旧）使用ZD9A-2B地电场仪，采用“多方向、多极距”方法进行观测，布极方式为“L”型，分别为NS，EW和NE三个测向的长、短极距共六道测线，线路采用绝缘铠装电缆以地埋的方式架设。电极装置采用固体不极化电极，深度为3.5 m，电极和铠装电缆均埋设在地表冻土层以下的潮湿土壤中，地电场布极装置方式如图2所示。
![]()
图 2 大武地震台（旧）地电场观测布极图Figure 2. Electrode distribution of electric field （old site） of Dawu seismic station1.2 分析方法
大地电场具有日变波形特征，目前大地电场来源于地球日月潮汐作用和空间地球电离层两部分。而地壳中岩体内总存在含水裂隙，这些裂隙水或水中的电荷以日为周期沿裂隙往返渗流或移动，建立了大地电场的日变化的岩体裂隙水（电荷）渗流（移动）模型。即地下介质的应力积累会导致岩体裂隙结构的变化，从而使大地电场的强度或方位发生改变，因此岩体裂隙结构的优势方位基本就是大地电场的优势方位。
在实际地电场观测中，设潮汐地电场（大地电场主要成分）为ET，场地裂隙水主体渗流方向为α，也称为地电场优势方位角。当台站地电场的NS测向与NW测向之间的相关性最高时，大地电场ET的优势方位角α（北偏东）的计算公式如下（谭大诚等，2012，2013，2019）：
$$ \alpha {\text{≈}} {180}－ \frac{180}{\pi} {{\rm{arctan}}}\left(\sqrt{2}\frac{\displaystyle\sum\limits _{i＝1}^{10}{A}_{{\rm{NW}} （ i ） }}{\displaystyle\sum\limits _{i＝1}^{10}{A}_{{\rm{NS}} （ i ） }}－1\right) ，$$ (1) 式中ANW（i）和ANS（i）分别为NW测向和NS测向的第i阶潮汐谐波振幅。应用谐波分析法对地电场分钟值数据（每天1 440个）开展信号处理，即对这些数据进行快速傅里叶变换，得到周期为24，12，8，6，4.8，4，3.4，3，2.7和2.4 h谐波的振幅（谭大诚等，2019）。
此外，本文分析数据采用原始分钟值数据，减少了预处理造成的干扰，真实地反映异常的变化情况。
2. 异常特征
利用大地电场优势方位角方法，基于2020年以来大武台地电场（新旧两台）观测数据，得到地电场优势方位角随时间的变化曲线如图3所示。由图可见：大武台（旧）得到的观测数据其优势方位角在2020年1—5月中旬变化不大，即Δα在20°以内跳变，而在2020年5月中旬以后，方位角跳变范围出现快速、持续偏转的趋势，其偏转范围最大接近45° （图3a），持续至2021年4月方位角跳变出现快速偏转回升，随后发生玛多MS7.4地震；大武台（新）得到的观测数据其优势方位角在2020年1—5月跳变范围稳定，即Δα跳动范围在10°以内，而在2020年5月中旬以后，方位角跳变范围也同步出现缓慢、持续偏转的趋势，其偏转范围增大接近90° （图3b），随后发震。从整个演化过程来看，两套地电场台的地电场优势方位角于2020年5月中旬同步出现显著的快速偏转现象，随后在震前两个月再次出现快速偏转回升至正常跳变范围，随后发震。
![]()
图 3 2020年以来大武台地电场旧（a）、新（b）台的地电场优势方位角变化Figure 3. The variation of the dominant azimuth of geoelectric field at old （a） and new （b） sites of Dawu station since 20203. 分析与讨论
为进一步探求玛多MS7.4地震前大地电场优势方位角异常的空间演化特征，选取2020年以来玛多地震震中500 km范围内的八个地电场台站的观测资料进行分析。首先查阅各观测台站的工作日志，对资料进行核实和排查，确保异常的可信度；然后使用地电场优势方位角方法得到这八个地电台的地电场优势方位角变化，详见表1，从表中可见四川省甘孜台、甘肃玛曲台同时期的地电场优势方位角均出现较为同步的异常变化，如图4所示。可见：甘孜台地电场优势方位角在2020年6月初开始出现显著、快速、持续偏转的现象，其偏转范围最大接近90°，收缩成近直线，随后在震前两个月出现跳变增大现象（图4a）；玛曲台在震前同时期也出现了异常现象，但因环境干扰较为严重，异常形态不显著，主要以持续突跳为主（图4b）。
表 1 玛多MS7.4地震前震中500 km内地电场优势方位角的变化特征Table 1. Variation characteristic of the dominant azimuth based on geoelectric field before MS7.4 Maduo earthquake within 500 km地电场台站 震中距/km 异常形态 异常出现日期
年-月-日开始观测日期
年-月-日大武（旧） 170 快速偏转，Δα≈45° 2020-05-10 2014-12-12 大武（新） 169 Δα增大至90° 2020-05-08 2020-01-01 玛曲 335 Δα增大至60° 2020-05-04 2014-03-22 甘孜 360 快速偏转，Δα≈90° 2020-05-15 2017-12-01 都兰 190 无异常 无 2014-12-01 白水河 355 无异常 无 2015-01-01 金银滩 256 无异常 无 2019-06-04 门源地电场 433 无异常 无 2020-01-01 ![]()
图 4 2020年以来甘孜（a）、玛曲台（b）的地电场方位角异常变化Figure 4. The variation of the dominant azimuth angle based on geoelectric field data at the seismic stations Garze （a） and Maqu （b） since 2020通过上述分析可以看出，地震前出现大地电场优势方位角异常现象是真实存在的，并且时间上具有一定准同步性，但场地具有选择性，例如：都兰台虽然距玛多地震震中仅190 km，但震前无异常现象，对该台有记录以来的数据分析也显示都兰台地电场使用该方法在以往地震前均无异常变化，目前认为这可能与其处于应力不敏感地区有关。根据詹艳（2008）利用大地电磁方法对巴颜喀拉地块探测的研究，大武台所处的巴颜喀拉地块内部的地壳上地幔电性结构具有成层性，中、强地震的震中多位于壳内低阻层的地层层段，埋深约20 km，这与本次地震的震源深度相似。此外，优势方位角未呈现异常现象的都兰台、白水河台、门源台和金银滩台处于柴达木地块以及祁连地块，而出现异常的玛曲台、甘孜台则与玛多地震共处于巴颜喀拉地块，可见，在玛多地震的孕育过程中，所处的巴颜喀拉地块存在电性结构的改变，而柴达木地块和祁连地块均未记录到明显的电性结构变化。综合分析认为，地电场优势方位角的变化显著受到构造的影响，异常响应多出现在同一地块内部或边缘。
大武台（旧）自2014年改造后，已经积累了近六年的地电场观测资料。通过梳理历史数据可知，大武台（旧）地电场对周边中强震以上的地震反应较为敏感，即在地震前出现明显的地电场优势方位角跳变异常现象。2015年以来，大武台（旧）500 km范围内共记录到MS5.5以上地震共四组，其相关的地电场优势方位角变化详见表2，可见：在四组中强地震前有三组出现了地电场优势方位角异常现象，异常比例为75%，说明大武地震台的地电场优势方位角预报效能较高。异常主要表现为方位角跳变范围出现准同步的大幅度突跳或收缩以及发生偏转等现象，且持续一段时间，在异常时间出现四个月或异常恢复后三个月内，震中附近500 km范围发生M6以上地震（或连续两次M5左右地震）的可能较大，若单台场地岩体出现剪裂（方位角跳变偏转45°或90°），异常信度则更高。
表 2 大武台周围500 km范围内MS5.5地震前地电场优势方位角的变化特征Table 2. Variation characteristic of dominant azimuth based on geoelectric field before MS5.5 earthquakes within 500 km of Dawu station台站 发震日期
年-月-日发震地点 MS 震中距/km 异常 异常出现月份
年-月备注 大武台 2016-01-21 门源 6.4 380 Δα增大 2015-10 多台同步 2016-10-17 杂多 6.2 500 震前2个月快速偏转恢复 2016-08 附近仅一台 2017-08-08 九寨沟 7.0 360 无 无 无 2019-10-28 夏河 5.7 240 快速偏转 2019-06 多台同步 2020-04-01 石渠县 5.6 200 震前3个月快速偏转恢复 2020-01 多台同步 孕震过程中应力不断变化，理论上岩体裂隙结构会因应力的变化而变化。在实际场地中，岩体结构差异会使其裂隙结构对应力变化的响应出现差异，这导致了不同场地裂隙优势方位角α异常具有场地选择性现象。在部分场地，应力积累过程会导致岩体裂隙结构发生剧烈变化，使得该场地的大地电场优势方位角α发生显著变化，例如：α出现大幅度持续突跳、范围收窄、偏转等，而当岩石受压破裂时，剪裂会导致Δα发生约45°的变化，共轭剪裂会导致Δα发生90°左右的变化。此次异常与以往异常特征相似，方位角跳变出现偏转45°或90°，且多台出现准同步异常现象。
4. 结论
通过大武台大地电场优势方位角计算和分析证实玛多MS7.4地震前多个场地出现准同步优势方位角跳变异常是真实存在的，真实反映了地下介质的变化，主要结论如下：
1） 大武台新、旧两套地电场在2020年1月至5月中旬大地电场优势方位角的跳变范围即Δα在10°—20°之内，至2020年5月中旬两套方位角出现或快或慢的持续偏转，其偏转范围最大接近45°或90°，2021年4月方位角跳变出现快速偏转回升过程。随后发生玛多MS7.4地震，两套资料的异常具有准同步性。
2） 对震中500 km范围内的八个地电场的优势方位角的计算表明，处于同一地块的甘孜台、玛曲台与大武台出现准同步异常现象，且异常形态相似，而其它地块的台站均无异常出现，说明地电场优势方位角异常受构造影响显著，异常响应多出现在同一地块内部或边缘。
3） 自2014年大武台地电场观测以来，台站周围500 km范围内共记录到MS5.5以上地震四组。地电场优势方位角计算显示，在四组中强地震前有三组出现了大武台地电场优势方位角异常现象，异常占比为75%，异常主要表现为方位角跳变范围出现准同步的大幅度突跳、或收缩以及发生偏转等现象。地震主要发生在异常出现后四个月或异常恢复后三个月内，台站附近400 km范围存在发生M6以上地震的可能（或连续两次M5地震）。
随着经济的发展和人类活动的加剧，地电场观测环境日益复杂，这使得从台站观测数据中识别相对较弱的地震电信号愈发困难，通过大地电场岩体裂隙水（电荷）渗流（移动）模型分析中强地震孕育前后附近场地岩体裂隙结构的变化特征，为异常提取提供了有效途径。但由于目前地电场场地稀疏且分布不均，导致样本量有限，不能完全真实地反映客观规律。今后随着地电观测的累积、震例的增加，应该加大对该方法的长期深入研究，从而更好地完善预测指标。
-
![]()
图 1 样品Ⅰ(a)和样品Ⅱ(b)的测线布设和加压方向示意图
Figure 1. Arrangement of measuring lines and directions of loading for sample Ⅰ(a) and sample Ⅱ(b)
![]()
图 3 温纳-α装置示意图
A, B为供电电极；M，N为测量电极；I为供电电流；ΔV为测量电极间的电位差；a为电极间距
Figure 3. Sketch diagram of Wenner-α array
I indicates the intensity of current between two current electrodes A and B. ΔVindicates the potential difference between two potential elec-trodes M and N, a indicates electrode spacing
![]()
图 4 样品在无压力状态下电阻率重复测量结果相对第一次测量结果的变化图
Figure 4. The relative resistivity change by comparing repeated measurement with the first measurement before loading
![]()
图 5 样品Ⅰ(a)和样品Ⅱ(b)的应力(单位为MPa)加载和卸载曲线
Figure 5. The stress (in MPa) loading and unloading curves of sample Ⅰ(a) and sample Ⅱ(b)
![]()
图 6 加压前样品Ⅰ和样品Ⅱ中布设的测线L1，L2和L3观测到的视电阻率图像
Figure 6. The apparent resistivity images observed by measuring lines L1, L2, L3 laid in sample Ⅰ and sample Ⅱ before loading
![]()
图 7 加压(左)，卸压(右)过程中不同应力下样品中布设的测线L1(a)，L2(b) L3(c)观测到的视电阻率相对变化图像
Figure 7. Relative resistivity change images observed by measuring lines L1 (a), L2 (b), L3(c) laid in sample Ⅰ. Left column represents loading process, and right column represents unloading process
![]()
图 7 加压(左)，卸压(右)过程中不同应力下样品Ⅰ中布设的测线L1 (a), L2 (b)，L3(c)观测到的视电阻率相对变化图像
Figure 7. Relative resistivity change images observed by measuring lines L1 (a), L2 (b), L3 (c) laid in sample Ⅰ. Left column represents loading process, and right column represents unloading process
![]()
图 8 加压(左)，卸压(右)过程中不同应力下样品Ⅱ中布设的测线L1 (a)，L2 (b)，L3 (c)观测到的视电阻率相对变化图像
Figure 8. Relative resistivity change images observed by measuring lines L1(a), L2(b), L3 (c) laid in sample Ⅱ. Left column represents loading process, and right column represents unloading process
![]()
图 8 加压(左)，卸压(右)过程中不同应力下样品Ⅱ中布设的测线L1(a), L2(b)，L3 (c)观测到的视电阻率相对变化图像
Figure 8. Relative resistivity change images observed by measuring lines L1(a), L2(b), L3(c) laid in sample Ⅱ. Left column represents loading process, and right column represents unloading process
![]()
图 9 温纳-α阵列视电阻率数据对选取示意图
Figure 9. Sketch diagram of apparent resistivity set extraction for Wenner-α array
![]()
图 10 视电阻率数据对对应的有效深度图
Figure 10. Effective depths corresponding to apparent resistivity sets
![]()
图 11 样品Ⅰ在不同有效深度h处的各向异性系数λ*和各向异性主轴方位角α随应力的变化图
空心圆、三角形和实心圆分别表示应力加、卸载过程及加压前(0 MPa)不同应力对应的λ*和α值，下同
Figure 11. The changes of λ* and α with stress for sample Ⅰ at different effective depths h
Open circles, triangles and dots represent values of λ* and α corresponding to different stresses under loading, unloading and before loading, the same below
![]()
图 12 样品Ⅰ在不同深度h处的视电阻率各向异性系数λ*随应力的变化(a)和视电阻率各向异性主轴方位角α的最大值与最小值之差随深度的变化(b)
Figure 12. Variations of λ* with stress curves at different depths of sample Ⅰ (a), variations of αmax-αmin with depth curve of sample Ⅰ (b)
![]()
图 13 样品Ⅱ在不同有效深度h处的各向异性系数λ*和各向异性主轴方位角α随应力的变化
Figure 13. The changes of λ*and α with stress for sample Ⅱ at different effective depths h
表 1 MIR-2007直流电法仪的主要技术指标
Table 1 Main technical indicators of MIR-2007 DC resistivity meter
输入阻抗/MΩ 电压测量范围/V 电流测量范围/A 供电电压/V 工频抑制/dB 自电补偿方式及范围 4104 [-4, 4] [-4, 4] ≤700 ≥80 全量程跟踪式 精度优于±0.5%；分辨率为1 μV 自动补偿 
下载: 导出CSV
-
安金珍, 修济刚, 陈峰, 陈大元. 1996.单轴压力下有补给水岩石电阻率变化各向异性研究[J].中国地震, 12(3): 300-306. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW199408001340.htm An J Z, Xiu J G, Chen F, Chen D Y. 1996. Anisotropy studies of rock resistivity changes under uniaxial pressure and water replenishment[J]. Earthquake Research in China, 12(3): 300-306 (in Chinese). http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW199408001340.htm
安金珍, 周平根, 马麦宁, 陈峰, 董颖, 钟沛林. 2008.用地电阻率探测和监测滑坡体实验[J].地震学报, 30(3): 254-261. http://www.dzxb.org/Magazine/Show?id=26580 An J Z, Zhou P G, Ma M N, Chen F, Dong Y, Zhong P L. 2008. Experiments on exploring and monitoring landslipmass using geoelectric resistivity observations[J]. Acta Seismologica Sinica, 30(3): 254-261 (in Chinese). http://www.dzxb.org/Magazine/Show?id=26580
陈大元, 陈峰, 王丽华. 1983.单轴压力下岩石电阻率的研究:电阻率的各向异性[J].地球物理学报, 26(增刊): 784-792. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW199408001340.htm Chen D Y, Chen F, Wang L H. 1983. Studies on electrical resistivity of rocks: Its anisotropy[J]. Acta Geophysica Sinica, 26(S): 784-792 (in Chinese). http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW199408001340.htm
陈峰, 修济刚, 安金珍, 廖椿庭, 陈大元. 2000.岩石电阻率变化各向异性与微裂隙扩展方位的实验研究[J].地震学报, 22(3): 310-318. http://www.dzxb.org/Magazine/Show?id=26868 Chen F, Xiu J G, An J Z, Liao C T, Chen D Y. 2000. Research on dependence of resistivity changing anisotropy on microcracks extending in rock with experiment[J]. Acta Seismologica Sinica, 22(3): 310-318 (in Chinese). http://www.dzxb.org/Magazine/Show?id=26868
陈峰, 安金珍, 廖椿庭. 2002.弹性约束承载岩石电阻率变化形态研究[J].北京大学学报:自然科学版, 38(3): 427-430. http://www.cnki.com.cn/Article/CJFDTOTAL-BJDZ200203026.htm Chen F, An J Z, Liao C T. 2002. Research on dependence of resistivity changing feature on axial load under elastic constraint with experiment[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 38(3): 427-430 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-BJDZ200203026.htm
陈峰, 安金珍, 廖椿庭. 2003a.原始电阻率各向异性岩石电阻率变化的方向性[J].地球物理学报, 46(2): 271-280. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200302022.htm Chen F, An J Z, Liao C T. 2003a. Directional characteristic of resistivity changes in rock of original resistivity anisotropy[J]. Chinese Journal of Geophysics, 46(2): 271-280 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200302022.htm
陈峰, 廖椿庭, 安金珍. 2003b.剪切和摩擦滑动大模型的视电阻率变化幅度和各向异性[J].地球物理学报, 46(5): 667-675. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200305013.htm Chen F, Liao C T, An J Z. 2003b. Amplitude and anisotropy of apparent resistivity change in big models during shear and frictional slipping[J]. Chinese Journal of Geophysics, 46(5): 667-675 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200305013.htm
陈峰, 马麦宁, 安金珍. 2013.承压介质电阻率变化的方向性与主应力的关系[J].地震学报, 35(1): 84-93. http://www.dzxb.org/Magazine/Show?id=28819 Chen F, Ma M N, An J Z. 2013. Relation between directional characteristics of resistivity changes and principal stress[J]. Acta Seismologica Sinica, 35(1): 84-93 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28819
杜学彬, 薛顺章, 郝臻, 张世中. 2000.地电阻率中短期异常与地震的关系[J].地震学报, 22(4): 368-376. http://www.dzxb.org/Magazine/Show?id=26877 Du X B, Xue S Z, Hao Z, Zhang S Z. 2000. On the relation of moderate-short term anomaly of earth resistivity to earthquake[J]. Acta Seismologica Sinica, 22(4): 368-376 (in Chinese). http://www.dzxb.org/Magazine/Show?id=26877
杜学彬, 阮爱国, 范世宏, 郝臻. 2001.强震近震中区地电阻率变化速率的各向异性[J].地震学报, 23(3): 289-297. http://www.dzxb.org/Magazine/Show?id=26951 Du X B, Ruan A G, Fan S H, Hao Z. 2001. Anisotropy of the variation rate of apparent resistivity near the epicentral region of strong earthquakes[J]. Acta Seismologica Sinica, 23(3): 289-297 (in Chinese). http://www.dzxb.org/Magazine/Show?id=26951
杜学彬, 马占虎, 叶青, 谭大诚, 陈军营. 2006.与强地震有关的视电阻率各向异性变化[J].地球物理学进展, 21(1): 93-100. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200601014.htm Du X B, Ma Z H, Ye Q, Tan D C, Chen J Y. 2006. Anisotropic changes in apparent resistivity associated with strong earthquakes[J]. Progress in Geophysics, 21(1): 93-100 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200601014.htm
杜学彬, 刘君, 崔腾发, 范莹莹, 安张辉, 闫睿, 王丽. 2015.两次近距离大震前成都台视电阻率重现性、相似性和各向异性变化[J].地球物理学报, 58(2): 576-588. doi: 10.6038/cjg20150220 Du X B, Liu J, Cui T F, Fan Y Y, An Z H, Yan R, Wang L. 2015. Repeatability, similarity and anisotropy changes in apparent resistivity recorded by station Chengdu at near distances before two great earthquakes[J]. Chinese Journal of Geophysics, 58(2): 576-588 (in Chinese). doi: 10.6038/cjg20150220
冯锐, 郝锦绮, 周建国. 2001.地震监测中的电阻率层析技术[J].地球物理学报, 44(6): 833-842. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200106012.htm Feng R, Hao J Q, Zhou J G. 2001. Resistivity tomography in earthquake monitoring[J]. Chinese Journal of Geophy-sics, 44(6): 833-842 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200106012.htm
郝锦绮, 冯锐, 周建国, 钱书清, 高金田. 2002.岩石破裂过程中电阻率变化机理的探讨[J].地球物理学报, 45(3): 426-434. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200203013.htm Hao J Q, Feng R, Zhou J G, Qian S Q, Gao J T. 2002. Study on the mechanism of resistivity changes during rock cracking[J]. Chinese Journal of Geophysics, 45(3): 426-434 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200203013.htm
金耀, 张天中, 华正兴, 徐明发, 黄平章. 1983.单轴压缩下多裂隙含水岩样电阻率变化与体积应变[J].地震学报, 5(1): 99-106. http://www.dzxb.org/Magazine/Show?id=28004 Jin Y, Zhang T Z, Hua Z X, Xu M F, Huang P Z. 1983. The electrical resistivity changes and volumetric strain of water-bearing cracked rock samples under uniaxial compression[J]. Acta Seismologica Sinica, 5(1): 99-106 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28004
康云生, 安海静, 马可兴, 谭大诚. 2013.天水地电阻率地表与井下多种观测方式的试验分析[J].地震工程学报, 35(1): 190-195. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201301029.htm Kang Y S, An H J, Ma K X, Tan D C. 2013. Test analysis on geoelectrical resistivity observation combining the surface and deep-well methods at Tianshui seismic station in Gansu Province[J]. China Earthquake Engineering Journal, 35(1): 190-195 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201301029.htm
毛桐恩, 王铁城, 姚家榴, 卢军, 张洪魁. 1995.唐山地震地电阻率各向异性变化特征研究[J].地震学报, 17(4): 505-510. http://www.dzxb.org/Magazine/Show?id=27812 Mao T E, Wang T C, Yao J L, Lu J, Zhang H K. 1995. The characteristics of anisotropy changes in earth resistivity of Tangshan MS7.8[J]. Acta Seismologica Sinica, 17(4): 505-510 (in Chinese). http://www.dzxb.org/Magazine/Show?id=27812
毛桐恩, 胥广银, 范思源, 赵明, 孙景芳. 1999.地电阻率各向异性度的动态演化图象与地震孕育过程[J].地震学报, 21(2): 180-186. http://www.dzxb.org/Magazine/Show?id=27374 Mao T E, Xu G Y, Fan S Y, Zhao M, Sun J F. 1999. Heterogeneity of dynamic evolution pattern of geoelectric resistivity and the seismogenic process[J]. Acta Seismologica Sinica, 21(2): 180-186 (in Chinese). http://www.dzxb.org/Magazine/Show?id=27374
钱家栋, 赵和云. 1988.地电阻率的数值模拟和多极距观测系统[J].地震学报, 10(1): 77-88. http://www.dzxb.org/Magazine/Show?id=28238 Qian J D, Zhao H Y. 1988. Earth resistivity modelling and observational system with multi-separation of electrodes[J]. Acta Seismologica Sinica, 10(1): 77-88 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28238
钱家栋, 马钦忠, 李劭秾. 2013.汶川MS8.0地震前成都台NE测线地电阻率异常的进一步研究[J].地震学报, 35(1): 4-17. http://www.dzxb.org/Magazine/Show?id=28812 Qian J D, Ma Q Z, Li S N. 2013. Further study on the anomalies in apparent resistivity in the NE configuration at Chengdu station associated with Wenchuan MS8.0 earthquake[J]. Acta Seismologica Sinica, 35(1): 4-17 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28812
苏鸾声, 王邦本, 夏良苗, 李验轩. 1982.井下电极观测地电阻率排除地面干扰的实验[J].地震学报, 4(3): 274-276. http://www.dzxb.org/Magazine/Show?id=27987 Su L S, Wang B B, Xia L M, Li Y X. 1982. Elimination of surface disturbances in earth-resistivity measurement by lowering the electrodes in shallow wells[J]. Acta Seismologica Sinica, 4(3): 274-276 (in Chinese). http://www.dzxb.org/Magazine/Show?id=27987
汪雪泉, 郑兆苾, 薛维龙. 2002.安徽嘉山台的多极距电阻率观测[J].地震学刊, 22(1): 1-4. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXK200201000.htm Wang X Q, Zheng Z B, Xue W L. 2002. Multiple polar distance resistivity observation at Jiashan, Anhui seismic station[J]. Journal of Seismology, 22(1): 1-4 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DZXK200201000.htm
王兰炜, 朱旭, 朱涛, 张世中, 刘大鹏, 胡哲, 张宇. 2011.地电阻率多极距观测系统及试验研究[J].地震, 31(1): 20-31. http://youxian.cnki.com.cn/yxdetail.aspx?filename=NMHB201706095&dbname=CJFDPREP Wang L W, Zhu X, Zhu T, Zhang S Z, Liu D P, Hu Z, Zhang Y. 2011. Multi-separation array geo-electrical resistivity observation system and its experimental observation[J]. Earthquake, 31(1): 20-31 (in Chinese). http://youxian.cnki.com.cn/yxdetail.aspx?filename=NMHB201706095&dbname=CJFDPREP
薛顺章, 温新民, 董永德, 梁子彬, 张庆渊, 赵和云. 1994.地电阻率预报地震新方法的研究[J].地震学报, 16(2): 227-234. http://www.dzxb.org/Magazine/Show?id=28607 Xue S Z, Wen X M, Dong Y D, Liang Z B, Zhang Q Y, Zhao H Y. 1994. Study on a new method of earthquake prediction using geoelectric resistivity[J]. Acta Seismologica Sinica, 16(2): 227-234 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28607
张金铸, 陆阳泉. 1983.不同三轴应力条件下岩石电阻率变化的试验研究[J].地震学报, 5(4): 440-445. http://www.dzxb.org/Magazine/Show?id=28041 Zhang J Z, Lu Y Q. 1983. An experimental study on the variation of rock resistivity under triaxially different stresses[J]. Acta Seismologica Sinica, 5(4): 440-445 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28041
张天中, 华正兴, 徐明发. 1985. 1.2千巴围压下岩样破裂和摩擦滑动过程中电阻率变化[J].地震学报, 7(4): 428-433. http://www.dzxb.org/Magazine/Show?id=28132 Zhang T Z, Hua Z X, Xu M F. 1985. Electrical resistivity changes in rock samples during fracture and frictional sliding at 1.2 Kb confine pressure[J]. Acta Seismologica Sinica, 7(4): 428-433 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28132
张学民, 李美, 关华平. 2009.汶川8.0级地震前的地电阻率异常分析[J].地震, 29(1): 108-115. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZN200901014.htm Zhang X M, Li M, Guan H P. 2009. Anomaly analysis of earth resistivity observations before the Wenchuan earthquake[J]. Earthquake, 29(1): 108-115 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DIZN200901014.htm
赵和云, 钱家栋. 1987.郫县台多极距观测资料的反演和分析[J].中国地震, 3(增刊): 73-78. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD1987S1008.htm Zhao H Y, Qian J D. 1987. Inversion and analysis of data of temporal change in apparent resistivity with multi-separations in Pixian stations, Sichuan Province[J]. Earthquake Research in China, 3(S): 73-78 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD1987S1008.htm
周启友, 杭悦宇, 刘汉乐, 戴水汉, 徐建平. 2009.饱水和排水过程中岩石电阻率各向异性特征的电阻率成像法研究[J].地球物理学报, 52(1): 281-288. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200901036.htm Zhou Q Y, Hang Y Y, Liu H L, Dai S H, Xu J P. 2009. The anisotropic properties of rock resistivity during saturation and desaturation processes revealed by electrical resistivity tomography[J]. Chinese Journal of Geophysics, 52(1): 281-288 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200901036.htm
朱涛. 2013.汶川MS8.0地震前区域性地电阻率异常初步研究[J].地震学报, 35(1): 18-25. http://www.dzxb.org/Magazine/Show?id=28813 Zhu T. 2013. Preliminary study on regional geo-resistivity anomaly before the Wenchuan MS8.0 earthquake[J]. Acta Seismologica Sinica, 35(1): 18-25 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28813
Brace W F, Orange A S, Madden T R. 1965. The effect of pressure on the electrical resistivity of water-saturated crystalline rocks[J]. J Geophys Res, 70(22): 5669-5678. doi: 10.1029/JZ070i022p05669
Brace W F, Orange A S. 1966. Electrical resistivity changes in saturated rock under stress[J]. Science, 153(3743): 1525-1526. doi: 10.1126/science.153.3743.1525
Brace W F, Orange A S. 1968a. Electrical resistivity changes in saturated rocks during fracture and frictional sliding[J]. J Geophys Res, 73(4): 1433-1444. doi: 10.1029/JB073i004p01433
Brace W F, Orange A S. 1968b. Further studies of the effects of pressure on electrical resistivity of rocks[J]. J Geophys Res, 73(16): 5407-5420. doi: 10.1029/JB073i016p05407
Diaferia I, Barchi M, Loddo M, Schiavone D, Siniscalchi A. 2006. Detailed imaging of tectonic structures by multiscale earth resistivity tomographies: The Colfiorito normal faults (central Italy)[J]. Geophys Res Lett, 33(9): L09305. doi: 10.1029/2006GL025828.
Jackson D B, Kauahikaua J, Zablocki C J. 1985. Resistivity monitoring of an active volcano using the controlled-source electromagnetic technique: Kilauea, Hawaii[J]. J Geophys Res: Solid Earth, 90(B14): 12545-12555. doi: 10.1029/JB090iB14p12545
Kayal J R, Banerjee B. 1988. Anomalous behaviour of precursor resistivity in Shillong area, NE India[J]. Geophys J Int, 94(1): 97-103. doi: 10.1111/gji.1988.94.issue-1
Kurite K. 1986. How can we identify location of a fracture plane? Anisotropy of electrical conductivity and seismic velocity in dilatancy[J]. Earthq Predict Res, 4(1/2): 39-45.
Parkhomenko E I, Bondarenko A T. 1960. Effect of uniaxial pressure on electrical resistivity of rock[J]. Bull Acad Sci USSR, Geophys Ser, 2: 326.
Rikitake T, Yamazaki Y. 1969. Electrical conductivity of strained rocks, the 5th paper, residual strains associated with large earthquakes as observed by a resistivity variometer[J]. Bull Earthq Res Inst, 47(1): 99-105. https://www.researchgate.net/publication/29772651_6_Electrical_Conductivity_of_Strained_Rocks_The_Fifth_Paper_Residual_Strains_Associated_with_Large_Earthquakes_as_Observed_by_a_Resistivity_Variometer
Stopinski W, Teisseyre R. 1982. Precursory rock resistivity variations related to mining tremors[J]. Acta Geophys Pol, 30(4): 293-320. https://www.researchgate.net/publication/292878216_Precursory_rock_resistivity_variations_related_to_mining_tremors
Utada H. 2003. Interpretation of time changes in the apparent resistivity observed prior to the 1986 eruption of Izu-Oshima volcano, Japan[J]. J Volcanol Geoth Res, 126(1/2): 97-107. http://www.sciencedirect.com/science/article/pii/S0377027303001197
Yamazaki Y. 1965. Electrical conductivity of strained rocks, the first paper, laboratory experiments on sedimentary rocks[J]. Bull Earthq Res Inst, 43(4): 783-802. https://www.researchgate.net/publication/29772438_50_Electrical_Conductivity_of_Strained_Rocks_The_First_Paper_Laboratory_Experiments_on_Sedimentary_Rocks
Zhu T, Zhou J G, Hao J Q. 2012. Experimental studies on the changes in resistivity and its anisotropy using electrical resistivity tomography[J]. Int J Geophys, 2012: 142069. doi: 10.1155/2012/142069.
-
期刊类型引用(9)
1. 张丽琼，高曙德，李娜. 积石山M_S6.2地震前地电场异常特征. 大地测量与地球动力学. 2025(04): 367-372 . 
百度学术
2. 刘伟，柴宁娇，王耀临，李东，冯雪东. 乌加河地电场优势方位角地震前后变化特征. 地震地磁观测与研究. 2025(01): 71-77 . 百度学术
3. 周瀚琳，赵玉红，徐恺晖，李国英. 格尔木地磁观测数据与中强地震关系研究. 高原地震. 2024(01): 39-45 . 百度学术
4. 刘海洋，饶文，徐衍刚，艾萨·伊斯马伊力. 2024年1月23日乌什M_S7.1地震前新疆地电优势方位角变化特征分析. 内陆地震. 2024(02): 182-193 . 百度学术
5. 孙召华，李军辉，李君，孙亮亮，张洋，张钧琪. 2015年安徽阜阳4.3级地震前地电场方位角异常分析. 中国地震. 2024(03): 690-699 . 百度学术
6. 郭雨帆，杜晓辉，董磊，汤兰荣，赵爱平，王甘娇. 基于张衡一号卫星监测的2021年青海玛多7.4级地震前电离层效应. 地震. 2023(02): 85-102 . 百度学术
7. 席继楼，赵家骝，高尚华，王晓蕾，李国佑，孟凡博. 长周期地电场变化特征及机理——以都兰地震台为例. 地震地质. 2023(05): 1092-1111 . 百度学术
8. 赵玉红，李霞，冯丽丽，刘磊，张朋涛，卢嘉沁，孙玺皓. 2次门源地震前地电场优势方位角异常特征研究. 地震地磁观测与研究. 2023(S1): 199-202 . 百度学术
9. 格根，张帆，陈立峰，梁沙沙，王磊. 内蒙古乌加河地电场异常特征研究. 地震地磁观测与研究. 2023(S1): 226-228 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 788
- HTML全文浏览量: 306
- PDF下载量: 32
- 被引次数: 9
