共查询到20条相似文献，搜索用时 625 毫秒
1.
2.
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探. 相似文献
3.
Raghu K. Chunduru Mrinal K. Sen Paul L. Stoffa R. Nagendra 《Geophysical Prospecting》1995,43(8):979-1003
The inversion of resistivity profiling data involves estimation of the spatial distribution of resistivities and thicknesses of rock layers from the apparent resistivity data values measured in the field as a function of electrode separation. The drawbacks of using traditional curve-matching techniques to solve this inverse problem have been overcome by iterative linear techniques but these require good starting models even if the shape of the causative body is asssumed known. In spite of the recent developments in inversion techniques, no robust method exists for the inversion of resistivity profiling data for the simple model of dikes and spheres which are the classical models of geophysical prospecting. We apply three different non-linear inversion schemes to invert synthetic resistivity profiling data for the classical models embedded in a uniform matrix of contrasting resistivity. The three non-linear algorithms used are called the Metropolis simulated annealing (SA), very fast simulated annealing (VFSA) and a genetic algorithm (GA). We compare the performance of the three algorithms using synthetic data for an outcropping vertical dike model. Although all three methods were successful in obtaining optimal solutions for arbitrary starting models, VFSA proved to be computationally the most efficient. 相似文献
4.
为充分利用微动信号中基阶和高阶模式瑞雷波,本文研究了基于多阶瑞雷波SPAC系数直接反演的方法.该方法首先基于地层介质响应计算多阶瑞雷波的能量占比,考虑实际观测台阵有限台站个数对SPAC系数影响,正演计算多阶瑞雷波SPAC系数,再采用快速模拟退火算法对其反演以获得地下介质横波速度结构.在此基础上,本文通过数值模拟验证该方法的可靠性,分别选取三种典型地质模型,基于模式叠加算法合成理论微动信号,采用本文方法计算其理论多阶瑞雷波SPAC系数并反演,给出反演结果与真实模型对比.我们将该方法应用于上海中心城区的地质调查中,通过与钻探结果对比,进一步验证该方法的有效性.本文理论与实际应用研究表明,基于多阶瑞雷波SPAC系数直接反演的微动探测方法有助于提高反演结果的可靠性,尤其对含软硬夹层的复杂地层介质,可提高探测精度. 相似文献
5.
A code for 3-D resistivity modelling and inversion of vertical electrical soundings has been developed based on the finite-element technique and regularisation method. Synthetic data were used to test the effectiveness of the code and to examine the resolving capability of the Schlumberger array in investigating 3-D resistivity distributions. The code was applied to experimental data set constituted by 35 Schlumberger soundings collected near the Cairo city in order to study the subsurface resistivity distribution. The results have shown that valuable imaging of the subsurface resistivity distribution can be constructed even when the vertical electrical soundings are acquired in a sparse field data set. 相似文献
6.
Yasala Srinivas Stanley A. Raj Oliver D. Hudson Durairaj Muthuraj Nainarpandian Chandrasekar 《Studia Geophysica et Geodaetica》2013,57(3):520-534
Soft computing tools play a vital role in fixing certain non-linear problems related to the earth. More specifically, digging out the mysteries of subsurface of the earth, the nonlinearity can be converging to assemble an approximate solution which resembles the real characteristics of the earth. Adaptive Neuro Fuzzy Inference System (ANFIS) tool is one of the best soft computing tools to estimate the complex data analysis. ANFIS was applied to estimate the subsurface parameters of earth using the Vertical Electrical Sounding (VES) data. Classifying the lithology based on the resistivity values by ANFIS is employed here in this paper. As the resistivity of each formation varies in range of values, ANFIS tool thus approximates the subsurface features based on effective training. In this study, ANFIS performance was checked with training data, and successively it has been tested with the field data. Optimized ANFIS algorithm provides the necessary tool for predicting the non-linear subsurface features. The best training performance of this soft computing tool efficiently predicts the subsurface lithology. Also the interpreted results show the true resistivity and thickness of the subsurface layers of the earth. The proposed technique was represented in Graphical User Interface (GUI), and the lithological variables are predicted in texture format and linguistic variables. 相似文献
7.
研究和提出了一种坑道直流电阻率超前聚焦探测新方法.该方法通过在掌子面上设立不同功能的环状电极组,使一次场电流具有像探照灯一样的聚焦功能,可有效探测坑道掘进前方不良地质体的存在,达到超前预报的目的.本文首先介绍了直流电阻率超前聚焦探测的基本原理,并给出了两种聚焦观测方式.为了验证所提方法的可行性,通过轴对称电性介质二维异常电位有限元数值模拟方法,对聚焦观测条件下几例坑道模型进行了模拟计算.结果表明,该方法对异常反映明显,可进行即时解释,建议在坑道超前预报中逐步推广使用. 相似文献
8.
The interpretation of electrical sounding data for a subsurface with monotonic continuous variation of the resistivity with depth is becoming increasingly necessary. The contribution of this article is the derivation of the solution for the Wenner and the Schlumberger apparent resistivity functions for a resistivity varying as a real power of a linear positive function of the depth. The interpretation of sounding data in these cases can be used to estimate the variation of the porosity or the salt content of the pore water with depth. 相似文献
9.
Interpretation of a single geophysical data set is not sufficient to get complete subsurface information. Cooperative or joint inversion of geophysical data sets is the preferred method for most case studies. In the present study, we present the results of the cooperative inversion approach of direct current resistivity and gravity data. The algorithm uses fuzzy c-means clustering to determine the petrophysical relationship between density and resistivity to obtain the similarity. Synthetic data set has demonstrated that the cooperative inversion approach can produce more reliable and better resistivity and density models of the subsurface as compared to those obtained through individual inversions. To utilize the presented cooperative inversion algorithm, the number of geologic units (number of clusters) in the study region must be known a priori. As a field study, the cooperative inversion approach was used to identify the extension of uranium-bearing target rock around the Beldih open cast mine. We noted the inconsistencies in both resistivity and density models obtained from the individual inversions. However, the presented cooperative inversion approach was able to produce similar resistivity and density models while maintaining the same error level of the respective individual inversions. We have considered four geologic units in the presented cooperative inversion as a field case study. We have also compared our cooperative results with drilled borehole and found to be a reliable tool to differentiate between the target rocks (kaolinite and quartz–magnetite–apatite rocks) and the ultramafic rock (host rock quartzite/alkaline granite). However, this study is subject to certain limitations such as the inability to differentiate between closely spaced kaolinite and quartz–magnetite–apatite rocks. 相似文献
10.
Surface geophysical surveys provide an effective way to image the subsurface and the ground water zone without a large number of observation wells. DC resistivity sounding generally identifies the subsurface formations-the aquifer zone as well as the formations saturated with saline/brackish water. However, the method has serious ambiguities in distinguishing the geological formations of similar resistivities such as saline sand and saline clay, or water quality such as fresh or saline, in a low resistivity formation. In order to minimize the ambiguity and ascertain the efficacy of data integration techniques in ground water and saline contamination studies, a combined geophysical survey and periodic chemical analysis of ground water were carried out employing DC resistivity profiling, resistivity sounding, and shallow seismic refraction methods. By constraining resistivity interpretation with inputs from seismic refraction and chemical analysis, the data integration study proved to be a powerful method for identification of the subsurface formations, ground water zones, the subsurface saline/brackish water zones, and the probable mode and cause of saline water intrusion in an inland aquifer. A case study presented here illustrates these principles. Resistivity sounding alone had earlier failed to identify the different formations in the saline environment. Data integration and resistivity interpretation constrained by water quality analysis led to a new concept of minimum resistivity for ground water-bearing zones, which is the optimum value of resistivity of a subsurface formation in an area below which ground water contained in it is saline/brackish and unsuitable for drinking. 相似文献
11.
Ahmad Neyamadpour 《Studia Geophysica et Geodaetica》2011,55(4):755-769
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. 相似文献
12.
Use of VFSA for resolution, sensitivity and uncertainty analysis in 1D DC resistivity and IP inversion 总被引：1，自引：0，他引：1
We present results from the resolution and sensitivity analysis of 1D DC resistivity and IP sounding data using a non-linear inversion. The inversion scheme uses a theoretically correct Metropolis–Gibbs' sampling technique and an approximate method using numerous models sampled by a global optimization algorithm called very fast simulated annealing (VFSA). VFSA has recently been found to be computationally efficient in several geophysical parameter estimation problems. Unlike conventional simulated annealing (SA), in VFSA the perturbations are generated from the model parameters according to a Cauchy-like distribution whose shape changes with each iteration. This results in an algorithm that converges much faster than a standard SA. In the course of finding the optimal solution, VFSA samples several models from the search space. All these models can be used to obtain estimates of uncertainty in the derived solution. This method makes no assumptions about the shape of an a posteriori probability density function in the model space. Here, we carry out a VFSA-based sensitivity analysis with several synthetic and field sounding data sets for resistivity and IP. The resolution capability of the VFSA algorithm as seen from the sensitivity analysis is satisfactory. The interpretation of VES and IP sounding data by VFSA, incorporating resolution, sensitivity and uncertainty of layer parameters, would generally be more useful than the conventional best-fit techniques. 相似文献
13.
由于可控源音频大地电磁法（CSAMT）采用人工场源在大大增加信号强度的同时也带来了在近区产生非平面波效应的问题,其表现特征为在双对数坐标中卡尼亚电阻率曲线呈45°上升,即所谓的近场效应.本文首先提出了基于牛顿迭代法求解视电阻率方程的近场效应校正方法,通过对二层、三层理论模型试算验证了该方法的有效性;为进一步增强对噪声的压制能力,本文进一步引入了遗传反演的思想,将全频段误差最小作为目标函数,采用自适应正则化方法引入模型参数梯度最小作为稳定因子,由遗传算法求解得到校正后视电阻率,通过对理论模拟数据加随机噪声的校正结果,表明该方法在实现近场校正的同时能剔除噪声干扰,得到的结果与理论模拟曲线吻合度较高.通过对内蒙古曹四夭钼矿区实测CSAMT数据处理结果表明,牛顿法和遗传算法均能明显校正CSAMT因近场效应引起的假高阻异常,能更好地反映地电结构特征. 相似文献
14.
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. 相似文献
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. 相似文献
15.
David Beamish 《Journal of Applied Geophysics》1994,32(4)
Very low frequency electromagnetic (EM) methods using VLF transmitters have found many applications in subsurface geophysical investigations. Surface measurements involving both the vertical component of the magnetic field (VLF-EM or VLF-Z) and of the apparent resistivity (VLF-R) are increasingly common. Although extensive VLF data sets have been successfully used for mapping purposes, modelling and interpretation techniques which asess the third (i.e. depth) dimension appear limited.Given a profile of VLF-R measurements the main purpose of the present study is to demonstrate an automatic method for the construction of a resistivity cross-section. The technique used is one of a new generation of regularised inversion methods. These techniques attempt to overcome the problem of equivalence/non-uniqueness in EM sounding data by constructing the resistivity distribution with the minimum amount of structure that fits the data.VLF data represent a special case of plane-wave EM sounding in that they conform, in practice, to a single-frequency technique. This fact imposes a limitation in the amount of vertical resolution that we can expect using such data. In the case of two-dimensional modelling and inversion, resolution through the cross-section is a resultant attribute from both vertical and lateral resistivity gradients within the subsurface. In order to provide insight into the practical application of regularised inversion techniques to VLF data, both synthetic and field examples are considered. Both sets of examples are primarily concerned with VLF data applied to near-surface fault mapping where the main aim is to assess the location, dip and depth extent of conductive subsurface features. 相似文献
16.
Automatic extraction of geometrical characteristics hidden in ground‐penetrating radar sectional images using simultaneous perturbation artificial bee colony algorithm 
下载免费PDF全文
下载免费PDF全文 Ground‐penetrating radar is one of the most effective methods of detecting shallow buried objects. Ground‐penetrating radar radargram is a vertical map of the radar pulse reflections that returns from subsurface objects, and in the case of cylindrical objects, it would be a hyperbola. In order to get clear and accurate information on the presence, location, and geometry of the buried objects, the radargrams need to be interpreted. Interpretation of the results is a time‐consuming task and needs an expert with vast knowledge. Development of an automatic interpretation method of B‐scan ground‐penetrating radar images would be an effective and efficient solution to this problem. A novel automatic interpretation method of ground‐penetrating radar images, based on simultaneous perturbation artificial bee colony algorithm using tournament selection strategy, simultaneous perturbation stochastic approximation method, and new search equations, is introduced in this paper. The proposed algorithm is used to extract geometrical parameters, i.e. depth, location, and radius, of buried cylindrical objects in order to assess its accuracy. Synthetic data, simulated using GprMax2D forward modelling program, and real data, surveyed in the campus of Isfahan University of Technology, are used in the assessment. The performance of the proposed method in detecting synthetic hyperbolas is compared with that of the original artificial bee colony algorithm, genetic algorithm, and modified Hough transform. The results show superiority of the proposed algorithm, in detecting synthetic hyperbolas. Furthermore, the performance of the proposed method in estimating depth and radius of pipes in real ground‐penetrating radar images is compared with that of the modified Hough transform. The results indicate higher accuracy of the proposed method in estimating geometrical parameters of the buried cylindrical objects. 相似文献
17.
《Journal of Applied Geophysics》1999,41(1):31-47
The Idaho National Laboratory (INEL) Cold Test Pit (CTP) has been carefully constructed to simulate buried hazardous waste sites. An induced polarization (IP) survey of the CTP shows a very strong polarization and a modest resistivity response associated with the simulated waste. A three-dimensional (3-D) inversion algorithm based on the simultaneous iterative reconstruction technique (SIRT) and finite difference forward modelling has been applied to generate a subsurface model of complex resistivity. The lateral extents of the waste zone are well resolved. Limited depth extent is recognized, but the bottom of the waste appears too deep. With a modelling experiment, the intrinsic polarizability of the waste material is determined. Since IP is a technique for detection of diffuse occurrences of metallic material, this method holds promise as a method to distinguish buried waste from conductive soil material. 相似文献
18.
实际地层可能出现2种电阻率的宏观各向异性:微观的统计平均效应和裂缝走向等因素的构造效应.三分量感应测井方法可能是确定地层电阻率各向异性的最好方法,但是,该仪器目前在中国还没有应用.已有的数据一般是梯度和双感应曲线,梯度和双感应曲线单独使用都难以揭示各向异性,但是,将二者联合使用则有可能.在反演模型中考虑了电阻率各向异性的存在,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法.通过直井水平层的二维人工模型表明,用感应方法可确定地层的水平电阻率和厚度.而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值,梯度方法的视厚度是各向异性参数和真实厚度的乘积.因而不能用梯度方法单独确定各向异性参数和真实厚度2个参数中的任何1个.然而,联合梯度和双感应的方法则可能确定3个参数:各向异性参数,水平电阻率和地层厚度.人工数据的例子表明了这种联合的可行性. 相似文献
19.
多通道瞬变电磁法由于采用伪随机编码发射、最佳匹配相干解码技术,所以相比于传统时间域电磁方法勘探深度大且抗干扰能力强,是目前地面资源勘查领域的研究热点.针对内蒙兴安盟某铅锌银矿埋藏深度大等特点,首先根据矿区资料建立一维以及二维模型,数值计算脉冲响应并分析其特征,说明多通道瞬变电磁对于大埋深低阻异常体有良好的反映,随后通过多通道瞬变电磁技术对矿区进行实例勘查,数据处理过程中利用了反卷积的方法提取大地脉冲响应,根据轴向大地脉冲响应的峰值时刻计算视电阻率,通过对视电阻率进行全波形反演得到解释剖面,与测区钻孔资料对比分析证明,多通道瞬变电磁法是反映深部低阻异常体的有效手段. 相似文献
20.
与可控源音频大地电磁（CSAMT）相比,广域电磁法通过采用全区视电阻率定义,突破了卡尼亚视电阻率所需的远区条件限制,极大拓展了可控源电磁观测区域和探测深度.考虑到电偶源激发场的三维特征以及地下复杂三维结构,为提高广域电磁数据解释精度,本文实现了基于二次耦合势的广域电磁法三维正演计算.该算法利用Helmholtz定理将麦克斯韦方程转化为库伦规范下的磁矢势和电标势耦合方程,有效改善了离散所得大型线性方程组的谱性质,并通过强加散度条件来消除电场伪解的影响.此外,采用散射场方法,其中一次场使用准解析法求解,二次场使用有限体积法求解,克服了局部激发场源奇异性问题.通过与一维层状模型下电偶源产生的电磁场准解析解对比,验证了本文算法的正确性.在此基础上,利用本文的正演算法对比分析了广域电磁法与CSAMT对典型三维目标体的探测能力,结果表明在相同的观测条件下,广域电磁法能够更准确地反映地下目标体信息,拥有更优的分辨能力. 相似文献