共查询到20条相似文献，搜索用时 46 毫秒
1.
Bispherical coordinates are used to derive an exact mathematical solution for the potential field generated by direct current electric conduction in an earth model consisting of two spherical inclusions in a uniform whole-space. The solution takes the form of a spherical harmonic expansion in bispherical coordinates; coefficients in the expansion are obtained by solving sets of linear equations. Rapid forward modelling of numerous interesting situations in d.c. resistivity prospecting is facilitated by the generality and computational efficiency inherent to this new solution. For example, the accuracy of image (or superposition) methods for calculating potential solutions can be quantified. Similarly, the ability of d.c. conduction methods to resolve two distinct bounded bodies in three-dimensional space can be examined by repeatedly calculating the secondary potential or apparent resistivity response of an earth model as a selected parameter is varied. Synthetic mise à la masse, crosshole, or areal potential data sets can be generated for subsequent use in inversion studies. Improvements in solution technique derived here also apply to a simpler model consisting of a single sphere buried in a half-space. 相似文献
2.
J. Sun 《Geophysical Prospecting》2005,53(3):411-421
We use the integral equation for a d.c. electric field, published in the literature, to introduce the concept of the electrical reflectivity tensor into d.c. electric field modelling. It is shown that in d.c. electric field modelling, the electric reflectivity tensor can be obtained in exactly the same way as in electromagnetic modelling. As a result, for a d.c. electric field, the quasi‐linear and the quasi‐analytical approximations, as well as the quasi‐analytical series, can be constructed in exactly the same way as in electromagnetic modelling. If the primary field is uniform, and if the anomalous body is a uniform circular cylinder or a uniform sphere, the reflectivity tensor is zero order (constant), relating to the free surface charge density. Thus, for some homogeneous bodies that have simple shapes and are embedded in a uniform primary field, the electrical reflectivity tensor is not only a mathematical mechanism for obtaining approximate solutions, but also a physical reality. Indeed, the free surface charge density is defined as the change of the electric displacement vector across the boundary surface under consideration. If the primary field is caused by a point source, and if the anomalous body is a uniform sphere, the reflectivity tensor is second order, varying slowly within the sphere. The relationship to the free surface charge density can be established only when both the reflectivity tensor and the free surface charge density are approximated by the first terms of their series solutions. If the point source is far from the centre of the sphere, the corresponding reflectivity tensor reduces to zero order, and is independent of the observation position within the sphere, i.e. it is a constant. Therefore, the basic idea of the quasi‐analytical approximation, i.e. taking the reflectivity tensor outside the integral operator, is justified in the case considered here. 相似文献
3.
In any numerical solution of the DC resistivity experiment, care must be taken to deal with strong heterogeneity of electrical conductivity. In order to examine the importance of conductivity contrasts, we develop a scattering decomposition of the DC resistivity equation in the sparse differential domain as opposed to the traditional dense integral formulation of scattering‐type equations. We remove the singularity in the differential scattered series via separation of primary and secondary conductivity, thereby avoiding the need to address the singularity in a Green's function. The differential scattering series is observed to diverge for large conductivity contrasts and to converge for small contrasts. We derive a convergence criterion, in terms of matrix norms for the weak‐form finite‐volume equations, that accounts for both the magnitude and distribution of heterogeneity of electrical conductivity. We demonstrate the relationship between the differential scattering series and the Fréchet derivative of the electrical potential with respect to electrical conductivity, and we show how the development may be applied to the inverse problem. For linearization associated with the Fréchet derivative to be valid, the perturbation in electrical conductivity must be small as defined by the convergence of the scattered series. The differential scattering formulation also provides an efficient tool for gaining insight into charge accumulation across contrasts in electrical conductivity, and we present a derivation that equates accumulated surface charge density to the source of scattered potential. 相似文献
4.
S. V. Anisimov S. V. Galichenko K. V. Aphinogenov A. A. Prokhorchuk 《Izvestiya Physics of the Solid Earth》2018,54(5):764-774
The physical mechanisms determining the variability of the vertical profiles of electrical conductivity, space charge density, and electric field in the undisturbed midlatitude lower atmosphere are discussed. The influence of the global and local mesoscale processes on the variability of electrical conductivity and the main component of the atmospheric electric field is estimated. The sunrise effect is studied, estimates are obtained for the charge accumulation rate in the column of the lower atmosphere and the corresponding growth rate of the field strength close to the ground. It is shown that the increase in the average charge density is mainly due to the breakdown of the stable stratification of the atmospheric boundary layer and transformation of the vertical profile of electrical conductivity following the convective mixing of a radon and its daughter products. 相似文献
5.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
6.
首先给出柱坐标系中电导率分段线性变化的水平层的点源电场的二维边值问题,然后用变分法将边值问题转变为变分问题。用有限单元法解变分问题,将区域剖分成矩形单元,在单元中进行双线性函数插值,将变分方程化为线性代数方程组。解方程组,得各节点的电位值,由此可计算地表的视电阻率。 算例表明,本方法计算结果与精确解十分符合。本文还举了一个定量分析视电阻率年变化的例子。 本方法占用计算机内存约100K数量级。在MV/6000超小型计算机上计算一条电测深曲线的时间为几十秒钟。 相似文献
7.
E.H. ELORANTA 《Geophysical Prospecting》1986,34(6):856-872
An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite. 相似文献
8.
A. A. KAUFMAN 《Geophysical Prospecting》1985,33(2):171-184
The main features of the distribution of volume and surface charges in a conducting medium can be described separately for direct and alternating electromagnetic fields. The density of charges depends on the conductivity of a medium and on the electrical field. The relation is particularly simple for the quasi-stationary field, i.e., when the influence of displacement currents is negligible. Conditions are formulated under which electrical charges arise in a conducting medium: electrical charges are shown to exist for direct and quasi-stationary fields when there is a component of electric field parallel to the gradient of conductivity. The density of these charges is proportional to the applied electric field. 相似文献
9.
Following a previous paper in which the principles of a 3D ground-surface tomographic processing of self-potential data were established, we extend the method to active source geoelectric surveying. The main purpose of the new tomographic approach is to obtain a physical image reconstruction of the induced electric charges distributed over buried resistivity discontinuities. The information is produced in a probabilistic sense, as the mathematical formulation underlying the method treats only the intrinsic physical nature of the generated electric field underground and the method of its ground-surface detection, independently of the geometry of the unknown structures. In practice, a 3D tomography is realized by cross-correlating a set of distributed electric-field ground-surface data with a scanning function, representing a unit positive point charge located anywhere in the lower half-space. The resolution of the method is tested on the synthetic response of a 3D structural simulation of an archaeological target, consisting of an infinitely resistive prismatic body immersed in a half-space, including surface inhomogeneities and layering. Finally, the field response of a 3D structure consisting of a hypogeal dromos-chamber tomb inside the Sabine Necropolis at Colle del Forno, close to Rome, is presented and discussed. 相似文献
10.
电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解. 相似文献
11.
煤矿井下矿井瞬变电磁法(MTEM)探测中,电磁场呈全空间分布,全空间瞬变电磁反演是复杂的非线性问题,目前反演计算中全空间响应主要由半空间响应乘以全空间响应系数来得到,导致反演结果中顶板和底板异常(或前方和后方异常)叠加在一起难以分离,造成分辨率下降.论文提出采用粒子群优化算法(PSO)进行全空间MTEM反演,通过理论分析,在常规的粒子群算法基础上提出了一种新的进化公式改进策略,提高了粒子群算法的寻优能力.基于全空间瞬变电磁场理论,编写了粒子群算法反演程序,进行全空间条件下五层含巷道的复杂模型的反演计算.结合某矿井巷道顶板、底板岩层及断层含水性的探测实例,对实测数据进行反演计算和解释,探测结果得到钻探证实.研究表明,改进的粒子群优化算法对理论模型和实际资料的反演拟合程度较高,实现了矿井顶板、底板视电阻率异常的分离,提高了全空间瞬变电磁勘探资料的解释精度和分辨率. 相似文献
12.
New Methods for Modeling Laterolog Resistivity Corrections 总被引：1，自引：1，他引：0
Jadwiga Anna Jarzyna Adam Cichy Dezső Drahos Attila Galsa Maria Joanna Bała Andrzej Ossowski 《Acta Geophysica》2016,64(2):417-442
The paper presents methods for laterolog response modeling. In Coulomb’s charges method, Laplace’s equation is solved for the electric field distribution in rock medium with internal boundaries between different resistivity layers. There, the boundary problem is reduced to Fred-holm integral equation of the second kind. The second method uses a finite element array to model apparent resistivity from laterolog. The task is treated as DC problem and the Laplace equation is solved numerically. The presented methods were applied to borehole data covering a typical stratigraphie section of the Fore-Sudetic Monocline in southwestern Poland. Apparent resistivity was calculated using the Coulomb’s charges method and alternatively modeled using a finite element method which gave similar results. Then, a series of linear corrections for borehole, shoulder bed, and filtration effects for apparent resistivity obtained by the Coulomb’s charges method demonstrated the feasibility of calculating true resistivity of virgin and invaded zones. The proposed methods provide a flexible solution in modeling which can be adapted to other logs. 相似文献
13.
为研究电测井和油气层的关系,根据电场理论,将砂岩电导率定义为地层水电导率、含量、相对电场强度的积.提出电介质颗粒面电荷与电偶极子等效的认识,进一步得出砂岩的电场方程.解砂岩电场方程得出相对电场强度表达式,将其带入砂岩电导率定义得电测井和油气层关系的一个理论公式,并用实验资料验证了公式的有效性. 相似文献
14.
Cross-well electrical measurement as known in the oil industry is a method for determining the electrical conductivity distribution between boreholes from the electrostatic field measurements in the boreholes. We discuss the reconstruction of the conductivity distribution of a three-dimensional domain. The measured secondary electric potential field is represented in terms of an integral equation for the vector electric field. This integral equation is taken as the starting point to develop a non-linear inversion method, the so-called contrast source inversion (CSI) method. The CSI method considers the inverse scattering problem as an inverse source problem in which the unknown contrast source (the product of the total electric field and the conductivity contrast) in the object domain is reconstructed by minimizing the object and data error using a conjugate-gradient step, after which the conductivity contrast is updated by minimizing only the error in the object. This method has been tested on a number of numerical examples using the synthetic 'measured' data with and without noise. Numerical tests indicate that the inversion method yields a reasonably good reconstruction result, and is fairly insensitive to added random noise. 相似文献
15.
We report electric potential gradient measurements carried out at Sakurajima volcano in Japan during: (1) explosions which generated ash plumes, (2) steam explosions which produced plumes of condensing gases, and (3) periods of ashfall and plume-induced acid rainfall. Sequential positive and negative deviations occurred during explosions which generated ash plumes. However, no deflections from background were found during steam explosions. During periods of ashfall negative electric potential gradients were observed, while positive potential gradients occurred during fallout of plume-induced acid rain from the same eruption. These results suggest that a dipole arrangement of charge develops within plumes such that positive charges dominate in the volcanic gas-rich top and negative charges in the following ash-rich part of the plume. The charge polarity may be reversed for other volcanoes (Hatakeyama and Uchikawa 1952). We suggest that charge is generated by fracto-emission (Donaldson et al. 1988) processes probably during magma fragmentation within the vent, rather than by frictional effects within the plume. 相似文献
16.
A three‐dimensional (3D) electrical resistivity modelling code is developed to interpret surface and subsurface data. Based on the integral equation, it calculates the charge density caused by conductivity gradients at each interface of the mesh, allowing the estimation of the potential everywhere without the need to interpolate between nodes. Modelling generates a huge matrix, made up of Green's functions, which is stored by using the method of pyramidal compression. The potential is compared with the analytical and the numerical solutions obtained by finite‐difference codes for two models: the two‐layer case and the vertical contact case. The integral method is more accurate around the source point and at the limits of the domain for the potential calculation using a pole‐pole array. A technique is proposed to calculate the sensitivity (Jacobian) and Hessian matrices in 3D. The sensitivity is based on the derivative with respect to the block conductivity of the potential computed using the integral equation; it is only necessary to compute the electrical field at the source location. A direct extension of this technique allows the determination of the second derivatives. The technique is compared with the analytical solutions and with the calculation of the sensitivity according to the method using the inner product of the current densities calculated at the source and receiver points. Results are very accurate when the Green's function that includes the source image is used. The calculation of the three components of the electric field on the interfaces of the mesh is carried out simultaneously and quickly, using matrix compression. 相似文献
17.
Domenico Patella 《Geophysical Prospecting》1997,45(4):653-681
A new approach to self-potential (SP) data interpretation for the recognition of a buried causative SP source system is presented. The general model considered is characterized by the presence of primary electric sources or sinks, located within any complex resistivity structure with a flat air-earth boundary. First, using physical considerations of the nature of the electric potential generated by any arbitrary distribution of primary source charges and the related secondary induced charges over the buried resistivity discontinuity planes, a general formula is derived for the potential and the electric field component along any fixed direction on the ground surface. The total effect is written as a sum of elementary contributions, all of the same simple mathematical form. It is then demonstrated that the total electric power associated with the standing natural electric field component can be written in the space domain as a sum of cross-correlation integrals between the observed component of the total electric field and the component of the field due to each single constitutive elementary charge. By means of the cross-correlation bounding inequality, the concept of a scanning function is introduced as the key to the new interpretation procedure. In the space domain, the scanning function is the unit strength electric field component generated by an elementary positive charge. Next, the concept of charge occurrence probability is introduced as a suitable function for the tomographic imaging of the charge distribution geometry underground. This function is defined as the cross-correlation product of the total observed electric field component and the scanning function, divided by the square root of the product of the respective variances. Using this physical scheme, the tomographic procedure is described. It consists of scanning the section, through any SP survey profile, by the unit strength elementary charge, which is given a regular grid of space coordinates within the section, at each point of which the charge occurrence probability function is calculated. The complete set of calculated grid values can be used to draw contour lines in order to single out the zones of highest probability of concentrations of polarized, primary and secondary electric charges. An extension to the wavenumber domain and to three-dimensional tomography is also presented and discussed. A few simple synthetic examples are given to demonstrate the resolution power of the new SP inversion procedure. 相似文献
18.
L. ESKOLA 《Geophysical Prospecting》1979,27(3):616-627
A numerical method is proposed for solving the problem of steady current flow. The electrodynamic model is replaced by the equivalent stationary charge distribution obtained by Poisson's analysis, in which the surface integral equation for field intensity is reduced to a set of simultaneous linear algebraic equations by means of the method of sub-areas. The solution of the set allows the calculation of an approximation for the charge density distribution on the discontinuity surfaces of conductivity. The method is valid for complex conductivities, whereby the apparent phase shift of IP can be calculated from the complex potential or field intensity. The phase shift anomaly calculated as an application is very similar to the corresponding frequency effect anomaly. The method allows the calculation of the mise-à-la-masse effect as a solution to a potential problem, in which the primary current electrode is located within the body to be surveyed. 相似文献
19.
M. OKABE 《Geophysical Prospecting》1982,30(5):653-672
A Fredholm integral equation of the first kind with respect to the surface charge density is presented via a weighted residual formulation in the standard isotropic problem. The surface charge densities are numerically obtained by solving the well-known Fredholm integral equation of the second kind, where the new equation can be regarded as a constraint. The accuracy of boundary element solutions is examined in connection with the violation of such a constraint and the “modified reciprocal averaging techniques” are proposed. 相似文献
20.
利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。 相似文献