原始电阻率各向异性岩石电阻率变化的方向性

在3种非含岩成分造成的岩石原始电阻率各向异性标本(在采集和加工过程中，标本完整性遭 到破坏的标本)，和一种由含岩成分形成的 岩石原始电阻率各向异性标本上，布设多个电极，并将标本用水饱和. 然后用多电极组合法 ，将己布电极组合成不同方向、不同极距的电阻率各向异性测线. 采用单轴压缩、低围压 三 轴压缩和剪切3种加载方式，对标本进行动态电阻率变化实验. 观测标本电阻率随承载力的 变化. 实验结果表明∶对于岩石原始电阻率各向异性标本，在电阻率变化各向异性方面，与 岩石 原始电阻率各向同性标本相似，即∶裂隙和破碎带通过区域的测点，视电阻率变化各向异性 结果好，４种组合求得的４个各向异性主轴方向趋向一致，且与破碎带方向基本吻合；裂隙 和破碎带不经过区域的测点，４个视电阻率变化各向异性主轴方向不一致，或者根本求不出 各向异性解. 这后一种情况，在裂隙面平行测量面时，表现最为明显.     

Anisotropy of the apparent resistivity variation rate near the epicentral region for strong earthquakes

IntroductionThe geo-electrical anisotropy before earthquake is of special significance for the prediction of earthquake location and for the research on crustal stress status near the epicentral region. The anisotropy has been studied from various aspects (CHEN, et al, 1983; QIAN, et al, 1985, 1996; LU, et al, 1990a; MAO, et al, 1995; RUAN, et al, 1999). But the anisotropic phenomena presented by geo-electrical data are not practically satisfactory and not clear as compared with the res…     

强震近震中区地电阻率变化速率的各向异性

用归一化月速率方法处理了６次强震前近震中区11个台的地电阻率数据，得到在孕震中短期至短临阶段，与主压应力方向正交（或近于正交）测向的地电阻率下降变化的速率大于平行（或近于平行）主应力方向的变化速率．产生变化速率各向异性的原因可能是:在孕震后期的扩容阶段，裂隙走向沿主压应力方向优势排列，导电流体活动产生真电阻率变化速率的各向异性，表现为地表视电阻率变化速率的各向异性．本文为视电阻率各向异性的实验结果提供了震例支持，物理解释比较清晰，可能成为研究某些强震孕育晚期震源区及其附近地壳应力状态的参考依据．      

APE理论与地电阻率前兆

介绍了含裂隙介质的ＡＰＥ理论（ａｎｉｓｏｔｒｏｐｉｃ ｐｏｒｏ－ｅｌａｓｔｉｃｉｔｙ）。根据ＡＰＥ理论,对裂隙膨胀提出了新的看法,认为当岩石不产生大量破裂时,裂隙的扩张和闭合是互为依存、彼此平衡的。用笔者建立的裂隙方面性集中体积模型的各向异性电导率张量表达式,研究了应力场导致裂隙纵横比变化与地壳岩石电阻率变化的关系。结果表明,即使地壳岩石的平均体应变很小或不产生新的裂隙,也可以出现电阻率大的突变。     

剪切和摩擦滑动大模型的视电阻率变化幅度和各向异性

用石英砂、河砂和水泥模压制成1m×1m×0.3m的均匀和非均匀介质模型，对模型进行剪切和摩擦滑动实验，观测模型不同位置和不同方位的视电阻率变化幅度和各向异性特征. 均匀介质模型实验结果为：电阻率变化与测线的位置有关，距裂隙近的测线，电阻率变化幅度大，从百分之几到百分之几十；距裂隙远的测线，电阻率变化幅度小，从百分之几到不变. 对于非均匀介质模型，电阻率变化幅度要小一些，最大也只有百分之几. 电阻率变化除与测线的位置有关外，还与测线的方位有关，同一测点不同方位的测线，电阻率变化幅度不同，有的差别很大. 剪切与摩擦滑动两种加载方式的电阻率变化幅度数量级相同. 裂隙穿过部位及其附近测点的电阻率变化各向异性主轴方位解与剪切和摩擦滑动的实际裂隙方位吻合较好.     

承压介质电阻率变化的方向性 与主应力的关系

