平凉台井下地电阻率观测影响系数分析

根据地电阻率影响系数理论,以平凉台4层电性结构为例,分析了井下对称四极地电阻率观测影响系数随深度和极距的变化。结果表明:对于固定的观测极距,影响系数与电极埋深之间的关系复杂;通过计算各层影响系数的大小,认为平凉台井下观测对地表及浅层干扰有较好的抑制作用,其分析结果可为在类似台址电性结构中实施井下地电阻率观测时选择电极埋深和供电极距提供参考。     

A SIMPLE METHOD OF MAGNETOTELLURIC INTERPRETATION *

A simple method of magnetotelluric interpretation is derived using the property that the earth currents flow in horizontal sheets. It is shown that when the depth is taken as two-thirds of the Cagniard's depth of penetration (Cagniard, 1953), the mean resistivity-over the depth is the same as the apparent resistivity of the medium. From a mean resistivity versus depth of penetration plot, resistivity can be easily computed at all depths of the sounding. This method gives satisfactory results in a short time and makes it possible to interpret the soundings over media of several layers. However, the results are inaccurate near maximum and minimum points of the plot as well as in the zone of thin layers.     

THE INFLUENCE OF LATERAL INHOMOGENEITIES ON THE DIPOLE METHODS*

The theory of electrical dipole soundings proved that this method can produce resistivity measurements, which are comparable with those obtained by electrical soundings of the Wenner or Schlumberger type. Their main advantage is the use of short cable lengths, which is important if the depth of penetration should be large. A considerable disadvantage of the dipole method is the great sensitivity to lateral discontinuities. Though these have an influence on the Schlumberger arrangement as well, they can disturb a dipole sounding to such an extent than an interpretation based on a horizontal layer case is no more possible. There are six different dipole arrays, which differ from each other with respect to the angle enclosed by the two dipole orientations-the current dipole AB and the measuring dipole MN. The theoretical comparison of the dipole arrays with the Schlumberger array concerning their sensitivity to lateral discontinuities is a useful basis for the choice of the most suitable configuration. Considering geological subsurface conditions the right choice of a dipole array can give an optimal result, i.e. a dipole sounding for which the sensitivity to lateral discontinuities is as small as possible under the given circumstances.     

DIGITAL LINEAR FILTER FOR COMPUTING TYPE CURVES FOR THE TWO-ELECTRODE SYSTEM OF RESISTIVITY SOUNDING*

In this paper a technique for computing type curves for the two-electrode resistivity soundings is presented. It is shown that the apparent resistivity due to the system can be represented by a convolution integral. Thus, it is possible to apply the principle of digital linear filtering and compute the desired type-curves. The filter function required for the purpose is found to be identical with that used to compute the EM sounding curves for the two coplanar horizontal loop system. It is further shown that from the two-electrode apparent-resistivity expression one can easily derive the apparent resistivities for other configurations. A comparison of depths of investigation for various systems reveals that the two-electrode system has greater depth of investigation than other conventional systems. This is also supported by the field example presented in the end, which illustrates the relative performance of the two-electrode system vis-a-vis the Wenner system.     

Geophysical surveys for identifying saline groundwater in the semi‐arid region of the central Altiplano,Bolivia

In the central part of the Bolivian Altiplano, the shallow groundwater presents electrical conductivities ranging from 0·1 to 20 mS/cm. In order to study the origin of this salinity pattern, a good knowledge is required of the geometry of the aquifer at depth. In this study, geophysics has been used to complement the sparse data available from drill holes. One hundred time‐domain electromagnetic (TDEM) soundings were carried out over an area of 1750 km². About 20 geological logs were available close to some of the TDEM soundings. Three intermediate results were obtained from the combined data: (i) the relationship between the electrical conductivity of the groundwater and the formation resistivity, (ii) geoelectrical cross‐sections and (iii) geoelectrical maps at various depths. The limited data set shows a relationship between resistivity and the nature of the rock. From the cross‐sections, a conductive substratum with a resistivity of less than 1 Ω·m was identified at most of the sites at depths ranging from 50 to 350 m. This substratum could be a clay‐rich formation containing brines. Using derived relationships, maps of the nature of the formation (sandy, intermediate and clayey sediments) were established at depths of 10 and 50 m. Discrimination between sand and clays was impossible where groundwater conductivity is high (>3 mS/cm). In the central part of the area, where the groundwater conductivity is low, sandy sediments are likely to be present from the surface to a depth of more than 200 m. Clayey sediments are more likely to be present in the south‐east and probably constitute a hydraulic barrier to groundwater flow. In conclusion, the study demonstrates the efficiency of the TDEM sounding method to map conductive zones. Copyright © 2001 John Wiley & Sons, Ltd.     

