A GENERAL TECHNIQUE FOR THE DIRECT INTERPRETATION OF RESISTIVITY DATA OVER TWO-DIMENSIONAL STRUCTURES*

The indirect method of interpreting resistivity data is capable only of limited success because of the difficulty of calculating type curves for complex structures. Consequently a need arises for a direct method of interpretation for complex generalized structures. Such a technique for the direct interpretation of apparent resistivity data, obtained by electrical soundings carried out over two-dimensional structures, has been developed from an examination of the Hankel transform of such data. The method is based on the observation that the asymptotic expansion of the Hankel transform is critically dependent upon the minimum distance between the measuring device and the surface of discontinuity in resistivity. The variation in the across strike direction may be mapped by a sequence of depth soundings made parallel to each other and separated in the across strike direction.     

视电阻率法的地形改正

文本提出用理论计算来消除复杂地形对视电阻率影响的方法。借助数学公式解析单元地形的线源畸变电场和点源畸变电场,计算并编制出量板,再将单元地形的视电阻率畸变曲线组合迭加,就能近似地取得复杂地形的改正曲线。文中举出实例,阐明了这种地形改正方法的应用和效果。     

二维频率测深边界单元法正演计算

在二维层状介质中含有横向非均匀体的地电断面情况下，对电偶源的偶极长度作了延长，并选取了适当的基本解及无穷远边界条件，因而将相应的含源谐变电磁场的边值问题转化为较简单的边界积分方程．用边界单元法求出了相应的电磁场分量，进而计算阻抗视电阻率，并对这一算法作了相应的检验．     

THE DETERMINATION OF FILTER COEFFICIENTS FOR THE COMPUTATION OF STANDARD CURVES FOR DIPOLE RESISTIVITY SOUNDING OVER LAYERED EARTH BY LINEAR DIGITAL FILTERING *

The technique of linear digital filtering developed for the computation of standard curves for conventional resistivity and electromagnetic depth soundings is applied to the determination of filter coefficients for the computation of dipole curves from the resistivity transform function by convolution. In designing the filter function from which the coefficients are derived, a sampling interval shorter than the one used in the earlier work on resistivity sounding is found to be necessary. The performance of the filter sets is tested and found to be highly accurate. The method is also simple and very fast in application.     

电偶源频率电磁测深二维阻抗视电阻率计算的源效应校正法

研讨了频率域电磁法中不同源装置的大地电磁测深、线源频率电磁测深和偶极源频率电磁测深阻抗视电阻率的源效应影响特征。在唯象分析的基础上，提出了几种电磁测深法阻抗视电阻率的相互换算法──源效应校正法（大地电磁测深二维ＴＥ极化视电阻率和其它两种电磁法的赤道装置二维阻抗视电阻车）。模型试验表明，利用这一源效应校正法可以由大地电磁二维视电阻率近似地计算出线源频率电磁二维阻抗视电阻率。这一方法被尝试应用于由线源频率电磁二维阻抗视电阻率估算偶极源频率电磁二维阻抗视电阻率。     

Quantitative interpretation of self potential anomalies of some specific geometric bodies

Summary Some direct and quantitative methods of SP anomalies caused by some specific geometric bodies have been developed in this paper. The models of current sources which have been considered are i) single pole, ii) a doublet, iii) a pair of single point poles separated by a horizontal distance, iv) single finite line pole, v) single infinite line pole and vi) two similar double infinite vertical line poles separated by a horizontal distance.     

Improved imaging of a karst aquifer using focused source electromagnetic and differentially normalized method: a qualitative analysis

Geophysical investigations using conventional techniques applied to groundwater exploration can often present strong limitations involving high financial costs, complex acquisition logistics and high ambiguity in results. Dispersion of the electric current flow, induced polarization) effects, cultural noises and shallow lateral heterogeneities represent the main problems faced by geoelectric methods in these types of surveys. Moreover, elements such as intrusions and mineralization at different depths may be responsible for signal attenuation as well as high resistivity in unsaturated zones and complex three-dimensional formations or clayey zones cause variations in the electric current. The focused source electromagnetic and differentially normalized method approaches can help to solve some these issues. Aiming at a higher signal-to-noise ratio, the focused source electromagnetic method and approaches of the differentially normalized method, first applied to petroleum exploration, are tested on a groundwater target, in a karst environment sectioned by a diabase dyke. We performed the processing and analysis on real IP resistivity profiling data acquired with two-way dipole-dipole array, guided by magnetic data acquired on the same profile, mapping a diabase dyke. The inversion of focused source electromagnetic method/differentially normalized method was not performed, instead that we converted the induced polarization–resistivity data to a differential signal to qualitatively prove the presence of aquifer. Joint interpretation of focused source electromagnetic method curves and inverted two-dimensional induced polarization–resistivity sections allowed for precise delineation of a conductive zone associated with the karst aquifer, le magnetics allowed for the definition of a neighbour dyke. The techniques have great potential in the aid of groundwater exploration, contributing substantially to the reduction of interpretation ambiguity. Focused source electromagnetic method/differentially normalized method/ approaches show that a simple linear combination of the conventional geoelectric data is able to remove the geological noise and provide the vertical focusing of the electric current.     