介质电阻率变化的方向性与主应力方向的关系研究一直受到岩石电阻率实验和地震学等学科的关注. 本文在电极面4 cm×8 cm, 6 cm×12 cm, 8 cm×16 cm和4.6 cm×10 cm的10块分别为花岗岩、 凝灰质粗砂岩(夹砾石)和细晶花岗岩标本上及100 cm×100 cm的1块人工合成模型上, 采用对称四极法布设多位置、 多极距、 多方位电阻率测线, 用单轴压缩、 二维约束差应力压缩、 低围压三轴差应力压缩和真三轴差应力压缩方式加载. 其中在单轴压缩和低围压三轴压缩加载方式下, 7块岩石标本中的4块被压破, 3块没有被压破, 只压到电阻率出现明显的下降异常时卸压. 另在上述10块岩石标本中, 还有3块是电阻率原始各向异性标本. 将这样11块标本所获得的多位置、 多极距、 多方位测线电阻率资料, 按从开始加载到破裂的过程, 分成30%附近破裂应力段、 50%附近破裂应力段、 80%附近破裂应力段和100%破裂应力段, 研究同一测点多方位测道中该4个破裂应力段视电阻率变化最大(最小)幅度测线方位与力源最大主压应力方位的关系, 结果未发现二者之间有确定的关系.      

地震地电阻率各向异性变化数值模拟

将APE模式中差应力与裂隙的演化理论与前期研究中建立的张量形式的各向异性电性模型相结合,推导出了裂隙演化参数(密度、孔压、纵横比、闭合角)与电导率关系的动态方程,用数值模拟方法研究了饱水岩石在水平向应力作用下电阻率的各向异性变化．并结合实验和观钡5结果,讨论了电阻率各向异性变化的可能机理。     

含裂隙介质中的视电阻率各向异性变化

我国50多年的视电阻率连续观测结果表明,大地震前近震中区域的视电阻率呈现出与主压应力方位有关的各向异性变化,即:垂直于主压应力方向观测的变化幅度最大,平行方向最小或不明显,斜交方向介于二者之间.目前我国定点台站视电阻率观测的探测范围主要在浅层沉积层以内,通常含有较多的含水裂隙.本文将地下岩土介质简化为由固体基质和含流体/气体裂隙组成的固液气三相介质,且基质、流体和气体具有标量形式的电阻率,推导出了包含基质和流体电阻率、裂隙率、饱和度和裂隙面积率因子的电阻率张量表达式.以裂隙的扩展/闭合表示应力作用下裂隙的变化,得到了电阻率随裂隙变化的微分形式,电阻率变化对裂隙体积变化放大系数的表达式和裂隙横向变化对纵向电阻率影响的横向权系数的表达式.在此基础上得到了介质电阻率和视电阻率的各向异性变化特征:对于含水裂隙介质,无论裂隙如何变化,均是最小主轴方向电阻率的变化幅度大于其他方向;对于含水孔隙介质,沿孔隙主要变化方向的主轴电阻率变化幅度大于其他方向.对于各向异性变化,视电阻率和介质电阻率存在π/2的方向差异.相较于含水岩石,无水岩石介质电阻率的各向异性变化不显著.本文提出的电阻率表达式可以对实验室和野外实际观测的许多结果做出合理的解释.     

2013年芦山MS7.0地震前甘孜台地电阻率变化分析

甘孜台地电阻率N30°E测道观测资料自2011年7月开始出现趋势上升变化,N60°W测道地电阻率则从2012年出现趋势上升变化.2011年测区原317国道实施扩建工程,N30°E测道测量电极分别向供电极方向移动10m.采用甘孜台电测深曲线以水平层状模型反演了测区的电性结构,理论计算表明,测量电极的移动将会引起N30°E测道4Ω·m的上升变化,扣除这部分变化后,甘孜台两测道观测值于2012年同步上升.以水平层状模型计算了甘孜台两测道各层介质的影响系数,两测道浅层两层介质影响系数均为负,能合理地解释甘孜台地电阻率在雨季降水量增加时观测值上升、旱季降水量减少时观测值下降这一“夏高冬低”的年变现象.建立三维有限元模型计算了317国道拓宽部分对观测的影响,计算结果表明,拓宽部分仅能引起N60°W测道约0.15Ω·m的下降变化和N30°E测道约0.1Ω ·m的上升变化,其对观测的影响非常小.同时2013年1月甘孜台两测道年变低值显著高于2008年以来各年的年变低值,在芦山地震前呈同步的上升变化,但是与汶川地震前的下降变化相反,因此甘孜台自2012年的趋势上升变化是不是芦山地震的前兆异常还难以确定.     

二维大地电磁各向异性参数对视电阻率的影响研究