A comparison of peak and plate electrodes in electrical resistivity tomography: application to the chalky groundwater of the Beauvais aquifer (northern part of the Paris basin,France)

Electrical resistivity tomography was used in order to explore an experimental site of the LaSalle Beauvais Polytechnic Institute (France). The test was conducted along a profile line of 315 m length, using 64 electrodes deployed at an inter‐electrode spacing of 5 m, and the data were recorded using gradient, Wenner and pole–dipole arrays. The performance of plate electrodes (non‐conventional flat‐based) is compared with the performance of peak electrodes (conventional spike). The hydrogeophysical investigation of the chalk aquifer system of Beauvais shows that the performance of plate electrodes is satisfactory, leading to inversions of small root‐mean‐square errors. Peak and the plate electrodes were tested before and after injection of a salt tracer in the piezometer of the experimental site. The study demonstrates the usefulness of plate electrodes (efficient, less time consuming) and the possibility of aquifer characterization by a salt tracer. Copyright © 2010 John Wiley & Sons, Ltd.     

DC resistivity methods for determining resistivity in the earth's crust

The use of Schlumberger and dipole arrays for crustal-scale resistivity soundings is considered. Advantages and disadvantages of the two methods are described. The depth to which resistivity may be determined from field measurements is discussed as well as the determination from the sounding curves of various parameters associated with layered structure. The interpretation of experimental data using reference curves as well as two approaches used in computer assisted interpretation are discussed.     

Marine transient electromagnetic sounding of deep buried hydrocarbon reservoirs: principles,methodologies and limitations

The possibility of a time‐domain electromagnetic sounding method using excitation and measurement of vertical electric fields to search for and identify deeply buried reservoirs of hydrocarbons offshore is investigated. The method operates on source–receiver offsets, which are several times less than the depth of the reservoir. Geoelectric information is obtained from the transient responses recorded in the pauses between the pulses of electric current in the absence of the source field. The basics of the method, as well as its sensitivity, resolution, and the highest accessible depth of soundings for various geological conditions in a wide range of sea depths, are analyzed. For the analysis, 1D and 3D geoelectric models of hydrocarbon reservoirs are used. It is shown that under existing technologies of excitation and measurement of vertical electric fields, the highest accessible depth of soundings can be up to 4 km. Technology for the inversion and interpretation of transient responses is demonstrated on experimental data.     

On the Locality of Transient Electromagnetic Soundings with a Single-Loop Configuration

The possibilities of reconstructing two-dimensional (2D) cross sections based on the data of the profile soundings by the transient electromagnetic method (TEM) with a single ungrounded loop are illustrated on three-dimensional (3D) models. The process of reconstruction includes three main steps: transformation of the responses in the depth dependence of resistivity ρ(h) measured along the profile, with their subsequent stitching into the 2D pseudo section; point-by-point one-dimensional (1D) inversion of the responses with the starting model constructed based on the transformations; and correction of the 2D cross section with the use of 2.5-dimensional (2.5D) block inversion. It is shown that single-loop TEM soundings allow studying the geological media within a local domain the lateral dimensions of which are commensurate with the depth of the investigation. The structure of the medium beyond this domain insignificantly affects the sounding results. This locality enables the TEM to reconstruct the geoelectrical structure of the medium from the 2D cross sections with the minimal distortions caused by the lack of information beyond the profile of the transient response measurements.     

深埋电极的地电阻率观测研究

首先研究了四极观测系统装置系数与电极埋深的关系;然后给出了点电流源在3层地壳模型的地表和第二层时,电源所在层的电位的解析表达式;最后将天津宝坻地区的电性结构简化成一个3层模型,计算给出了当地表层和基岩中的电阻率出现变化时,在地表和基岩上层开展四极地电阻率观测结果与供电极距和深度的关系.     

Geoelectromagnetic measurements across the southern Senegal basin (West Africa)

Magnetotelluric and differential geomagnetic sounding surveys, consisting of nine soundings, were performed in 1984 along a 200-km profile across the southern Senegal basin. They were intended to obtain information concerning the resistivity structure of the crust and upper mantle and the distribution of the induced electric currents. Magnetotelluric data indicate that two-dimensional resistivity models are appropriate for the region. The zone above the basement is highly inhomogeneous in geoelectrical structure. Very conductive sediments (0.6-3 ohm m) appear in the Mesozoic-Cenozoic Senegal basin. These sediments lie at depths of up to 4500 m on the west end of the profile. Below this, a modest resistivity material (10–30 ohm m) extends to a maximum depth of about 3000 m. The material at depth on the cast part of the traverse line is thought to be Palaeozoic sediments of the Bove basin. The depth of the magnetotelluric basement lies between about 250 m (in the east) and 4800 m (in the west). The crust is characterised by a drop in electrical resistivity at a depth of 15 km below the east part of the profile. Considering the total section, we observe a general trend towards lower resistivities at depths in excess of 100 km, the transition from 2000 ohm m to about 2 ohm m occurs in the depth range 100 to 175 km. An analysis of the geomagnetic variation field has identified a concentration of telluric current flow beneath the deep basin. It appears that the additional currents flowing in the striking direction of the Senegal basin are largely controlled by sedimentary rocks of high conductivity lying at depths less than 5 km. Model studies show that the local conductivity distribution is able to explain the currents circulating in the thick well-conducting sediments.     

