钻井-地表电极联合电阻率观测装置的异常特征研究

本文利用有限单元法对多种钻井-地表电极联合电阻率观测装置在球体、直立板和水平板体等典型地质体模型上的异常响应进行正演计算,并分析其异常特征和分布规律.计算结果表明利用井中和地表电极的联合观测方式进行地下介质探测,可在兼顾地质体的水平分辨能力的同时,提高电阻率测深在垂向上分辨率.就水平分辨率而言,四极（偶极）装置对地质体的水平定位能力最强,井-地三极装置次之,地-井三极的分辨能力最差.就垂向分辨能力来说,各种观测装置的分辨率相差不大.     

Field test to compare 3D imaging capabilities of three arrays in a site with high resistivity contrast regions

The applicability of three kinds of electrode configurations used to delineate a buried horizontal pipe was studied. A 3D resistivity imaging survey was carried out along eight parallel lines using pole-pole, pole-dipole, and dipole-dipole arrays with 1m minimum electrode spacings. Roll-along measurements were carried out to cover a rectangular grid. The 2D and 3D least squares algorithms based on the robust inversion method were used in the inversion of the apparent resistivity data sets. The 2D inversion of data sets could not delineate the orientation and dimension of the subsurface anomalies clearly. To obtain more accurate results, a 3D joint inversion of the pole-pole and pole-dipole data sets was performed, as well as of pole-pole and dipole-dipole data sets. In this case, both horizontal and vertical dimensions of subsurface structures were resolved. The resulting model obtained from each array was compared to those of joint inversion method. The result showed that the horizontal resolution does not improve so much as that in the vertical direction when joint inversion is applied.     

Method to obtain standard pseudosections from pseudo pole-dipole arrays

This study deals with electrical resistivity tomography data and it is addressed at obtaining, from linear combinations of data, other datasets not directly measured. The method presented here allows performing the fieldwork systematically, without deciding the type of the more suitable common electrode array (e.g. Wenner, Wenner-Schlumberger, dipole-dipole or multiple-gradient) until or even after the interpretation time. The electrode configuration used by this method is denoted as pseudo pole-dipole array, because it is based on the standard pole-dipole one, but avoiding the disadvantage of locating the remote electrode far away from the profile (to “infinity”). The pseudo pole-dipole datasets can be acquired with common equipment using standard pole-dipole recording sequences. Once the desired datasets have been calculated, they can be inverted using standard interpretation software. The procedure used allows a data quality control to be introduced that is similar to the one based on normal and reciprocal measurements. To assess the method we considered noise-contaminated model responses as well as field data. It has been applied to build dipole-dipole and Wenner-Schlumberger datasets. Results show the suitability of both, the proposed method and the quality control.     

DATA QUALITY ASSESSMENT OF TIME-DOMAIN IP DECAY CURVES*

Data quality assessments have been made on time-domain induced polarization decay curves taken over sulphide ore bodies of Cyprus. The data were taken with the dipole-dipole, pole-dipole and symmetrical quadripole arrays. The results indicate that at larger array separations decay curves taken with the dipole-dipole and pole-dipole arrays are often seriously affected by noise, but that curves taken with the symmetrical quadripole array are not so affected. Optimum pseudo-sections designed to minimize the noise show that the noise results in spurious highs and lows. This work was carried out as part of a NERC research studentship at the University of Leicester.     

Use of four-electrode arrays in three-dimensional electrical resistivity imaging survey

The objective of this paper is to investigate the applicability of four-electrode arrays in 3D electrical resistivity imaging survey. A 3D resistivity imaging survey was carried out along fourteen parallel lines using dipole-dipole, Wenner-Schlumberger, and Wenner arrays with 2 m minimum electrode spacings. Roll-along measurements using a line spacing of 1 m were carried out covering a grid of 20 × 14 electrodes. The 3D least squares algorithm, based on the robust inversion method, was used in the inversion of the 3D apparent resistivity data sets. The results show that the 3D electrical resistivity imaging survey using the Wenner-Schlumberger and the dipole-dipole arrays, or the Wenner and the dipole-dipole arrays, in combination with an appropriate 3D inversion method, can be highly useful when the site conditions do not allow using the pole-pole or pole-dipole arrays.     