有限元直接迭代算法及其在线源频率域电磁响应计算中的应用

采用有限元直接迭代算法实现了线源频率域测深电磁响应的二维正演计算. 首先给出了线源正演问题的有限元直接迭代格式,然后由迭代法进行求解. 在处理奇异源问题上,采用向内递推的组合网格技巧,在源点附近可进行局部加密,并实现粗细网格的对接,从而较好地解决了奇异源附近的计算问题. 还提出一种迭代求取全区视电阻率的方法,避免了远近区的划分. 通过对均匀半空间、层状介质和二维模型电磁响应的计算,获得了与大地电磁测深相似的视电阻率曲线,验证了算法的正确性;通过对计算结果的分析,在理论上说明了线源频率域近区测深的可行性.     

Analysis of the earth resistivity response of buried cables

Current flow in homogeneous and layered conducting media is analyzed for the case where an axial conductor or cable is present. The cable is characterized by a specified axial impedance and is assumed to be infinite in length. Various configurations are chosen such as a current point source in an infinite, semi-infinite, and layered region where the cable is taken parallel to the interface(s). The resulting formulas for the potentials reduce to known cases in the absence of the cable. Using these formulations, we present some concrete calculated examples that are relevant to resistivity probing of perturbed homogeneous and layered structures. Only the two-electrode array is treated, but various cable orientations are considered. In general, it is found that a long axial conductor such as a bare cable will distort the potential distribution of the current in a major way. This leads to profound departures from the apparent resistivity curves calculated for idealized homogeneous and layered structures.U.S. Government Work not subject to U.S. copyright.     

GEOELECTRIC AXIAL DIPOLE SOUNDING CURVES FOR A CLASS OF TWO-DIMENSIONAL EARTH STRUCTURES*

With the aim of studying the behaviour of geoelectric axial dipole vertical soundings over complex geology, a systematic theoretical approach is presented for a class of earth structures characterized by horizontal and vertical parallel boundary planes. The two-dimensional cylindrical bodies of infinite length and rectangular cross-section are constrained to have resistivities satisfying Alfano's condition at every intersection line of the graticule, in order to adopt the image-point theory. A detailed analysis is performed for models with any number of horizontal boundaries and two vertical discontinuities. The apparent resistivity formulas are obtained and selected apparent resistivity curves are drawn for different parameter combinations and various directions of the sounding expansion axis. The class under consideration contains as a particular case the HVC model elaborated in Alpin's monograph, where only a small collection of master curves is available for the axial array. The reconstruction of those curves by the present formulation shows the existence of large discrepancies. A test based on the transformation to equivalent half-Schlumberger sounding curves supports the conclusion that an unidentified error must exist in some part of the theoretical approach of the Russian researchers. Finally, some field sounding curves based on geothermal and volcanological surveys are presented and interpreted by complete curve matching, essentially to show the applicability of the theoretical solutions.     

THE QUANTITATIVE INTERPRETATION OF DIPOLE SOUNDINGS BY MEANS OF THE RESISTIVITY TRANSFORM FUNCTION*

In this paper a theorem is demonstrated which allows—after the introduction of a suitable dipole kernel function or dipole resistivity transform function—to write the apparent resistivity function as an Hankel transformable integral expression. As a practical application of the theorem a procedure of quantitative interpretation of dipole soundings is suggested in which the dipole resistivity transform function obtained after inversion of the original dipole apparent resistivity data is used to control the goodness of the set of layering parameters which have been derived with our previous method of transformation of dipole sounding curves into equivalent Schlumberger diagrams.     

THE AUTOMATIC FITTING OF A RESISTIVITY SOUNDING BY A GEOMETRICAL PROGRESSION OF DEPTHS *

