RESISTIVITY SOUNDING ON A MULTI-LAYERED EARTH WITH TRANSITIONAL LAYERS. PART I: THEORY*

The electrical potential generated by a point source of current on the ground surface is studied for a multi-layered earth formed by layers alternatively characterized by a constant conductivity value and by conductivity varying linearly with depth. The problem is accounted for by solving a Laplace's differential equation for the uniform layers and a Poisson's differential equation for the transitional layers. Then, by a simple algorithm and by the introduction of a suitable kernel function, the general expression of the apparent resistivity for a Schlumberger array placed on the surface is obtained. Moreover some details are given for the solution of particular cases as 1) the presence of a infinitely resistive basement, 2) the absence of any one or more uniform layers, and 3) the absence of any one or more transitional layers. The new theory proves to be rather general, as it includes that for uniform layers with sharp boundaries as a particular case. Some mathematical properties of the kernel function are studied in view of the application of a direct system of quantitative interpretation. Two steps are considered for the solution of the direct problem: (i) The determination of the kernel function from the field measurements of the apparent resistivity. Owing to the identical mathematical formalism of the old with this new resistivity theory, the procedures there developed for the execution of the first step are here as well applicable without any change. Thus, some graphical and numerical procedures, already published, are recalled. (ii) The determination of the layer distribution from the kernel function. A recurrent procedure is proposed and studied in detail. This recurrent procedure follows the principle of the reduction to a lower boundary plane, as originally suggested by Koefoed for the old geoelectrical theory. Here the method differs mainly for the presence of reduction coefficients, which must be calculated each time when passing to a reduced earth section.     

视电阻率法的地形改正

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

CSAMT全区电阻率法数值模拟及应用探讨

常规的可控源电磁法理论在计算视电阻率公式上,多半采用其电磁场的渐近特征,难以直接反映全区视电阻率的值,及直观地显现地下介质的地质构造.文中采用水平偶极子激发的电磁场,提出了电场的全区精确表达式,直接计算出大地电阻率.利用汉克尔数值滤波算法和逆样条插值算法对水平层状电磁场进行正演计算,并与计算的卡尼亚视电阻率的对比和野外试验结果表明：该方法的结果在远区等价卡尼亚电阻率,在近区和过渡带则明显地改善了卡尼亚电阻率的非波场区场畸变,从而能更好地接近基底的真电阻率,更形象地反映了地下介质的垂向电性变化.     

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.     

井间电阻率层析成像研究新进展

由于井间电阻率层析成像方法具有施工方便、成本低、对钻孔无破坏作用且穿透深度也较大等优点，因而得到越来越广泛的应用，在过去几年内对该方法研究的重点在于提高计算效率和改善成像质量，本文从井间电阻率层析成像的基本原理出发，对该方法近年的新进展进行了比较和综述，还介绍了概率层析成像算法的原理，认为今后几年内在此领域内的理论研究重点将是探索有效的反演方法和提高成像质量。     

FAST AUTOMATIC PROCESSING OF RESISTIVITY SOUNDINGS*

The difficulty to use master curves as well as classical techniques for the determination of layer distribution (e_i, ρ_i) from a resistivity sounding arises when the presumed number of layers exceeds five or six. The principle of the method proposed here is based on the identification of the resistivity transform. This principle was recently underlined by many authors. The resistivity transform can be easily derived from the experimental data by the application of Ghosh's linear filter, and another method for deriving the filter coefficientes is suggested. For a given theoretical resistivity transform corresponding to a given distribution of layers (thicknesses and resistivities) various criteria that measure the difference between this theoretical resistivity transform and an experimental one derived by the application of Ghosh's filter are given. A discussion of these criteria from a physical as well as a mathematical point of view follows. The proposed method is then exposed; it is based on a gradient method. The type of gradient method used is defined and justified physically as well as with numerical examples of identified master curves. The practical use for the method and experimental confrontation of identified field curves with drill holes are given. The cost as well as memory occupation and time of execution of the program on CDC 7600 computer is estimated.     

New advances in regularized inversion of gravity and electromagnetic data

The interpretation of potential and electromagnetic fields observed over 3D geological structures remains one of the most challenging problems of exploration geophysics. In this paper I present an overview of novel methods of inversion and imaging of gravity and electromagnetic data, which are based on new advances in the regularization theory related to the application of special stabilizing functionals, which allow the reconstruction of both smooth images of the underground geological structures and models with sharp geological boundaries. I demonstrate that sharp-boundary geophysical inversion can improve the efficiency and resolution of the inverse problem solution. The methods are illustrated with synthetic and practical examples of the 3D inversion of potential and electromagnetic field data.     

Higher-order statistics correlation stacking for DC electrical data in the wavelet domain