大地电磁法是广泛应用于深部地质结构探测、油气和矿产资源勘查等领域的一种地球物理方法.电性各向异性对电磁观测数据有很大影响,但介质各向异性参数对不同模式视电阻率的影响还较少有较为系统的研究.本文基于Maxwell方程组,推导了二维大地电磁场在任意各向异性介质中电场和磁场相互耦合的变分方程,结合有限单元法及并行计算编写了二维大地电磁任意各向异性正演程序,采用三角形网格剖分.验证程序正确性后,以倾斜板状体作为模型来研究三个主轴电阻率及三个旋转欧拉角和四种模式的视电阻率之间的关系.结果表明,主轴各向异性时,xy模式视电阻率几乎只受x方向电阻率影响,yx模式视电阻率主要受y方向电阻率影响,但同时也受z方向电阻率一定影响;三个欧拉角中只有倾角不为零时,yx模式视电阻率受倾角大小的影响较大,xy模式视电阻率几乎不受倾角的影响;只有走向角不为零时,四种模式的视电阻率同时受x、y两个主轴电阻率和走向角的大小的影响.     

THE EFFECTS OF ANISOTROPY ON SQUARE ARRAY RESISTIVITY MEASUREMENTS*

An account is given of the use of the square array technique in investigating the surface effects of rotational anisotropy when the axes of anisotropy are inclined to the surface. It is shown that, as with other arrays, two anisotropy parameters and n can be derived by varying the array orientation. On the basis of these considerations, the effects of such anisotropy on the values of the mean apparent resistivity and azimuthal inhomogeneity ratio normally obtained in square array measurements is reviewed. Particular attention is paid to the variation of resistivity with orientation and it is noted that, in areas of moderate anisotropy, this variation is lower for the square than for the Schlumberger array. In addition to this advantage, the azimuthal inhomogeneity ratio obtained from square array measurements may be used to indicate the severity of anisotropy in an area and two field examples of this use are given. Where anisotropy is severe, gross variations of apparent resistivity with orientation are obtained with either square or collinear arrays. In these circumstances, the use of crossed measurements is considered and the particular stability of the crossed square array demonstrated.     

九寨沟M_S7.0地震地电阻率变化时空演化分析

本文应用归一化变化速率(NVRM)方法,研究了2017年九寨沟M_S7.0地震周边800km范围内25个地电阻率台站2012—2018年的连续观测数据,分析了震中区周围的地电阻率时、空演化过程.结果显示:(1)25个台站中共有14个台站在震前出现了地电阻率变化,结合本次地震震源机制等分析其时空变化,认为邻近地震破裂带的临夏、通渭、武都和宝鸡台地电阻率在震前出现了下降—折返回升型的变化,符合前人已总结的强地震前地电阻率的变化过程.(2)地震发生前震中区周围大区域内出现了地电阻率负异常变化,且以震中区为中心呈条带分布,异常区的长轴方向与地震断裂带走向或地震烈度分布的长轴方向几乎垂直,与地震主压应力轴方向吻合.本文认为九寨沟M_S7.0强震发生之前,震中区周围出现的区域性地电阻率异常空间丛集现象是很好的中期和短临前兆;地电阻率震前异常动态演化、各向异性等特征对认识本次强震发震断层活动以及震源区应力场分布有启示意义.     

裂隙非均匀分布岩石中的电位及视电阻率计算

计算了在均匀高阻介质中含有低阻裂隙的岩石的电场分布，并由此推导出裂隙非均匀分布的岩石中视电阻率的计算公式．结果表明，裂隙的埋深、规模、走向（平行或垂直视电阻率测线）以及真电阻率等，对视电阻率的观测和计算结果都有不同程度的影响．本文结果对视电阻率法用于地震分析预报提供了一种有效的计算方法．      

岩石电阻率的实验研究

本文首先概述岩石在室内各种受载条件下电阻率的实验研究工作,这对于认识和理解电阻率变化特征、规律以及探索其机理起着重要的作用。但是在实验中主要采用4个固定电极的测量系统,它只能获得从某个深度到岩样表面的整个深度层的视电阻率变化的综合信息,而无法获得不同深度层的视电阻率变化特征,不利于进一步研究电阻率及其各向异性变化的机理。我们利用电阻率层析成像方法获得了多个方向的电阻率数据,据此可以获得随深度变化的电阻率及其各向异性曲线以及随应力变化的电阻率图像,可能为深入理解和研究电阻率变化机理起到积极的帮助。      

基于各向异性介质模型讨论地电阻率各向异性变化

讨论了均匀方位各向异性介质模型,运用电阻率变化的视各向异性与真各向异性的理论关系。研究了20个强震近震中区31个地电台的视电阻率变化各向异性与主应力的关系发现:与主压应力方位正交的测向(包括近于正交的测向)月速率变化幅度大于平行主压应力方位测向(包括近于平行的测向)的变化幅度;视各向异性变化与地震震级、震中距等均有一定的关系。     

Experiments on exploring and monitoring landslip-mass using geoelectric resistivity observations