The problem of interpretation of geoelectrical resistivity soundings has been studied and a compromise method of interpretation developed. A simple depth sequence of increasing bed thickness is taken and resistivities assigned to the various layers so that a good fit is obtained to the observed apparent resistivity curve. An initial estimate of the resistivities is made by a direct method but the model is then altered to give an improved fit by the indirect method of curve matching. The whole method is computerized and input data to the computer program consists of spacings and apparent resistivities. The method is extremely rapid, about one third of the computation time is devoted to obtaining the resistivity transform by the method of convolution. A normal type of solution, incorporating additional information if necessary, can also be obtained.     

多层地电断面视电阻率的高精度计算

目前地电的方法在我国已应用到地震预报的研究中,为了了解与地下力学场变化相联系的地电变化,需要了解台址下勘探体积内的电阻率分层情况,从而正确地认识各种影响因素对地电阻率的作用。     

THREE-DIMENSIONAL INTERPRETATION OF MULTIPLE-SOURCE BIPOLE-DIPOLE RESISTIVITY DATA USING THE APPARENT RESISTIVITY TENSOR1

The bipole-dipole resistivity technique, which uses a single current source (bipole) to map variations in (apparent) resistivity has been much criticized in the past. A series of 3D models are used to show that the use of two distinct current bipoles in the same location but with different orientations, combined with analysis in the form of a previously defined tensor apparent resistivity, can greatly improve many aspects of bipole-dipole mapping. The model study shows that, for measurement stations more than a few bipole lengths from the current source, the apparent resistivity tensor behaves, to a close approximation, as though the current bipoles are idealized dipoles, and hence is independent of the orientation of the individual current sources used. Any pair of current bipoles (in the same location but with different orientations) can therefore be used to determine the tensor resistivity properties. The invariants of the apparent resistivity tensor have considerable advantages over the normal scalar apparent resistivities. Modelling shows that although the electric field vector corresponding to a single current bipole can be highly perturbed by a local inhomogeneity for some considerable distance beyond the inhomogeneity itself, the tensor invariants are virtually unperturbed beyond the extent of the inhomogeneity. Hence false anomalies, which are a characteristic of apparent resistivity measured using only single current bipole models, are almost completely eliminated by the use of tensor invariants. Of the possible tensor invariants, the invariant given by the square root of the determinant gives the best representation of a buried 3D body. Resistivity anomalies are localized, and occur only over the causative body. Even with complex models involving several buried bodies, the tensor invariants clearly delineate the extent of each body. Outside the bounds of perturbing bodies, the tensor data can be analysed by conventional techniques, for example, to determine layered structure.     

电性源地-井瞬变电磁全域视电阻率定义

地-井瞬变电磁(TEM)方法是在地表发射,探头沿钻孔(井)测量瞬变响应的一种井中物探方法.由于接收探头沿钻孔深人地下,因此能够获得更加可靠的地下目标体信息,尤其当存在低阻覆盖层、浅部矿化等地质干扰,或者勘查深部规模不大之良导矿体时,地-井TEM方法的优势更加明显.相对于磁性源系统,电性源系统发射功率大,辐射面积广,更适合地形复杂地区之深部找矿.本文的目的是研究电性源地-井TEM的近似解释方法,首先给出了层状介质电偶极子在地下的TEM响应,进而通过电偶极子叠加的方式获得了电性源条件下的TEM响应.利用感应磁场与均匀半空间电阻率的单调关系通过反函数定理进行了全域视电阻率定义.理论模型的视电阻率计算结果显示,不同深度测点的视电阻率曲线首支不同,而尾支基本一致.这说明了全域视电阻率是测点周围有限范围内介质的综合反映,所以不同深度测点的视电阻率曲线首支所反映的范围是不同的.绘制了不同时间道视电阻率剖面曲线,以期显示地下电性分布规律.计算结果表明全域视电阻率定义能够基本可靠地反映地下信息,从而为该技术在矿区深部找矿的应用提供有力支持.     

层状介质中膨胀球模拟力源的应力——应变场与视电阻率的关系

本文研究了三层介质中膨胀球模拟震源的应力、应变场以及由这个应力、应变场产生的电阻率变化量和地表视电阻率变化的关系。结果表明,对于一定强度的“震源”,在地电装置探测范围内的应力场随具有不同弹性参数的层位有显著的变化;而应变场却看不出有明显的突变;由地表视电阻率变化量的计算结果揭示了地电异常受台站介质电性结构与力学结构的综合影响。因此作为预报地震的地电阻率法一方面要寻找具有高电阻率——应变灵敏度岩层的台址,另一方面台址的视电阻率变化对各层真电阻率变化要具有好的响应。     