DC (direct current) electrical and shallow seismic methods are indispensable to the near surface geophysical exploration, but the near surface areas are very difficult environments for any geophysical exploration due to the random noise caused by near surface inhomogeneities. As a new algorithm based on higher-order statistics theory, the higher-order correlation stacking algorithm for seismic data smoothing in the wavelet domain has been developed and applied efficiently to filter some correlation noise that the conventional second-order correlation stacking could not inhibit. In this paper, this higher-order statistics correlation stacking technology is presented for DC electrical data in wavelet domain. Taking into account the single section and multiple section data, we present two new formulations of correlation stacking for DC electrical data. Synthetic examples with Gaussian noise are designed to analyze the overall efficiency of the new algorithm and to determine its efficacy. Meanwhile, comparison with the traditional least-squares optimization inversion method for field examples from electrical imaging surveys and time-domain IP measurement in China shows its significant advantages. The quality of the new algorithm also has been assessed by physical simulation experiments. This new technology in DC electrical exploration measurements provides a new application in engineering and mining investigation.     

The Application of Geotomography in Engineering

It seems that electromagnetic tomography is quite suitable to be used in metal and groundwaterexploration. especial for limestone areas where the absorption coefficient of electromagnetic wave is smalland geological body has a large property contrast with country rock it should be pointed out that there exitsome uncertainties in geological interpretation of EM tomogram, corresponding laboratory tests for rockphysical property are always necessary.Resistivity tomography provides us a new useful…     

各向异性层状介质中视电阻率与磁场响应研究

针对任意各向异性地层,利用极向型和环向型标量位函数,导出相应的直流电视电阻率和磁电阻率的磁场响应关系.计算了各向异性地层的直流电视电阻率和磁电阻率响应特征,重点分析了电阻率测深方法与磁电阻率探测方法对地下各向异性介质的探测能力.文中采用状态矩阵的分析方法,首先采用极向型和环向型标量位构造了各向异性层状介质电场与磁场的通解,利用各层界面电场、磁场的连续性及地面激励源的耦合条件,推导了不同层之间电磁场的状态矩阵,建立了空间电场与磁场的递归计算关系.其次,针对递归计算中指数项数值计算的不稳定性,借用状态矩阵的性质,导出了将不稳定指数计算项转化为稳定的指数项的转换关系.针对横向各向同性(TI)介质中极向位与环向位解耦的特点,导出了电磁场的直接积分解.最后,采用解析解验证了算法的正确性,给出了多层各向异性地层模型的视电阻率和磁场响应曲线,分析了直流电法探测裂缝性地层、估计裂缝分布性状的可能性.     

Resistivity exploration for altered zone at otake geothermal area — Japan

Geothermal electric power plant with a capacity of 10 mw has been now constructed by the Kyushu Electric Power Co. Ltd. at Otake, northern part of the Kyushu Island. Resistivity exploration, one of geophysical prospecting applied to this field in order to investigate the subsurface structure accompanied with natural steams, indicated low resistivity layer corresponding to altered zone of andesite due to hydrothermal agents. For the rapid interpretation of resistivity sounding curves, a method of calculating the standard curve for the multiple layered earth is developed in the utilization of an electronic digital computer. Some examples to show the application of this method are given, and tinally the result is described of resistivity measurement at Otake geothermal area.     

三维可视化在CSAMT勘探中的应用

CSAMT法是一种强有力的地球物理勘探手段,具有测量迅速、探测深度大、信噪比高等优点,已被广泛应用于金属矿、油气田、地下水、以及地热田的勘探中.在CSAMT勘探资料解释方面,传统的二维电阻率断面很难清晰、直观地揭示3D(三维)地电结构的分布特征,为此应加强3D可视化应用技术的研究.本文详细的介绍了3D可视化的建模过程,将3D构造的建模分为空间曲面剖分、空间曲面拟合插值、空间曲面拓扑关系生成、3D结构网格剖分、属性建模等过程.并以一个CSAMT勘探实例,展现出3D可视化方法技术的突出优势.三维可视化可以多角度、多层面展现地下3D电性结构的分布特征,使我们更加清晰、直观地了解3D地质体的空间范围与几何形态.     

AIRBORNE RESISTIVITY MAPPING USING A MULTIFREQUENCY ELECTROMAGNETIC SYSTEM*