COMPARATIVE FIELD PERFORMANCE OF ELECTRODE ARRAYS IN TIME-DOMAIN INDUCED POLARIZATION PROFILING *

Five test profiles in time-domain induced polarization measurements over conducting sulphide and graphite mineralisations were run with the two-electrode, three-electrode, modified unipole, pole-dipole, dipole-dipole, Schlumberger, and Wenner α and β configurations. The results show that, compared to the other electrode systems, the simplest two-electrode array produces the largest anomalies with the smallest of spacings.     

不同电极阵列联合反演在古墓探测中的应用

电阻率层析成像是一种广泛应用在水文、考古和地质等浅地表勘探领域的地球物理方法。为了增强电阻率层析成像的分辨率、应对复杂的地质问题,本文提出基于雅可比矩阵的不同电极阵列直流电阻率数据的加权联合反演算法,并以温纳和偶极-偶极电极阵列数据为例,在理论模型和古墓探测的野外实例中测试该算法的有效性。结果表明,加权联合反演结果的横向和纵向分辨率都优于单一电极阵列的反演结果,并在实例中缓解“U形”电极阵列的固有缺陷、减少反演模糊性、更好地约束墓室宽度的反演结果。      

2D电阻率正反演成像在水平和垂直模型上的异常响应研究

2D电阻率成像技术近年来被广泛应用于工程、水文、环境和矿产等领域,在实际应用中它具有多种的装置类型,不同的装置类型对特定的地质情况有不同的应用效果.本文通过采用Wenner、Wenner-Schlumberger和dipole-dipole三种排列装置在一个水平和一个垂直模型上正演模拟和块反演,揭示了不同的排列装置在水平结构和垂直结构的异常响应,并对高阻体和低阻体进行了模拟.结果显示当采用Wenner,Wenner-Schlumberger和dipole-dipole数据采集技术时,不同的模型结构异常响应有明显的不同,三种排列类型对地下低阻体的分辨能力均高于高阻体,Wenner装置和dipole-dipole装置分别对水平层状结构和垂直结构有较好的分辨力,wenner-Schlumberger装置对水平层状结构有好的分辨力,对垂直结构有中等的分辨能力.     

环工作面三维直流电阻率法研究及应用

本文在矿井直流电法理论基础上,推导环工作面三维直流电法测量深度公式,得到环工作面三维直流电法记录点坐标位置.基于以上理论,建立了环工作面三维直流电法异常体模型,分析了环工作面三维直流电法二极装置、三极装置、四极装置异常体响应特征,得出环工作面三维直流电法各装置敏感度低的主要原因是缺少必要的数据采集点.为了提高环工作面三维直流电法分辨率,最大限度增加工作面内侧的采集数据点.基于以上原因本文提出了对角偶极装置采集模式及多装置的数据叠加反演处理方法,以增加工作面内采集数据点,提高三维勘探的精度.本文通过某矿8802工作面底板富水性探查实例,说明该方法的有效性.     

Pilot‐scale field validation of the long electrode electrical resistivity tomography method

A validation experiment, carried out in a scaled field setting, was attempted for the long electrode electrical resistivity tomography method in order to demonstrate the performance of the technique in imaging a simple buried target. The experiment was an approximately 1/17 scale mock‐up of a region encompassing a buried nuclear waste tank on the Hanford site. The target of focus was constructed by manually forming a simulated plume within the vadose zone using a tank waste simulant. The long electrode results were compared to results from conventional point electrodes on the surface and buried within the survey domain. Using a pole‐pole array, both point and long electrode imaging techniques identified the lateral extents of the pre‐formed plume with reasonable fidelity but the long electrode method was handicapped in reconstructing vertical boundaries. The pole‐dipole and dipole‐dipole arrays were also tested with the long electrode method and were shown to have the least favourable target properties, including the position of the reconstructed plume relative to the known plume and the intensity of false positive targets. The poor performance of the pole‐dipole and dipole‐dipole arrays was attributed to an inexhaustive and non‐optimal coverage of data at key electrodes, as well as an increased noise for electrode combinations with high geometric factors. However, when comparing the model resolution matrix among the different acquisition strategies, the pole‐dipole and dipole‐dipole arrays using long electrodes were shown to have significantly higher average and maximum values within the matrix than any pole‐pole array. The model resolution describes how well the inversion model resolves the subsurface. Given the model resolution performance of the pole‐dipole and dipole‐dipole arrays, it may be worth investing in tools to understand the optimum subset of randomly distributed electrode pairs to produce maximum performance from the inversion model.     