USE OF THE RESISTIVITY METHOD IN GEOLOGICAL MAPPING—CASE HISTORIES FROM RANIGANJ COALFIELD,INDIA*

The effectiveness of the electrical resistivity method has been studied using various configurations with different spacings over the Southern Boundary Fault in the northwestern part, and across a dolerite dyke named “Salma dyke” in the central part, of the Raniganj Coalfield, India. It has been observed that the delineation of the fault and the dyke was made possible under differential tropical weathering conditions. Geoelectric sections across the fault and the dyke have been prepared on the basis of Schlumberger sounding results. In profiling, Wenner, two-electrode, half-Schlumberger and part of Al-Chalabi's configurations were applied over the Southern Boundary Fault. Across the dyke, alpha-, beta-, and gamma-Wenner, Schlumberger, half-Schlumberger and two-electrode configurations were employed. Azi-muthal two-electrode sounding was also studied over the fault. The interpretation of the results of apparent resistivity profiles across the Southern Boundary Fault suggests that the Wenner and two-electrode configurations possess certain diagnostic features which help in mapping a single lithological contact, provided sufficient resistivity contrast exists. Although Schlumberger configuration seems to be quite suitable, other configurations may also be usefully employed over the dyke.     

A NUMERICAL METHOD OF CALCULATING THE KERNEL FUNCTION FROM SCHLUMBERGER APPARENT RESISTIVITY DATA*

A method to calculate the resistivity transform of Schlumberger VES curves has been developed. It consists in approximating the field apparent resistivity data by utilizing a linear combination of simple functions, which must satisfy the following requirements: (i) they must be suitable for fitting the resistivity data; (ii) once the fitting function has been obtained they allow the kernel to be determined in an analytic way. The fitting operation is carried out by the least mean squares method, which also accomplishes a useful smoothing of the field curve (and therefore a partial noise filtering). It gives the possibility of assigning different weights to the apparent resistivity values to be approximated according to their different reliability. For several examples (theoretical resistivity curves in order to estimate the precision of the method and with field data to verify the practicality) yield good results with short execution time independent of shape the apparent resistivity curve.     

RESISTIVITY TYPE CURVES OVER OUTCROPPING VERTICAL DYKE—I *

The theoretical horizontal resistivity profiles over an outcropping vertical dyke for various parameters-electrode spacing, vein-width and resistivity contrast—with inline alpha, beta and gamma-Wenner electrode systems are described. The resistivity profiles present a most bewildering variety of shapes as compared to those in resistivity soundings. The analysis of type curves suggests suitable electrode configuration for detection of wide, moderately wide, and thin veins. The negative apparent resistivity values on the gamma anomalies for resistive vein over certain vein-widths of higher positive values of resistivity reflection factor contradict the concept of apparent resistivity hitherto held.     

FINITE-DIFFERENCE EVALUATION OF APPARENT RESISTIVITY CURVES*

The problem of numerical evaluation of apparent resistivity curves is treated by finite difference modeling. The models proposed are set up in cylindrical coordinates and yield the potential field due to a point source located in a radially symmetric environment. The Schlumberger configuration, widely used for surface measurements, is emphasized. However, the treatment is equally applicable to other similar situations such as the computation of synthetic electric logs when the resistivity of the borehole fluid is different from that of the surrounding uniform or stratified medium. Moreover, the individual layers may not necessarily be isotropic. The medium under investigation is discretized by using a very coarse system of horizontal and vertical grid lines whose distance from the source increases logarithmically; consequently, the physical dimensions of the medium can be made “infinite” without affecting the numerical size of the model. Finer features such as a thin but anomalously resistive or conductive bed which would ordinarily be missed in coarse discretization are accurately taken into account, since the calculations are done in terms of the Dar Zarrouk parameters derived from the exact resistivity distribution of the model. This enables one to compute the potential field by inverting a small sparse matrix. When the medium comprises only a few layers, the efficiency of the finite-difference model is comparable to that of the known analytical methods; for more complicated structures, however, the finite-difference model becomes more efficient. The accuracy of finite-difference results is demonstrated by comparing them with the corresponding analytically obtained data.