Airborne electromagnetic methods are most commonly used in mineral exploration. However, new developments, such as multifrequency capability and digital on-board field recording, as well as improvements in instrumentation resulting in high signal-to-noise ratios in recorded data, have made their application in geological mapping possible. A three-frequency airborne EM survey carried out over an area northwest of Timmins, Ontario, was interpreted in terms of thickness and resistivity of the layers of a two-layer earth section. Since both in-phase and quadrature components are measured, this provides six independent parameters at each point in space. Based on prior geological information and a preliminary interpretation of the field records, two two-layer models of the subsurface seemed to be appropriate for most of the survey area. An automatic computerized interpretation procedure was devised to interpret the field data at each point in terms of thickness and resistivity parameters of those two models. When the geology is more complex, the data do not fit the models and no interpretations are made. Two maps illustrating the variation of the resistivity and the thicknesses of the layers were constructed from the interpreted data. These maps agree with the known geological information about the distribution of glacial clay in the area. Areas where the layered models do not fit are known to be areas where the geology is complex with a large number of dykes and other lateral inhomogeneities. The study shows that multifrequency airborne EM surveys can be very useful in geological mapping over inaccessible terrain and can significantly help the mapping geologist where outcrops are scarce.     

电容耦合式电阻率法在城市遗址考古勘探中的应用

城市密集柏油公路网和坚硬水泥地面给地下遗址考古勘探造成了很大困难,传统的高密度电阻率法需要插置电极,无法在城市环境中使用。电容耦合式电阻率法采用非侵入式的方式,通过特制电极与大地形成电容器产生电流进行电法勘探,特别适合城市硬地表介质的条件。为了检验电容耦合式电阻率法在城市遗址考古勘探中的应用效果,在地处杭州市主城区的吴越捍海塘遗址进行了电容耦合式电阻率法探测试验,采集多个收发距的电法数据,对数据进行反演处理,将反演结果与考古发掘结果进行对比分析。结果表明,反演的电阻率图像推断结果与考古发掘结果基本吻合,电容耦合式电阻率法是城市遗址考古勘探的有效方法,具有广泛的应用前景。      

发震断层面上潮汐有效剪切应力增量的计算与加卸载的确定

发震断层面上潮汐有效剪切应力增量的数值计算是应用响应比理论的一个关键性基础问题。本文根据潮汐形变理论及弹性力学理论，导出了直接、严密和实用的潮汐有效剪切应力增量的数值计算公式，给出了判断加载、卸载区间的直接公式——潮汐有效剪切应力增量对时间-阶微分的解析表达式，最后还给出了一些应用实例。     

基于NMM随钻电磁波传输信道特性分析

基于数值模式匹配法(NMM)与源等效原理,建立了随钻电磁波传输理论模型,可以同时考虑介质的径向和轴向非均质性,便于分析套管、泥浆、高导层、高阻层对信号传输的影响.通过理论计算结果与试验结果对比,验证了理论模型的正确性.在此基础上对地层电阻率、工作频率、套管、泥浆、钻柱导电性、高导层和高阻层对信号传输的影响进行了计算分析.     

MT方法中Basokur视电阻率定义的研究

电法勘探中不同的视电阻率定义可以带来不同的反演效果,好的视电阻率定义可以提高反演分辨率。传统的Cagniar视电阻率(1953)只是利用了波阻抗的模值,而Basokur 1994年提出的视电阻率利用波阻抗的实部和虚部。本文做了4个不同模型的Basokur视电阻率响应及反演,结果显示新的视电阻率对薄层有很好的反应,对常规模型也有不错的表现。     

高密度电阻率法在西藏日喀则地区隐伏断裂探测中的应用

以跨谢通门—青都断裂的两条高密度电阻率法探测资料为基础, 对高密度电阻率法在青藏高原日喀则地区隐伏断裂探测中的首次应用进行了详细介绍． 所获取的高密度电法剖面显示, 该断层的电阻率异常特征清晰, 其上断点埋深可达20—30 m, 较浅层人工地震探测所揭示的断层上断点埋深(50 m)更浅, 结合地层年代资料推测该断裂的最新活动时期为早—中更新世． 探测结果表明: 高密度电法剖面清晰地显示了断层在浅部松散层的延伸, 适用于日喀则地区的隐伏断层探测; 相较于浅层人工地震探测, 该方法对浅部松散层的探测具有明显优势, 一定条件下能够更好地揭示断层上断点埋深, 可与浅层人工地震探测形成互补． 需要指出的是, 在应用中需重视测区水文地质及地层发育情况对探测的影响．      

综合物探技术在堤坝隐患探测中的应用

本文给出笔者参与的浙江温岭东浦新塘隐患探测、椒江外沙海塘渗漏探测、缙云白马水库和玉环里墩水库坝体渗漏探测等几个实例,详细说明了综合物探技术在堤坝隐患探测中的作用、优点和需要注意的问题,具体分析了堤坝隐患的不同特征和综合物探方法的特点.工作实践表明：高密度电阻率法、探地雷达法、瑞雷波法及其综合运用能有效探测坝体缺陷,查明渗漏原因.同时,由于每种方法各有应用前提及局限性,故应针对可能的隐患特征,选用合适的探测方法及其组合,并努力提高信号采集和处理解释水平.     

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

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