A new method for geoelectrical investigations underwater1

A new electrical method is proposed for determining the apparent resistivity of multi-earth layers located underwater. The method is based on direct current geoelectric sounding principles. A layered earth model is used to simulate the stratigraphic target. The measurement array is of pole-pole type; it is located underwater and is orientated vertically. This particular electrode configuration is very useful when conventional electrical methods cannot be used, especially if the water depth becomes very important. The calculated apparent resistivity shows a substantial quality increase in the measured signal caused by the underwater targets, from which little or no response is measured using conventional surface electrode methods. In practice, however, different factors such as water stratification, underwater streams or meteorological conditions complicate the interpretation of the field results. A case study is presented, where field surveys carried out on Lake Geneva were interpreted using the calculated apparent resistivity master-curves.     

Pole-dipole resistivity sounding technique for shallow investigations in hard rock areas

Theory of the pole-dipole resistivity sounding technique and its application in the hard rock areas for shallow groundwater exploration is presented in this paper. The different components of electric field produced by the point source of current, situated over the ground surface, are measured by a dipole placed at a large distance from the source. The theory of the method is rather simple, suggesting two configurations, namely radial pole-dipole and axial pole-dipole. Theoretical expressions derived for the apparent resistivity over layered Earth are directly related to the Schlumberger apparent resistivity, whereas expressions for geometrical factor for pole-dipole and Schlumberger configuration are different. The proposed technique has been tested in actual field conditions having different rock types. A few examples are presented along with Schlumberger sounding curves which confirm the applicability of the proposed sounding technique.     

