Three measuring lines were arranged on one of free planes of magnetite cuboid samples. Apparent resistivity data were acquired by MIR-2007 resistivity meter when samples were under uniaxial compression of servocontrol YAW-5000F loadingmachine in laboratory. Then we constructed the residual resistivity images using electrical resistivity tomography (ERT) and plotted the diagrams of apparent resistivity anisotropy coefficient (ARAC) λ∗ and the included angle α between the major axis of apparent resistivity anisotropy ellipse and the axis of load with pressure and effective depth. Our results show that with increasing pressure, resistivity and the decreased (D region) and increased (I region) resistivity regions have complex behaviors, but when pressure is higher than a certain value, the average resistivity decrease and the area of D region expand gradually in all time with the increase of pressure, which may be significant to the monitoring and prediction of earthquake, volcanic activities, and large-scale geologic motions. The effects of pressure on λ∗ and α are not very outstanding for dry magnetite samples.
1. Introduction
The striking changes in electrical resistivity and its anisotropy were widely observed before earthquake, mine earthquake, volcanic activities, and geologic motions [1–12]. In order to address the change characteristics of resistivity and its anisotropy, many investigators conducted a plenty of experiments on rock samples under compression and during fracture, shear, and frictional sliding in lab [13–25] and theoretic studies [26–30]. Experimental results reported by Brace et al. [13] and Brace and Orange [14–16] in lab showed that striking changes in electrical resistivity were observed for a variety of igneous and sedimentary crystalline rocks when water-saturated crystalline rocks were stressed to fracture. Resistivity decreased by as much as an order of magnitude. Most of the change occurred at stresses above half the fracture stress. The electrical resistivity of hematite decreased gradually with increasing shock pressure to 440 kb and decreased discontinuously to less than 1 Ω cm from several tens of Ω cm at pressure of 440–520 kb [18], and that of a saturated granite during frictional sliding showed clear variations premonitory to the occurrence of stick-slip: resistivity decreased by a few percent with increasing shear stress; the minimum coincided with a sudden release of shear stress of a few bars, accompanied by an immediate return of resistivity [17]. Morrow and Brace [19] and Yamazaki [31, 32] reported that the electrical resistance of tuffs changed very largely, particularly for small strains and could be applied to earthquake prediction [1]. Studies on the changes in resistivity anisotropy were conducted by Chen et al. [33] firstly in lab. They used the water-saturated crystalline rock samples and arranged three measuring lines parallel to, perpendicular to, and intersecting at the angle of 45° with maximum principle stress axis on one of its free planes, and studied the relations between apparent resistivity anisotropy and pressure. Considering the effects of groundwater, An et al. studied the changes in resistivity and its anisotropy with pressure in the presence of feed water [34]. In order to recognize the expansion of fissure with stress, Chen et al. [22–25] studied the water-saturated granite samples and man-made samples under uniaxial compression, triaxial compression with low confining pressure, and frictional sliding. This previous work played very important role in recognizing and understanding the change behaviors of resistivity and its anisotropy with pressure, but they had to analyze the curves of resistivity and its anisotropy with pressure because of only a few electrodes designed on one measuring line. In order to study the changes in resistivity and its anisotropy with pressure, we use MIR-2007 resistivity meter to acquire apparent resistivity data of dry magnetite samples in lab. In resistivity measurements, 37–120 electrodes and 12–30 “n” factors [35] and Wenner-α array are used, and 210–1205 apparent resistivity measurements are acquired. We construct the resistivity images at each pressure value using ERT and plot the diagrams of λ* and α at each effective depth. Then the changes in resistivity and its image and anisotropy with pressure and depth are studied. We find that when the pressure is higher than a certain pressure, average resistivity will decrease, and the area of D region will expand gradually, and most of λ*and α change slightly with the increase of pressure.
2. Experiments2.1. Samples
Two cuboid magnetite samples used in this study were collected in Yinan, Shandong province, China. Physical properties and geometric dimensions of specimens were listed in Table 1.
Physical properties and geometric dimensions of specimens.
Sample no.
Density (g/cm3)
Dimension (cm)
Modal analysis
YN1
5.04
34.0 × 31.8 × 22.0
64 Magnetite, 29 Quartz, 7 Feldspar
YN2
5.07
34.0 × 29.0 × 26.0
68 Magnetite, 27 Quartz, 5 Feldspar
They mainly consist of magnetite (larger than 60%), which leads to their good conductivity; so it is possible to conduct resistivity measurements easily. Their compression planes are the plane of 31.8×22.0 cm for sample YN1 and of 29.0×26.0 cm for sample YN2 in dimension.
2.2. Experimental System
The experimental system used is shown in Figure 1. The uniaxial servocontrol YAW-5000F loading machine (http://shijin.shuoyi.com/shtml/shijin/product/57fd78bf28f0244d.html) manufactured by Jinnan Shijin Group Co. Ltd. was used for axial load. Three measuring lines are arranged on one of free planes. Line 1 and 2 are perpendicular and parallel to the axis of load, and Line 3 intersects this axis at the angle of 45° (Figure 1). Three measuring lines all intersect at the point O which almost locates the center of the free plane. Copper wires spot welded on the free planes are used as electrodes of resistivity measurements. The parameters of measuring lines of two samples are listed in Table 2. MIR-2007 resistivity meter used in this experiment was specially designed for resistivity measurement of rock sample in lab by Beijing Geological Instrument Factory (http://bjdzyqc.shuoyi.com/shtml/bjdzyqc/index.htm). The input impedance of the meter is about 4104 MΩ, and its accuracy of potential measurement is about 0.15%. Wenner-α array (Figure 2) was used in data acquisition.
Parameters of measuring lines.
ML
NE
SE
NF
Sample YN1
Line 1
78
3 mm
24
Line 2
83
3 mm
26
Line 3
120
3 mm
30
Sample YN2
Line 1
37
5 mm
12
Line 2
52
5 mm
17
Line 3
66
5 mm
21
Note: ML, NE, SE, and NF indicate measuring line, number of electrode, spacing of electrode, and “n” factors, respectively.
Sketch of test system.
Sketch of Wenner-α array. I indicates the intensity of current between two current electrodes A and B. ΔV indicates the potential difference between two potential electrodes M and N.
2.3. Loading Curves and Data Acquisition
The samples were compressed by uniaxial load (Figure 1). Their loading curves are shown in Figure 3. Resistivity measurements were conducted along measuring lines one by one when pressure increased to 0, 0.86 MPa, 2.86 MPa, 7.15 MPa, 14.29 MPa, 21.44 MPa, and 28.29 MPa at uniform rate of 1.43×10-2 MPa/s for sample YN1 and to 0, 2.65 MPa, 6.63 MPa, 13.26 MPa, 19.89 MPa, 33.16 MPa, and 41.11 MPa at 2.65×10-2 MPa/s for sample YN2 and remained constant at each value. It took about 3 hours for sample YN1 and 1.5 hours for sample YN2 for resistivity measurement at each pressure value. In the meantime, time-lapse strain was recorded. To our regret, the strain curves remained straight lines throughout the experiment. We thought that the pressure was not enough high to lead to the measurable deformation recorded by our strain gauge.
The loading curves of samples.
Sample YN1
Sample YN2
3. Results and Analysis3.1. Residual Apparent Resistivity Image (RARI)
Residual apparent resistivity values were calculated viaρNR=ρN-ρ0,
where ρNR and ρN represent the residual and measured apparent resistivity value at the pressure value of N MPa, and ρ0 indicates the apparent resistivity measurement at zero pressure. N takes 0.86 MPa, 2.86 MPa, 7.15 MPa, 14.29 MPa, 21.44 MPa, and 28.29 MPa for sample YN1 and 2.65 MPa, 6.63 MPa, 13.26 MPa, 19.89 MPa, 33.16 MPa, and 41.11 MPa for sample YN2.
As are shown in Figures 4 and 5, the RARIs along three measuring lines have different change behaviors with pressure when the pressure is smaller than 7.15 MPa for sample YN1 (Figure 4) and 6.63 MPa for sample YN2 (Figure 5). Figure 4 indicates that the RARIs could be divided simply into D (blue to green zone) and I (yellow to red zone) region. With the increase of pressure, the average electrical resistivities of D region in Figure 4(a) increase (i.e., the average decreased magnitude becomes smaller) (compared the RARIs at 0.86 MPa and 2.86 MPa) firstly and then decrease (i.e., the average decreased magnitude becomes larger), (compared the RARIs at 2.86 MPa and 7.15 MPa), while those of I region in Figure 4(a) decrease firstly and then increase. The area of D region expands strikingly firstly and then almost remains constant, while that of I region reduces outstandingly firstly and then almost keeps constant. The electrical resistivities of D region in Figure 4(b) decrease in all time, and those of I region decrease firstly and then increase. The area of D region expands firstly and then reduces. The electrical resistivities of D and I region in Figure 4(c) both increase firstly and then decrease. The area of D region reduces firstly and then almost remains constant while that of I region increases firstly and then almost remains constant. With the increase of pressure from 2.65 MPa to 6.63 MPa, the electrical resistivities of D and I regions in Figures 5(a) and 5(b) all increase slightly. The area of D region expands slightly, and that of I region reduces slightly in Figure 5(a) while it is going in the opposite direction in Figure 5(b). The electrical resistivities of D and I region in Figure 5(c) both decrease slightly. The area of D region expands slightly while that of I region reduces slightly. When the pressure is larger than 7.15 MPa for sample YN1 (Figure 4) and 6.63 MPa for sample YN2 (Figure 5), the electrical resistivities of D regions in Figures 4 and 5 all decrease gradually while those of I regions probably have complex change behaviors with the increase of pressure. For instance, they decrease firstly (As shown in the RARIs at 6.63 MPa and 13.26 MPa in Figure 5(a)), then increase (As shown in the RARIs at 13.26 MPa and 19.89 MPa in Figure 5(a)) and decrease gradually, at last (As shown in the RARIs at 19.89 MPa, 33.16 MPa, and 41.11 MPa in Figure 5(a)). However, the areas of D regions in Figures 4 and 5 all expand gradually and those of I regions are in opposite direction.
The residual resistivity images of sample YN1. (a) Line 1, (b) Line 2, and (c) Line 3.
The residual resistivity images of sample YN2. (a) Line 1, (b) Line 2, and (c) Line 3.
Our results show, as the previous investigators [14–19] reported, that electrical resistivities of rocks have complex behaviors with the increase of pressure. However, due to only two electrodes used in (volume) resistivity measurements, insufficient resistivity measurements could be used to construct resistivity image, and only the curves of resistivity versus pressure could be used to study the changes in resistivity with pressure in these previous studies. In this paper, 37–120 electrodes and 12–30 “n” factors are used, and then 210–2205 resistivity measurements are acquired, which allow us to construct resistivity images easily we have the chance to study the changes in resistivity and its image with pressure in laboratory. Our results indicate that when pressure increases to a certain value, with the increase of pressure, the electrical resistivities will decrease, and the area of D region will expand gradually in all time, which may be significant to the monitoring and prediction of earthquake, volcanic activities, and large-scale geologic motions.
3.2. Apparent Resistivity Anisotropy3.2.1. Data Extraction
Studies on resistivity anisotropy require that resistivity measurements are conducted along at least two measuring lines arranged in different directions and that they must intersect at a point. For example, Mao et al. [7, 8] used the resistivity data measured along E-W and N-S trending measuring lines in studying the changes in the degree of earth resistivity anisotropy during the course of earthquake preparation. Chen et al. [22–25] designed four measuring lines in studying the changes of apparent resistivity anisotropy in lab. In this paper, we arranged three measuring lines intersecting at point O (Figure 1), which allowed us to study the changes in apparent resistivity anisotropy at point O with pressure and depth. According to the used minimum “n” factor and the characteristics of Wenner-α array in resistivity measurements, we could extract 12 apparent resistivity sets (Figure 6 illustrates how to extract the apparent resistivity set.) which are corresponding to 12 effective depths (Figure 7(a)) for sample YN1 and 6 apparent resistivity sets which are corresponding to 6 effective depths (Figure 7(b)) for sample YN2. It must be noted that one apparent resistivity set contains three apparent resistivity values.
The illustration of extracting apparent resistivity set for Wenner-α array. In the studies on apparent resistivity anisotropy at point O, we need the apparent resistivity data arranged along the vertical dash line. According to the characteristics of Wenner-α array, the apparent resistivity measurements could be extracted every one “n” factor. That is, the apparent resistivity measurements just corresponding to odd “n” factors and on the vertical dash line will be extracted out. Black solid circles represent the extracted data.
The diagram of apparent resistivity sets corresponding to effective depths. Because the effective depth of Wenner-α array was evaluated to be about 0.5 time of electrode spacing [36, 37], the 12 apparent resistivity sets were corresponding to the effective depths of about 1.5 mm, 4.5 mm, 7.5 mm, 10.5 mm, 13.5 mm, 16.5 mm, 19.5 mm, 22.5 mm, 25.5 mm, 28.5 mm, 31.5 mm, and 34.5 mm for sample YN1, and 6 ones were corresponding to the effective depths of about 2.5 mm, 7.5 mm, 12.5 mm, 17.5 mm, 22.5 mm, and 27.5 mm for sample YN2 according to Figure 6.
Chen et al. [25, 33] presented the formula to calculate ARAC λ*λ*2=2m2+Q2l-m2+Q2,l=ρs1-2+ρs2-2,m=ρs1-2-ρs2-2,Q=2ρs3-2-l,
where ρs1, ρs2, and ρs3 represent the apparent resistivities measured along three measuring lines perpendicular to, parallel to, and intersecting at the angle of 45° with the axis of load.
Seen from the curves of λ* and α with pressure and depth for sample YN1 (Figure 8), the absolute differences between the minimum and the maximum λ* at 10.5 mm (Figure 8(d)), 13.5 mm (Figure 8(e)), and 19.5 mm (Figure 8(g)) effective depths are about 1.6, 2.0, and 2.0, while those at the rest depths are all smaller than 0.5. The maximum change of α, only about 7°, occurs at 4.5 mm depth. For sample YN2, the absolute difference between the minimum and the maximum λ* at 17.5 mm (Figure 9(d)) is about 4.0, while those at the rest depths are all smaller than 1.0. The maximum change of α, about 80°, occurs at 22.5 mm depth, but those at the rest depths are all smaller than 5°. These results indicate that the effects of pressure on λ* and α are not very striking for dry magnetite samples.
The diagrams of λ* and α with pressure for sample YN1. (a)–(l) represent the diagrams at the effective depths of 1.5 mm, 4.5 mm, 7.5 mm, 10.5 mm, 13.5 mm, 16.5 mm, 19.5 mm, 22.5 mm, 25.5 mm, 28.5 mm, 31.5 mm, and 34.5 mm.
The diagrams of λ* and α with pressure for sample YN2. (a)–(f) represent the diagrams at the effective depths of 2.5 mm, 7.5 mm, 12.5 mm, 17.5 mm, 22.5 mm, and 27.5 mm effective depths.
4. Discussions and Conclusions
In the previous studies, the volume resistivity changes with pressure were actually studied [14–19]. It was impossible to construct resistivity image because of insufficient data; the curves of resistivity versus pressure were analyzed. In the present paper, resistivity measurements are conducted in the case of designed 37–120 electrodes and 12–30 “n” factors, and 210–2205 resistivity values acquired are sufficient to allow us to use ERT to construct resistivity images which may contain much aplenty information than curves. Thus we can study the changes in not only resistivity but also in the areas of D and I region with pressure. Our results indicate that resistivity and the areas of D and I region have complex behaviors with increasing pressure when pressure is lower than a certain value such as 7.15 MPa for sample YN1 and 6.63 MPa for sample YN2 in this paper, but when pressure is higher than this value, their behaviors become very simple with increasing pressure. For instance, when pressure is higher than 7.15 MPa for sample YN1 and 6.63 MPa for sample YN2, the average resistivity of D region will decrease, and its area will expand gradually in all time with the increase of pressure. These changes were ever observed in situ by Feng et al. [9] prior to Tangshan ML 5.0 and 4.4, which probably indicate that our experimental results may be helpful and significant to the monitoring and prediction of earthquake, even for volcanic activities and large-scale geologic motions.
In the previous studies, the changes in apparent resistivity anisotropy at one or two depths with pressure were studied [22–25, 33, 34] because of at most eight electrodes and two “n” factors used in resistivity measurements, which lead to the difficulty in the studies on its changes with depth. In this paper, as described above, at least 78 electrodes and 24 “n” factors for sample YN1 and 37 electrodes and 12 “n” factors for sample YN2 (Table 2) and Wenner-α array were used in resistivity measurements, which allow us to calculate λ* and α at 6 or 12 effective depths (Figure 7) and give us the chance to study their changes with pressure at more depths. Our results indicate that λ* and α do not change outstandingly with the increase of pressure at a certain depth. However, the striking changes in apparent resistivity anisotropy were observed before earthquake [7, 8, 38]. We thought that we did not use water-saturated samples while the underground rocks are often in saturation because of the presence of rich groundwater. Therefore, water-saturated samples would be considered in our future study.
Acknowledgments
The authors thank two anonymous reviewers and the editor Sabatino Piscitelli for their constructive suggestions for this paper. This research was supported by the Special Project for the Fundamental R&D of the Institute (DQJB07B09).
YamazakiY.Electrical conductivity of strained rocks, 5, residual strains associated with large earthquakes as observed by a resistivity variometer19694799FittermanD.MaddenT.Resistivity observations during creep events at melendy Ranch, California1977823354015408StopinskiW.TeisseyreR.Precursory rock resistivity variations related to mining tremors19823042933202-s2.0-0020386696JacksonD. B.KauahikauaJ.ZablockiC. J.Resistivity monitoring of an active volcano using the controlled- source electromagnetic technique: Kilauea, Hawaii19859014125552-s2.0-0022208660KayalJ. R.BanerjeeB.Anomalous behaviour of precursor resistivity in Shillong area, Northeast India198894197103YukutakeT.YoshinoT.UtadaH.WatanabeH.HamanoY.ShimomuraT.Changes in the electrical resistivity of the central cone, Miharayama, of Oshima Volcano observed by a direct current method19904231511692-s2.0-0025683507MaoT.WangT.YaoJ.LuJ.ZhangH.The variations of the degree of ground resistivity anisotropy during the tangshan earthquake1995846216272-s2.0-5124916583210.1007/BF02650530MaoT.XuG.FanS.ZhaoM.SunJ.Dynamic evolution patterns of the degree of ground resistivity anisotropy and the seismogenic process19991221992062-s2.0-0005308699FengR.HaoJ.ZhouJ.Resistivity tomography in earthquake monitoring20014468198302-s2.0-0035728715FengZ.YangJ.MeiW.GengJ.WangX.LiuY.Preliminary study on characteristics of seismic precursor of earth resistivity anisotropy degree in east China20041722442502-s2.0-5044220814UtadaH.Interpretation of time changes in the apparent resistivity observed prior to the 1986 eruption of Izu-Oshima volcano, Japan20031261-2971072-s2.0-004162932010.1016/S0377-0273(03)00119-7AnJ.ZhouP.MaM.ChenF.DongY.ZhongP.Experiments on exploring and monitoring landslip-mass using geoelectric resistivity observations20082132582662-s2.0-6274918092310.1007/s11589-008-0258-5BraceW. F.OrangeA. S.MaddenT. M.The effect of pressure on the electrical resistivity of water—saturated crystalline rocks1965702256695678BraceW. F.OrangeA. S.Electrical resistivity changes in saturated rock under stress19661533743152515262-s2.0-24944585311BraceW. F.OrangeA. S.Electrical resistivity changes in saturated rocks during fracture and frictional sliding196873414331444BraceW. F.OrangeA. S.Further studies of the effects of pressure on electrical resistivity of rocks1968731654075420WangC. Y.GoodmanR. E.SundaramP. N.MorrisonH. F.Electrical resistivity of granite in frictional sliding: application to earthquake prediction1975212525528KondoK.MashimoT.SawaokaA.Electrical resistivity and phase transformation of hematite under shock compression19808529779822-s2.0-0018911566MorrowC.BraceW. F.Electrical resistivity changes in tuffs due to stress198186429292934PearsonC.MurphyJ.HermesR.Acoustic and resistivity measurements on rock samples containing tetrahydrofuran hydrates: laboratory analogues to natural gas hydrate deposits198691141413214138SkagiusK.NeretnieksI.Diffusivity measurements and electrical resistivity measurements in rock samples under mechanical stress19862245705802-s2.0-0022696609ChenF.XiuJ.AnJ.LiaoC.ChenD.Detecting rupture precursors and determining the main fracture spread direction of rock with dynamic rock resistivity change anisotropy20001322342372-s2.0-15944382017ChenF.XiuJ.AnJ.LiaoC.ChenD.Research on dependence of resistivity changing anisotropy on microcracks extending in rock with experiment20001333313412-s2.0-0348007558ChenF.LiaoC.AnJ.Directional characteristic of resistivity changes in rock of original resistivity anisotropy20034623813952-s2.0-0041940315ChenF.LiaoC.AnJ.Amplitude and anisotropy of apparent resistivity change in big models during shear and frictional slipping20034659549652-s2.0-0346947073FittermanD.Theoretical resistivity variations along stressed strike—slip faults19768126490949152-s2.0-0016992885JohnsonD. L.ManningH. J.Theory of pressure dependent resistivity in crystalline rocks198691111161111617TeisseyreK.Anisotropy of electric resistivity related to crack process before fracturing1989372185192TeisseyreK.Simulation of anisotropic changes of resistivity in primarily anisotropic rocks under load1991413316333TeisseyreK.Modeling of crack induced resistivity changes-applications to earthquake studies1997242741632-s2.0-0031402581YamazakiY.Electrical conductivity of strained rocks, 1, Laboratory experiments on sedimentary rocks196543783YamazakiY.Electrical conductivity of strained rocks, 2, Further experiments on sedimentary rocks1966441553ChenD.ChenF.WangL.Studies on resistivity of rock samples under uniaxial pressure-anisotropy of resistivity198326784792AnJ.XiuJ.ChenF.ChenD.Anisotropy studies of rock resistivity changes under uniaxial pressure and water replenishment1996123300306DiaferiaI.BarchiM.LoddoM.SchiavoneD.SiniscalchiA.Detailed imaging of tectonic structures by multiscale Earth resistivity tomographies: the Colfiorito normal faults (central Italy)20063392-s2.0-3374568512610.1029/2006GL025828L09305RoyA.ApparaoA.Depth of investigation in direct current methods1971365943957EdwardsL. S.A modified pseudosection for resistivity and IP1977425102010362-s2.0-0017525952QianF.ZhaoY.HuangY.Calculation of the parameters of georesistivity anisotropy and case history of earthquake precursors1999946176272-s2.0-17944380361