基于电阻率敏感性函数的随钻电磁波测井探测特性研究

为研究地层电性变化时不同区域对电磁波测井响应的贡献分布,从响应信号对地层参数求偏导的角度,给出一种新的电阻率敏感性函数定义,引入模式匹配法对纵向成层、径向非均匀介质敏感性分布进行快速模拟;通过对敏感性函数纵、横向积分给出了单发双收线圈系纵、横向探测范围,研究了井眼、频率、背景地层电阻率等对探测特性的影响.结果表明:敏感性函数能够定量表征响应对地层纵向与径向微观与宏观敏感性,幅度比与相位差敏感性分布形态类似,幅度比较广较深,而相位差分辨率高,敏感范围小;敏感性函数进行径向积分后可表征仪器的探测深度,与伪几何因子对比达到了同样的效果;背景地层电阻率在1~100Ωm变化时,工作频率2 MHz下幅度比探测深度约为0.6~2.3m,相位差为0.3~0.8m,幅度比50%纵向积分敏感性层厚约为0.3~1.6m,相位差约为0.2~0.6m;异常体与背景地层电阻率对比度在1~50变化时,引起的探测深度与敏感性层厚差异约为0.1~0.2m,远小于地层电阻率的影响.     

THE USE OF A SQUARE CONFIGURATION IN RESISTIVITY PROSPECTING *

Consideration is given to the use of a configuration of four electrodes set in a square array for resistivity measurements. It is found that, by passing current successively between different pairs of electrodes, an apparent resistivity can be determined which is both more sensitive to the position of the array centre and less dependent on orientation than the measures usually obtained with colinear arrays of electrodes. At the same time the observations made enable the degree of the departure of local conditions from conditions of lateral homogeneity to be assessed. Theoretical and practical examples of the use of this electrode system are given and the use of the system both as a tool in mapping and in depth investigations is considered. It is shown that provided electrode spacings are suitably arranged the results of a probe carried out using the square array can be interpreted by conventional methods. The system is shown to have particular advantages in the investigation of lateral resistivity variations and the reduced dependence on orientation makes possible the recasting of interpretation data in an orientationally invariant form with a consequent drastic reduction in the number of type curves required for a particular problem.     

A MODIFIED GEOELECTRICAL PROCEDURE USING POLARDIPOLE ARRAYS—AN EXAMPLE OF APPLICATION TO DEEP EXPLORATION *

Dipole soundings are more sensitive to noise caused by lateral and superficial inhomogeneities than Schlumberger soundings. However, the former are preferable for deep explorations in view of the relatively short cables required. The simple solution of carrying out the field work by means of dipole spreads, and to transform the dipole resistivity diagrams into Schlumberger ones by means of proper formulae would be valid only for smooth and regular curves; but often, owing to the presence of lateral noises, the dipole data show a considerable scatter. For such cases a “continuous dipole sounding” method is proposed for which all successive dipoles are contiguous, so that all parts of the profiles are covered and interpolation is not necessary. Obviously the moving dipoles have lengths proportional to their distances, so that they appear equal in the usual bilogarithmic scale. It follows that only polar-dipole arrays may be used. The transition from a dipole to the corresponding Schlumberger apparent resistivity diagram requires an integration constant which is not unequivocally determined. Therefore, the solution is not unique, but all possible derived Schlumberger diagrams have a common part. Similarly, they have some common interpretative results, which may be referred to the original dipole diagram obtained in the field. A special measurement technique is required since the dipole-dipole voltages to be determined are noticeably smaller than the Schlumberger ones. This is true also because dipole soundings are used for great depths and for long distances between the two dipoles.     

THEORETICAL AND EXPERIMENTAL INVESTIGATIONS FOR CAVITY RESEARCH WITH GEOELECTRICAL RESISTIVITY METHODS*