A numerical comparison of 2D resistivity imaging with 10 electrode arrays

Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L₁‐norm) inversion and smoothness‐constrained least‐squares (L₂‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L₂‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.     

INDUCED POLARIZATION RESPONSE OF A HORIZONTALLY MULTILAYERED EARTH WITH NO RESISTIVITY CONTRAST*

The induced polarization response of a horizontally multilayered earth with no resistivity contrast can rapidly be calculated on a desk calculator or minicomputer for any electrode array. The formulation is a simple series summation of the products of weighting coefficients and the true induced polarization responses for each of the layers. The coefficients are directly derivable from the corresponding resistivity model. This series approach to IP formulation was originally described by Seigel but has not been treated extensively in the present-day geophysical literature. This method can be applied to either time or frequency domain induced polarization measurements. Once the coefficients are known, apparent induced polarization response can readily be obtained by judicious substitution of known, suspected, or assumed values of the true induced polarization of each layer. Basic formulation is presented for the IP potential coefficients (pole-pole or two array) with no resistivity contrast between the layers. From these coefficients, response of any number of layers for any electrode array can be obtained by suitable differentiation. Some examples of Wenner array for a three-layered earth and dipole-dipole array for a four-layered earth are used to illustrate the application. The results of this technique are valid for many natural situations of modest resistivity contrast. However, they definitely cannot be used if there are highly contrasting resistivity layers present. Such an approach is conceptually simple and is useful for survey planning, checking or setting the “depth-of-penetration”of a given array. For field induced polarization data that fits reasonably well to the no-resistivity-contrast model, this simple approach facilitates quantitative interpretation.     

超高密度电阻率法在土洞、煤窑采空区和岩溶勘探中应用研究

超高密度电阻率法是Australia ZZ Resistivity Imaging研发中心积近十年的研究成果提出的全新的电阻率法勘探理念,是电法勘探从手工时代向计算机时代的一次跨越式发展,采用一次布极,采集所有可能组合电极间的电位信息,信息量较常规电法在信息量上翻了40倍以上,并且避开视电阻率概念,利用处理软件直接反演获取断面真电阻率值,针对目前还没有行之有效的地面物探方法可以勘探清楚的土洞、煤窑采空区和岩溶等复杂地质问题,我们采用该方法分别对机场岩溶和土洞、煤矿运煤巷道、隐伏小煤矿采空区和复杂地质病害隧道进行了试验勘探,经钻探验证,准确率在60%以上,达到国内领先水平,取得了较好效果,对复杂地质的工程勘探具有一定参考价值.     

三维激发极化法反演研究及其在金属矿勘探中的应用

激发极化法在金属矿、硫化矿等资源勘探方面应用较广.随着勘探设备与计算机硬件的发展,所采集到的数据具有观测数据量大、布极方式多样等特点.针对实测数据的特点,我们研究完成了基于并行技术的激发极化法对数反演算法,该算法具有如下特点：（1）通过压缩存储技术和并行技术的集成实现了可以处理大数据特征的反演算法;（2）利用对数反演来约束每次计算得到新模型的电导率恒为正,充电率在0到1之间变化,从而保证反演的稳定性和可靠性.本文首先设计了两组对比模型进行试验,通过对不同区块数据采用不同加权的方式来减弱噪声对反演结果影响的效果;其次,采用并行技术提高了反演的计算速度,并利用理论模型分析了不同电极装置对反演分辨率的影响.最后,在甘肃某金属矿开发前景区利用激发极化法开展了中梯装置的采集工作,利用加权后的观测数据反演推断出了测区金属矿开发靶区的大致位置及分布特征.

     

The effect of subsurface pipes on apparent‐resistivity measurements

Subsurface conducting pipes can be either a target or a noise source in geophysical surveying. Their effect as a noise source in resistivity imaging can be so severe as to render the geophysical data uninterpretable. A method is developed here for identifying, locating and removing the effects of subsurface conducting pipes from image data, thus revealing the background resistivity structure.
A previously known analytic solution for the potential distribution produced by current injection in a uniform half‐space containing an infinitely long conducting cylinder is used to calculate apparent resistivities corresponding to electrode arrays on the surface of the half‐space. Most results concern the Wenner array and an examination is made of the effects produced by varying the electrode spacing and the depth, size and orientation of the pipe with respect to the array. A method is developed for locating pipes in resistivity image data by cross‐correlation of the analytic solution with the measured field data. Pipe effects are then removed by multiplying each datum point in the measurements by the reciprocal of the corresponding value in the analytic solution. The success of the method is demonstrated by applications to synthetic data sets involving one or two pipes embedded in non‐uniform half‐spaces.
In further examples, the method is applied to some measured resistivity images from an ex‐industrial site (a former oil distribution terminal), where an electromagnetic survey had previously revealed a labyrinth of underground pipes. The method is shown to be successful in removing the effects of the pipes to reveal the underlying geology.     

THE OFFSET SYSTEM OF ELECTRICAL RESISTIVITY SOUNDING AND ITS USE WITH A MULTICORE CABLE*†

One of the most troublesome problems in resistivity sounding is caused by the spurious effects of near-surface lateral resistivity variations. It has been found that the effects can be strongly reduced by measuring earth resistances at two electrode array positions such that in one position the lateral effect is a positive contribution to the total measured signal and in the other the contribution is negative. The subsequent combination of the two measurements virtually eliminates the lateral effect. The same technique can also result in a significant reduction in the total number of necessary electrode positions. Consequently the method can be used with a multicore cable system. In this way, in addition to improving the accuracy of the measured results, time and manpower are considerably reduced. Details of the method, examples and comparisons with traditional sounding techniques are given.     

Two‐dimensional joint inversions of cross‐hole resistivity data and resolution analysis of combined arrays

In this study, a new two‐dimensional inversion algorithm was developed for the inversion of cross‐hole direct current resistivity measurements. In the last decades, various array optimisation methods were suggested for resistivity tomography. However, researchers have still collected data by using classical electrode arrays in most cross‐hole applications. Therefore, we investigated the accuracy of both the individual and the joint inversion of the classical cross‐hole arrays by using both synthetic and field data with the developed algorithm. We showed that the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole electrode arrays gives inverse solutions that are closer to the real model than the individual inversions of the electrode array datasets for the synthetic data inversion. The model resolution matrix of the suggested arrays was used to analyse the inversion results. This model resolution analysis also showed the advantage of the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole arrays. We also used sensitivity sections from each of the arrays and their superpositions to explain why joint inversion gives better resolution than the any individual inversion result.     

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免, 其对直流电阻率法勘探影响巨大.近年来, 电阻率三维正演取得诸多进展, 特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟, 但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格, 并考虑到应用GPS/GNSS时, 区域地球物理调查中可非规则布设测网的实际特点, 实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性, 可应用于复杂野外环境下的三维电法勘探.