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.
-
引言
许多地震、矿震、火山活动和大型地质构造运动前均可观测到显著的电阻率异常变化(Rikitake，Yamazaki，1969；Stopinski，Teisseyre，1982；Jackson et al，1985；Kayal，Banerjee，1988；毛桐恩等，1999；杜学彬等, 2000, 2001, 2006, 2015；冯锐等，2001；Utada，2003；安金珍等，2008；钱家栋等，2013).为了研究电阻率异常变化的规律和机理并探索其作为地震预报参量的可行性，研究人员在室内和野外进行了大量的电阻率随应力变化的试验，并取得了丰硕的成果(Parkhomenko，Bondarenko, 1960；Brace et al，1965；Yamazaki，1965；Brace，Orange, 1966, 1968a, b；陈大元等，1983；金耀等，1983；张金铸，陆阳泉，1983；张天中等，1985；安金珍等，1996；陈峰等, 2000, 2002, 2003a, b，2013)，对理解地震前、后电阻率的异常变化，探索其形成机理及分析地震监测预报台站记录到的地电阻率数据等起到了重要作用.目前在地震地电阻率观测中，主要使用与室内实验研究装置相同的具有4个固定电极的对称四极测量系统，该系统只能获得地下某个深度到地表的整个深度层的视电阻率值.研究人员主要根据该视电阻率值的大小和变化规律来判断其是否存在异常，进而判断该异常变化是否与地震有关.由于地下电性结构存在非均匀性，可能导致测量系统即使安装在距离震中很近的台站，在某些测向上也难以观测到强震前显著的视电阻率异常变化，例如2008年汶川MS8.0地震前，成都台在NW向就未发现明显的地电阻率异常变化(张学民等，2009).另外，观测到的地电阻率变化中往往会包含由非震因素导致的异常变化，例如年变化和季节性变化，而在实际的电阻率数据处理中，很难将其从观测到的视电阻率变化中完全剔除，因为在剔除这些异常的同时，又会人为地引入新的异常变化，这就使得仅仅依靠测量系统观测到的地电阻率异常变化来判定其与地震活动的相关性的准确度降低、难度增大(朱涛，2013).因此，探索新的地震电阻率观测方式，不仅可以为目前的地震电阻率观测系统提供有益的补充，也对我国的防震减灾具有一定的实际意义.为此一些研究人员进行了新的尝试，例如井下观测(苏鸾声等，1982；康云生等，2013)和多极距观测(赵和云，钱家栋，1987；钱家栋，赵和云，1988；薛顺章等，1994；汪雪泉等，2002；王兰炜等，2011)，目前这些观测系统的效果仍未确定，因此有必要探索其它的地震电阻率监测方法. 冯锐等(2001)利用电阻率层析成像技术于1998年4月14日唐山ML4.4和ML5.0地震前，在河北省昌黎台均观测到了明显的电阻率图像的趋势性变化：震前约5个月的电阻率图像相对于震前9个月的电阻率图像，电阻率出现了区域性降低，而后随着地震的临近，电阻率降低的区域逐渐扩张，大约震后两个月，电阻率图像几乎恢复至震前9个月时的状态.这表明利用电阻率层析成像技术观测地电阻率图像的趋势性变化可以作为目前地震地电阻率观测方法的有益补充.
对于电阻率图像变化特征的实验研究目前还很少. 郝锦绮等(2002)开展了对磁铁矿石英岩样品进行应力反复加载和卸载的单轴压缩实验，研究了电阻率图像的变化特征，并探讨了其变化机理. 周启友等(2009)使用电阻率层析成像方法研究了X, Y, XY及-XY这4个方向上饱水和排水过程中岩石电阻率变化空间分布模式的各向异性特征.这些研究体现了岩石样品在二维和三维结构上的电阻率各向异性变化特征，为利用电阻率层析成像方法进行岩石视电阻率图像随应力变化特征的研究创造了条件.因此，本文拟从以下3个方面进行研究：① 在实验室验证冯锐等(2001)在地震前观测到的电阻率图像的趋势性变化是否可以重现；② 进一步认识和理解视电阻率图像及视电阻率各向异性随应力的变化特征；③ 探讨实验结果在地震、火山喷发以及大型地质构造运动中的意义.
1. 实验方法
本文实验使用两组人工样品，分别为样品Ⅰ和样品Ⅱ各两块，尺寸均为30 cm×30 cm×30 cm.样品制作过程如下：首先，将编号均为42.5R的普通硅酸盐水泥、白水泥、河砂、食盐按照11：2.2：11：1的重量比混合并充分搅拌均匀后浇灌到模具中；然后，在样品尚未固结时在任意一个表面以中心点为基准对称布设3条辐射状测线(图 1a)，其中测线L1垂直于加压方向，测线L2平行于加压方向，测线L3与加压方向夹角为45°.测线L1和L2的电极数目均为43个，隔离系数(深度因子)(Diaferia et al, 2006)均为14，测线L3的电极数目为57个，隔离系数为18.所有测线的电极间距均为5 mm.为了研究高阻体对岩石电阻率图像趋势性变化的影响，在样品Ⅱ中嵌入了一块高阻水泥长方体，尺寸为5 cm×2 cm×2 cm，放置于测线L1中心点右侧正下方2.5 cm深度处(高阻体上表面距测线L1表面2.5 cm)(图 1b)；最后，将样品放置于空气中至少一个月，以充分固结.经测试，高阻体的电阻率为32 Ω·mm，样品的抗压强度均约为44.4 MPa，相对密度均约为1.96 g/cm3.
![]()
图 1 样品Ⅰ(a)和样品Ⅱ(b)的测线布设和加压方向示意图Figure 1. Arrangement of measuring lines and directions of loading for sample Ⅰ(a) and sample Ⅱ(b)实验系统由加压系统和测量系统两部分组成(图 2).加压系统为济南时代试金集团制造的YAW-5000F单轴压力机.视电阻率数据采集采用专门为室内岩石电阻率实验设计的MIR-2007直流电法仪，其技术指标如表 1所示.数据采集过程中，使用了温纳-α观测装置(图 3).
表 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 自动补偿 ![]()
图 3 温纳-α装置示意图A, B为供电电极；M，N为测量电极；I为供电电流；ΔV为测量电极间的电位差；a为电极间距Figure 3. Sketch diagram of Wenner-α arrayI 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实验开始前需要测试仪器在连续工作状态下的稳定性.当MIR-2007直流电法仪处于无压力(0 MPa)状态下连续工作约72小时的过程中，取一块样品Ⅰ重复测量36次.结果表明，任意两次电阻率重复测量的误差均小于5%(图 4)，图中同一颜色的菱形块代表该次测量结果相对于第一次测量结果的相对变化.为了简便，共给出3次测量结果的相对变化，相距第一次测量的时间间隔分别约为12小时(黑色)、48小时(蓝色)和64小时(红色)，测量数据的相对变化基本分布在-5%—5%(两条绿色线)之间.说明在连续工作72小时的情况下，电阻率仪仍是稳定的，测量的数据是可靠的，亦即，当所获得的测量数据的相对变化的绝对值小于5%时，可以认为测量数据是可靠的.
![]()
图 4 样品在无压力状态下电阻率重复测量结果相对第一次测量结果的变化图Figure 4. The relative resistivity change by comparing repeated measurement with the first measurement before loading实验过程中，对两组样品实施单轴应力加载、卸载.样品Ⅰ的应力加、卸载曲线如图 5a所示，应力从0 MPa开始以2.0×10-2 MPa/s的速率均匀地、逐级地加载至11.11，22.22，27.77和33.33 MPa，然后以-2.0×10-2 MPa/s的速率逐级地卸载至0 MPa.样品Ⅱ的应力加、卸载曲线如图 5b所示，应力以2.0×10-2 MPa/s的速率从0 MPa开始逐级加载至11.11，22.22，33.33和40 MPa，然后以-2.0×10-2 MPa/s的速率逐级地卸载至0 MPa.在加载和卸载至不同应力时实施应力保持，直到3条测线的视电阻率数据均采集完毕，同时要求每块样品总的测量时间不超过20小时.
![]()
图 5 样品Ⅰ(a)和样品Ⅱ(b)的应力(单位为MPa)加载和卸载曲线Figure 5. The stress (in MPa) loading and unloading curves of sample Ⅰ(a) and sample Ⅱ(b)2. 实验结果及分析
2.1 视电阻率相对变化图像
前人的研究大多采用两个或者几个电极(Parkhomenko，Bondarenko, 1960；Brace et al，1965；Yamazaki，1965；Brace，Orange, 1966, 1968a, b；陈大元等，1983；金耀等，1983；张金铸，陆阳泉，1983；张天中等，1985；安金珍等，1996；陈峰等, 2000, 2002, 2003a, b，2013)，因此只能建立一点或同一深度上几点的视电阻率随应力变化的关系曲线，而无足够数据来构建视电阻率随应力变化的动态图像.在本研究的实验中，布设了两条均包含43个电极的测线L1和L2，隔离系数均为14，还有一条包含57个电极的测线L3，隔离系数为18.一个隔离系数代表一个深度层，因此测线L1和L2的视电阻率测量数目为287个，测线L3的视电阻率测量数目为513个，足以构建每条测线在不同应力下的视电阻率图像.为了更好地研究视电阻率图像随应力的变化特征，利用式(1) 计算不同应力下的视电阻率相对于加压前(无应力)视电阻率的变化量ΔρNR，即
(1) 式中ρN和ρ0分别为N MPa和0 MPa(加压前)下的视电阻率测量值.利用计算得到的ΔρNR构建视电阻率相对变化(relative resistivity change，简写为RRC)图像.
实验中分别测试了两块样品Ⅰ和两块样品Ⅱ.由于每种样品中的两块样品的RRC图像随应力变化的规律非常相似，文中仅给出每种样品中一块样品的实验结果. 图 6给出了加压前样品Ⅰ和样品Ⅱ各测线对应的视电阻率图像.可以看出，样品的电性结构基本上表现为分层均匀各向同性.在90 mm以浅的区域，样品Ⅰ的视电阻率约为4—54 Ω·mm，样品Ⅱ的视电阻率约为4—70 Ω·mm.样品Ⅰ和样品Ⅱ的视电阻率均随深度增加而增大.另外，样品Ⅱ中的高阻体似乎对其周围的视电阻率影响很小(虚线框标示范围)，可能是由于试验中设置的电极间距为5 mm，远远小于高阻体5 cm的长边长，且测线L1横越高阻体且平行于高阻体的长边所致.
![]()
图 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给出了样品Ⅰ在不同应力下各测线所对应的RRC图像.可以看出，对于垂直于加压方向的测线L1(图 7a)，当应力不超过33.33 MPa(约为破裂应力的75%)时，电阻率降低区域(resistivity-decreased region，简写为RDR)(蓝色至白色)只占图像的极小部分，且随着应力的增加，区域逐渐收缩；电阻率升高区域(resistivity-increased region，简写为RIR)(白色至红色)则占图像的绝大部分，且伴随着RDR收缩的同时，RIR逐渐扩张，这与冯锐等(2001)在两次唐山地震前观测到的电阻率图像的变化趋势非常相似.在应力逐级卸载至0 MPa的过程中，RRC图像中的RIR仍然占据着绝大部分，且继续扩张.对于平行于压力轴的测线L2(图 7b)和与压力轴成45°角的测线L3(图 7c)，当应力小于22.22 MPa(约为破裂应力的50%)时，RDR的分布范围比RIR大很多，但是随着应力增加至33.33 MPa时，RDR迅速收缩而RIR则迅速扩张.在应力逐级卸载至0 MPa的过程中，RDR继续收缩，直至消失，而RIR则占据了RRC图像的绝大部分，且区域内平均电阻率继续增大. 图 8给出了样品Ⅱ在不同应力下各测线所对应的RRC图像.从图中可以看出，RDR和RIR随应力的变化(加载和卸载过程)呈现出与样品Ⅰ相同的变化趋势.但是，当应力较小(如11.11 MPa)时，高阻体对其所在部位及其附近区域电阻率的影响很小，RRC图像随应力的变化趋势与图 7中的结果是一致的.当应力超过22.22 MPa时，横越高阻体且垂直于压力轴的测线L1所对应的RRC图像(图 8a)表明，高阻体所在部位及其附近区域的电阻率增幅相对较大，说明在应力较高时高阻体的存在可能使得RRC图像随应力的变化变得更为突出.
![]()
图 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 process2.2 视电阻率各向异性
计算岩石视电阻率各向异性系数需要测量相交于一点的不同方向测线上的电阻率值.电阻率的测量方法有多种，例如Brace和Orange(1966)在室内利用二极法测量垂直相交的两条测线的视电阻率，计算了横向和纵向电阻率变化的比值，研究了其随应力变化的特征; Kurite(1986)利用不同方向的4条测线的电阻率值，分析了单轴压力下电阻率各向异性与地震波速之间的关系；毛桐恩等(1995)利用地电阻率观测台站在EW和NS两个测向上的电阻率值探讨了地震孕育期的地电阻率各向异性变化特征；安金珍等(1996)和陈峰等(2000，2002，2003a，b)设计了与Kurite(1986)相同的共中心点的4条辐射状测线，讨论了不同加压条件下视电阻率各向异性的变化特征.在上述实验中，每条测线仅布设了4个或8个电极，只研究了一个或两个深度层的电阻率各向异性随应力的变化特征. Zhu等(2012)增加了每条测线的电极数目，研究了多个深度层的视电阻率各向异性随应力和深度的变化特征，但是由于对岩样所加应力相对于破裂应力来说较小，因此并未观测到视电阻率各向异性系数随深度和应力的明显变化.在本实验研究中，参照Zhu等(2012)的方法，在样品的自由表面布设了共中心点的3条测线(图 1)，利用温纳-α观测装置获得了人工岩石样品在单轴压缩至破裂以及卸载的过程中每条测线的视电阻率值，并以此来研究视电阻率各向异性随深度和应力的变化特征.根据温纳-α观测装置的特点，选取了隔离系数为偶数时的视电阻率数据对(图 9).每块样品可提取出6组视电阻率数据对(每个视电阻率数据对由分别来自测线L1，L2和L3的3个视电阻率值组成)，分别对应于6个有效深度(图 10).
![]()
图 9 温纳-α阵列视电阻率数据对选取示意图Figure 9. Sketch diagram of apparent resistivity set extraction for Wenner-α array![]()
图 10 视电阻率数据对对应的有效深度图Figure 10. Effective depths corresponding to apparent resistivity sets视电阻率各向异性系数λ*和裂隙优势方向与加压轴方向之间的夹角(各向异性主轴方位角)α采用陈峰等(2000)的公式进行计算，即
(2) 其中，
(3) 
(4) 
(5) 
(6) 式中，ρs1，ρs2，ρs3分别为沿测线L1，L2和L3测得的视电阻率值.
图 11给出了样品Ⅰ在不同深度处的λ*和α随应力的变化图，可以看出：① 虽然在岩样制作过程中对材料进行了充分的搅拌，但是样品仍然存在原始电性各向异性，且其基本上随深度的加深而减小；② 一般情况下，样品的视电阻率各向异性系数随着应力的增加而减小，即对于原始电性各向异性的样品，随着应力的增加，各向异性程度降低(图 12a)，原始电性各向异性岩石在单轴压缩、剪切、有围压的三轴压缩的过程中也观测到了相同的结果(陈峰等，2003a)；③ 在应力加载过程中，在有效深度为10 mm和20 mm处，样品Ⅰ在不同应力下的最大与最小λ*之差分别为1.26和1.35，而其它有效深度处的最大与最小λ*之差均小于0.55，表明随着应力的增加，裂隙主要在岩样的浅部(20 mm以浅的区域)产生和发展，而在较深的部位(20 mm以深的区域)裂隙产生和发展的速率相对较低.这种现象相应地出现在了卸载阶段，最大和最小λ*之差同样出现在岩样的浅部(有效深度为10 mm和20 mm)，表明样品在加载阶段末期的破裂主要发生在浅部；④ 随着应力的增加，视电阻率各向异性主轴方位角α无明显的趋势性变化，但随着深度的加深，α的最大值与最小值之差明显增大(图 12b)，表明样品浅部裂隙的产生和发展存在明显的方向性(陈峰等, 2000, 2002, 2003a, b)；⑤ 对于加载和卸载阶段，同一深度处、不同应力下的最大与最小λ*差值无固定趋势，但是卸载阶段的最大与最小α的差值几乎在所有有效深度上均明显大于加载阶段的，表明应力对各向异性的控制可能主要体现在视电阻率各向异性主轴方位角α上.
![]()
图 11 样品Ⅰ在不同有效深度h处的各向异性系数λ*和各向异性主轴方位角α随应力的变化图空心圆、三角形和实心圆分别表示应力加、卸载过程及加压前(0 MPa)不同应力对应的λ*和α值，下同Figure 11. The changes of λ* and α with stress for sample Ⅰ at different effective depths hOpen 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给出了不同深度下样品Ⅱ的λ*和α随应力的变化图.可以看出，样品Ⅱ的视电阻率各向异性与样品Ⅰ存在某些相似的变化特征：例如，随着应力的增加，样品Ⅱ的视电阻率各向异性系数λ*总体减小；在加载过程中，不同应力下的最大与最小λ*之差的最大值出现在浅部；卸载阶段的最大与最小α的差值几乎在所有有效深度上明显大于加载阶段的.但是，也存在与样品Ⅰ不同的变化特征，例如：随深度的加深，最大与最小α之差失去了明显的趋势性特征；卸载阶段的最大与最小λ*的差值几乎在所有有效深度(除10 mm深度外)上稍大于加载阶段的.这可能是受样品Ⅱ中高阻体的影响所致.
![]()
图 13 样品Ⅱ在不同有效深度h处的各向异性系数λ*和各向异性主轴方位角α随应力的变化Figure 13. The changes of λ*and α with stress for sample Ⅱ at different effective depths h3. 讨论与结论
前人对电阻率随应力变化的实验研究主要采用4个固定电极的对称四极测量系统，所获得的电阻率值仅能反映岩样电性结构的纵向(一维)变化(Parkhomenko，Bondarenko, 1960；Brace et al，1965；Yamazaki，1965；Brace，Orange, 1966, 1968a, b；陈大元等，1983；金耀等，1983；张金铸，陆阳泉，1983；张天中等，1985；安金珍等，1996；陈峰等, 2000, 2002, 2003a, b，2013).本次实验中，在每条测线上布设了至少43个电极，隔离系数至少为14，在同一应力下可以获得至少几百个视电阻率值，足以构建反映岩样横向和纵向电性结构的视电阻率图像.由此可知，研究重点应集中于视电阻率图像中应力所引起的视电阻率的区域性变化，这亦是前人的实验研究工作中所欠缺的.
本实验使用人工样品的主要原因为：① 干燥的天然岩石通常具有很高的电阻率值，例如花岗岩、辉绿岩及石英和长石的电阻率分别约为1010，108和1012 Ω·m (Parkhomenko，Bondarenko，1960)，导致直流电阻率仪重复测量的精度很差，难以获得可靠的测量数据；② 在天然岩石上布设电极非常困难，容易导致电极因接触电阻太大而无法获得有效、可靠的视电阻率测量数据；③ 本研究中虽然使用的是人工样品，但随着应力的增加，其表现出了与天然磁铁矿样品相似的视电阻率图像变化特征(Zhu et al，2012)，表明人工样品的结果具有较好的外延性，即本文的结果同样适用于天然岩石样品.
本文的实验结果表明：无论是样品Ⅰ，还是样品Ⅱ，所有测线对应的RRC图像随着应力变化均呈现出了相同的变化趋势，即随着应力的增加，RRC图像中的RDR逐渐收缩，而RIR逐渐扩张；另外，高阻体对其所在部位和附近区域的电阻率值增幅有较大影响，对视电阻率图像的趋势性变化没有影响. Zhu等(2012)曾获得了磁铁矿的视电阻率图像随着应力的增加，RDR逐渐扩张，而RIR逐渐收缩的结果.本文的结果看似与Zhu等(2012)的结果存在差异，但实际上二者是一致的，即视电阻率图像随着应力均出现了明显的趋势性变化.另外，本文实验结果与冯锐等(2001)于1998年唐山ML4.4和ML5.0地震前在河北省昌黎台所观测到的特征非常相似，这表明冯锐等(2001)在地震前观测到的电阻率图像的趋势性变化可以在实验室内重现，进而表明电阻率图像的趋势性变化可以与地震、火山活动和大型地质构造运动相关联.
在本次实验中，利用所获得的6个有效深度上的视电阻率数据对，计算并绘制了不同深度下视电阻率各向异性系数λ*和各向异性主轴方位角α随应力的变化图，结果与陈峰等(2003b)一致，即一般情况下，对于原始电性各向异性的样品，随着应力的增加，各向异性程度降低(图 12a，b).同时也获得了新的结果和认识：随着单轴应力的增加，裂隙主要在岩样的浅部产生和发展，而在较深的部位裂隙产生和发展的速率相对较低；在应力加载过程，不同应力下的最大与最小λ*之差最大值出现在浅部；卸载阶段的最大与最小α的差值几乎在所有有效深度上明显大于加载阶段的差值.一些实例研究结果表明，在唐山MS7.8地震前后，震中140 km范围内的地电阻率各向异性变化非常明显(毛桐恩等，1995). 2002年四川雅安陇西河峡谷的坡口滑坡发生前后测得的视电阻率各向异性系数λ*和各向异性主轴方位角α均发生了显著的变化(安金珍等，2008).这表明，本文的结果有助于解释和理解地震和大型地质构造运动引起的视电阻率及其各向异性的变化特征，且该方法可能成为目前地震电阻率观测方法的有益补充.
需要提及的是，本文的研究结果是非常初步的.这是因为在研究中，仅采用了固定速率的单轴加压方式而未考虑更加贴近地震孕育过程的加载方式，例如非固定速率和应力反复加载，而且加压间隔过大(约为破裂应力的25%)，这导致了难以据此更加深入地研究和认识电阻率图像的细微特征，这将在下一步工作中进行改善.
-
![]()
图 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. 刘海洋，艾萨·伊斯马伊力，徐衍刚，韩时雨. 2020年伽师M_S6.4地震前柯坪台地电阻率数据变化分析. 内陆地震. 2024(03): 273-280 . 
百度学术
2. 张宇，柯浩楠，娄晓宇，胡颢迪. 2023年12月18日甘肃积石山6.2级地震前地电阻率变化分析. 地震. 2024(04): 82-96 . 百度学术
3. 林井祥，张继忠，唐志超，张添新. 基于气驱水的煤系岩石电阻率变化实验研究. 煤炭技术. 2022(07): 129-133 . 百度学术
4. 段威，许巍，刘俊博，胡宇博，黄航. 高温高压煤复电阻率各向异性实验研究. 科学技术与工程. 2022(21): 9223-9230 . 百度学术
5. 朱涛. 地震电阻率实验研究新进展及展望. 地球与行星物理论评. 2021(01): 61-75 . 百度学术
6. 解滔，卢军. 含裂隙介质中的视电阻率各向异性变化. 地球物理学报. 2020(04): 1675-1694 . 百度学术
7. 何叶，赵明阶，汪魁，吉跃辉. 基于电阻率试验的真空饱水-干燥循环砂岩性能分析. 重庆交通大学学报(自然科学版). 2020(09): 127-131 . 百度学术
8. 高武平，闫成国，张文朋，王志胜. 电阻率层析成像在沉积区隐伏断层探测中的应用. 物探与化探. 2020(06): 1352-1360 . 百度学术
9. 解滔，卢军. 横向不均匀性对视电阻率各向异性变化的影响和地震前电阻率的变化深度. 地震地质. 2020(05): 1172-1187 . 百度学术
其他类型引用(7)
计量
- 文章访问数: 783
- HTML全文浏览量: 303
- PDF下载量: 32
- 被引次数: 16