In the experiments, a high-density resistivity method is used to explore the electric structure of landslip mass, and a resistivity-changing anisotropy method is used to monitor the orientation and speed of main fracture extending of landslip mass. The results are as follows. 1 The exploring experiments have verified a part of creep deformation borderline, the depth and thickness of groundwater horizon, and the property of superstrata in the landslip mass investigated formerly, which have proved that the landslip belts contain rich groundwater; 2 The main fracture extending orientation inferred from the resistivity-changing anisotropy accords with the strike of fracture belt of landslip mass deduced from GPS displacement. Moreover, the changing rates of resistivity-changing anisotropy coefficient matches with the changing speeds of deep displacement of landslip mass were measured by suing clinometer in the borehole.     

APPARENT RESISTIVITY TEMPORAL VARIATION CHARAC-TERISTICS AFFECTED BY THE FIXED DISTURBANCE SOURCE ON SURFACE OF MEASURING AREA

Current leakage,metallic conductor,and local anomalous resistivity body are main disturbance sources which affect the successive observation of apparent resistivity in stations,besides the observing system failure.We construct a finite element model using a 3-layered horizontal medium to discuss the dynamic characteristics of disturbances caused by metal conductor and local anomalous resistivity body in the measuring filed.The numerical results show that low resistivity source which is located in areas where the sensitivity coefficient is positive will cause decline on apparent resistivity observation.While low resistivity source will cause increase when it is located in areas where the sensitivity coefficient is negative.Disturbance caused by high resistivity source is opposite to the one from low resistivity source.The general dynamic feature of disturbance is that the disturbance amplitude increases as the resistivity of shallow layer decreases,while the amplitude declines when the shallow layer's resistivity increases.For the measuring direction which has normal annual variation form,low resistivity source which is located in area where the sensitivity coefficient is positive will increase the annual variation amplitude,while it will reduce annual amplitude when it is in a negative sensitivity coefficient area.Annual amplitude changes caused by high resistivity source are opposite to the changes caused by low resistivity source.For the measuring direction which has abnormal annual variation form,dynamic annual feature is opposite to the one in direction of normal annual variation form.If the dynamic feature is opposite to the annual variation and disturbance amplitude is also greater than annual amplitude,the annual variation will change direction.Disturbance amplitude from metallic conductor is affected by the resistivity and cross-section area,the lower of the resistivity and the larger of the cross-section area,the greater of the disturbance amplitude.     

Investigation of electrical anisotropy in the deep borehole KTB

Information about electrical anisotropy in and around the deep borehole of the German Continental Deep Drilling Program (KTB) has been acquired using three different methods: (1) Under surface temperature and pressure conditions, the azimuthal dependence of resistivity has been measured in horizontal planes perpendicular to the axes of drill cores. (2) Simulated in-situ conditions (high temperature and pressure) have been applied to mini cores drilled at different orientations to structural features such as foliation. The resistivity was then measured in the direction of the mini cores' long axes. (3) Under quasi-in-situ (downhole) conditions, the azimuthal dependence of resistivity has been measured in horizontal planes using the Azimuthal Resistivity Imager (ARI), a down hole logging tool.These investigations have shown that electrical anisotropy is strongly related to structural features. The anisotropy of gneisses is the highest due to the strong foliation (averaged anisotropy factor 2.8). Lamprophyres are most isotropic (factor 1.2). Most rocks from KTB are amphibolites with intermediate anisotropies (factor 1.3). The anisotropies measured under surface temperature and pressure conditions and quasi-in-situ conditions on amphibolites are approximately the same, whereas anisotropies under simulated in-situ conditions are partially higher.Furthermore, the anisotropies measured on the cm to m scale using the various methods outlined above shall be compared with large scale anisotropies of the km scale measured around KTB using Magneto Tellurics.     

磁化率张量的自益再采样法

岩石的磁化率张量的统计处理，传统的处理方法是用费歇尔法，但对于火成岩等磁化率各向异性较弱的岩石，这种方法存在着局限性。本文介绍了“自益再采样法”，并与传统的方法进行了对比，指出，自益再采样法不仅有助于消除测量过程中的误差，而且它是建立在数据的实际分布之上的，能更准确地反映磁化率轴的实际分布方位以及方位的分散程度。；     

Anisotropy of P-wave Velocity in the Upper Crust of the Central External Dinarides

The anisotropy of Pg-wave velocity in the area of the central External Dinarides is measured by using arrival time data as reported by local and regional seismological stations. The observed velocity varies between 5.73 km/s (in the ESE-WNWdirection) and 6.20 km/s (in the SSW-NNE direction), indicating azimuthal anisotropy with symmetry axis azimuths of 23°±n 90°. These closely match the orientation of the principal stress axes in the region, as revealed by analyses of available fault-plane solutions. The observed anisotropy may be modelled by assuming a system of vertical/subvertical cracks in the upper crust, aligned under the influence of the regional tectonic stress field.