The exact localization of subterranean cavities and the determination of their dimensions is very important for the planning of geotechnical and mining activities. It is a complicated geophysical task often at the limit of detection. Nevertheless geophysical investigation is the only alternative to a dense and expensive grid of boreholes. This report tests the usefulness of geoelectrical resistivity methods for cavity detection under some new aspects. The basis for evaluation was a theoretical analysis of different conventional and focussing measuring arrays and of special arrays for a geoelectrical research between two boreholes. The limit of detectability of a cylindrical cavity of defined cross-section and depth was calculated for the different measuring arrays on the basis of computation of the apparent resistivity ρa. Furthermore, the influence of possible errors (current supply of the electrodes and the distance between the electrodes) is discussed for focussed systems. The second part of the article is directed at the behaviour of the apparent resistivity ρa, the disturbing potential δV_d caused by the cavity and the normal potential δV₀ of the measuring array all in relation to a homogeneous earth. Some new results are presented. In the last part of the article theoretical results are compared with some field measurements.     

井下地电阻率观测中地表电流干扰影响计算

本文把地电阻率台站地下介质简化为水平三层均匀介质模型,以点电流模拟地表干扰电流,针对对称四极观测装置,计算了在不同电性结构中的不同深度观测时,地表电流干扰源对对称四极装置地电阻率观测的影响,计算得到:地表干扰源对电阻率观测的影响取决于电性结构的类型和层参数、供电电极和测量电极的埋深以及避开干扰源的距离.本文研究结果对实施井下电阻率观测中台址电性结构的选择、电极埋深、干扰源避让距离等有参考意义.     

A field test of imaging properties of rotational invariants of the magnetotelluric impedance tensor

A part of the Békés Basin (an extensional sub‐basin of the Pannonian Basin, where the basement under thick Pannonian sediments is well known from deep boreholes and from seismic measurements, and where many magnetotelluric (MT) soundings have been carried out for frequencies ranging from 1 to 10^?3 Hz) was selected as a test area to assess the imaging performances of various apparent‐resistivity definitions computed with rotational invariants of either the real part of the complex impedance tensor, or its imaginary part, or both. A comparison (based on earlier 3D numerical studies) has been made between the magnetotelluric images obtained in this way and the depths to the high‐resistivity basement, as known from boreholes and seismic investigations. The correlation coefficient between the series of basement depth values at 39 MT sites and the apparent‐resistivity values was found to be stronger and high correlation appeared at a shorter period when it was computed with apparent resistivities based on the real tensor rather than with apparent resistivities based on the imaginary tensor. In the light of our studies, ρ_{Re Z} and the impedance phase seem to be more informative than any other combination of magnetotelluric interpretation parameters.     

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Réunion)

 A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Réunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60–100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos, the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise. Received: 3 November 1999 / Accepted: 15 September 1999     

A Low‐Cost,In Situ Resistivity and Temperature Monitoring System

We present a low‐cost, reliable method for long‐term in situ autonomous monitoring of subsurface resistivity and temperature in a shallow, moderately heterogeneous subsurface. Probes, to be left in situ, were constructed at relatively low cost with an electrode spacing of 5 cm. Once installed, these were wired to the CR‐1000 Campbell Scientific Inc. datalogger at the surface to electrically image infiltration fronts in the shallow subsurface. This system was constructed and installed in June 2005 to collect apparent resistivity and temperature data from 96 subsurface electrodes set to a pole‐pole resistivity array pattern and 14 thermistors at regular intervals of 30 cm through May of 2008. From these data, a temperature and resistivity relationship was determined within the vadose zone (to a depth of ~1 m) and within the saturated zone (at depths between 1 and 2 m). The high vertical resolution of the data with resistivity measurements on a scale of 5‐cm spacing coupled with surface precipitation measurements taken at 3‐min intervals for a period of roughly 3 years allowed unique observations of infiltration related to seasonal changes. Both the vertical resistivity instrument probes and the data logger system functioned well for the duration of the test period and demonstrated the capability of this low‐cost monitoring system.     

LES COURBES DAR ZARROUK ET LES SONDAGES DIFFÉRENTIELS DE RÉSISTIVÉ*

Differential and difference soundings of resistivity for real, imaginary and complex values of the newly introduced parameters ν and μ are considered. The transverse and longitudinal differential soundings, known up to the present, correspond with the values ν= 1 and ν=— 1. The point of departure is the generalization of the concept of Dar Zarrouk parameters and the dimensional analogy between D.Z. resistivity and transverse and longitudinal differential resistivities. Properties of generalized D.Z. curves as well as the method of their construction are given. The examples of differential and difference curves are presented and the properties of these curves are discussed; particulary the effective spacings are determined. The problems of the realization of differential and difference soundings are considered; the bisymmetrical and asymmetrical arrays are proposed as a practical ones.