共查询到20条相似文献，搜索用时 31 毫秒
1.
Following the probability tomography principles previously introduced to image the sources of electric and electromagnetic anomalies, we demonstrate that a similar approach can be used to analyse gravity data. First, we give a coherent derivation of the Bouguer anomaly concept as a Newtonian-type integral for an arbitrary mass distribution buried below a non-flat topography. A discretized solution of this integral is then derived as a sum of elementary contributions, which are cross-correlated with the gravity data function in the expression for the total power associated with the Bouguer anomaly. To image the mass distribution underground we introduce a mass contrast occurrence probability function using the cross-correlation product of the observed Bouguer anomaly and the synthetic field due to an elementary mass contrast source. The tomographic procedure consists of scanning the subsurface with the elementary source and calculating the occurrence probability function at the nodes of a regular grid. The complete set of grid values is used to highlight the zones of highest probability of mass contrast concentrations. Some synthetic and field examples demonstrate the reliability and resolution of the new gravity tomographic approach. 相似文献
2.
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. 相似文献
3.
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. 相似文献
4.
Numerical computations using the integral equation method are presented for resistivity and IP responses due to arbitrarily shaped 3-dimensional bodies in a layered earth. The unknown surface charge density distribution is expressed as the solution of Fredholm's integral equation of the second kind. Use of moment method (with pulse basis function and point-collocation) yields the matrix equations for the unknowns. The contributions to Green's function are solved (a) analytically for the primary and (b) by convolution for the secondary contributions resulting in a fast algorithm. The further step of computing potential, apparent resistivity, chargeability etc., for any electrode system, is straightforward. Our results show a good agreement with those from finite difference methods and physical tank experiments. The CPU time is only 138 s on a super-minicomputer for an apparent resistivity pseudo-section, even with 96 elementary cells as used for discretization. A large number of models for different geological situations were studied; some are presented here. 相似文献
5.
基于地下电流场的积分公式,离散化的二次电流场被分解为电性不连续界面上的一系列点电荷电场的代数和.单位正点电荷电场被引入作为空间扫描函数（SDS）,积累电荷出现的概率（COP）函数定义为二次电流场与SDS的互相关.为对概率成像结果进行定量分析解释,提出了规范的积累电荷出现的概率（NCOP）函数.通过应用有限元算法对2D地电模型进行二次电流场合成,实现了二次电流场的多次叠加概率成像.结果表明对均匀半空间中赋存地质异常体的电性结构,概率成像方法对地下异常体的空间位置有较好的指示作用. 相似文献
6.
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. 相似文献
7.
When an electric current is introduced to the earth, it sets up a distribution of charges both on and beneath the earth's surface. These charges give rise to the anomalous potential measured in the d. c. resistivity experiment. We investigate different aspects of charge accumulation and its fundamental role in d. c. experiments. The basic equations and boundary conditions for the d. c. problem are first presented with emphasis on the terms involving accumulated charges which occur wherever there is a non-zero component of electric field parallel to the gradient of conductivity. In the case of a polarizable medium, the polarization charges are also present due to the applied electric field, yet they do not change the final field distribution. We investigate the precise role of the permittivity of the medium. The charge buildup alters the electric fields and causes the refraction of current lines; this results in the well-known phenomenon of current channelling. We demonstrate this by using charge density to derive the refraction formula at a boundary. An integral equation for charge density is presented for whole-space models where the upper half-space is treated as an in-homogeneity with zero conductivity. The integral equation provides a tool with which the charge accumulation can be examined quantitatively and employed in the practical forward modelling. With the aid of this equation, we investigate the effect of accumulated charges on the earth's surface and show the equivalence between the half-space and whole-space formulations of the problem. Two analytic examples are presented to illustrate the charge accumulation and its role in the d. c. problem. We investigate the relationship between the solution for the potential via the image method and via the charge density. We show that the essence of the image method solution to the potential problem is to derive a set of fictitious sources which produce the same potential as does the true charge distribution. It follows that the image method is viable only when the conductivity structure is such that the effect of the accumulated charge can be represented by a set of point images. As numerical examples, we evaluate quantitatively the charge density on the earth's surface that arises because of topography and the charge density on a buried conductive prism. By these means, we demonstrate the use of the boundary element technique and charge density in d. c. forward modelling problems. 相似文献
8.
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. 相似文献
9.
J. DOHERTY 《Geophysical Prospecting》1988,36(6):644-668
The theory by which the Surface Integral Equation method may be applied to the solution of electromagnetic transmission boundary value problems is presented. For a 3D target of arbitrary electrical property contrast with its host medium excited by an arbitrary time-harmonic source, two integral equations are derived which need to be simultaneously solved for tangential electric and magnetic source density on the target's surface. If the target is 2D, though still excited by an arbitrary source (the 2½ D case), the problem is best solved in the transform domain for a number of different wavenumbers in the target's strike direction. Then a set of four simultaneous scalar integral equations needs to be solved for the components of the surface source density transforms in the target's strike direction and in the direction of the tangent vector to the target's cross-sectional contour. Examples are presented in which the 2½D problem is solved numerically using the method of moments with piecewise linear basis functions. Although the results generally compare well with analytical solutions, or solutions obtained numerically by other means, errors appear in the calculation of the real response of these targets to excitation by a magnetic dipole source at low frequencies. This is attributed to ill-conditioning of the system resulting from a non-unique solution at zero frequency. 相似文献
10.
Summary The potential of the electric field of a stationary current in a two-layered Earth is calculated by applying Green's formula in the case where a three-dimensional inhomogeneity of different conductivity is located in the basement of the layer. It is proved that the potential outside and inside the perturbing body can be calculated from the potential of an electric double-layer distributed on the surface of this body. An integral equation of the Fredholm type is derived for the surface density of the double-layer, together with some of its integral properties. A similar procedure can be applied to computing the magnetic anomalies of three-dimensional magnetized bodies, geothermal anomalies due to three-dimensional inhomogeneities of different heat conductivity, as well as to potential problems of theoretical electrical engineering. 相似文献
11.
对比分析了随机结构动力可靠度计算的三种估计算法.渐进展开法是基于Laplace算法对概率积分进行渐进估计的,此法通过计算最大被积分式值对应点,并将其代入概率积分的渐进估计表达式求解失效概率.由于概率积分的主要贡献来自于最大被积分式值对应点的周围,因此本文的重要抽样法假定重要抽样函数的最大似然值等于最大被积分式值对应点值.极值分布-泰勒展开法首先通过结构随机参数的极值分布函数给出失效概率的表达式,随后利用泰勒展开法对失效概率进行估计,其中采用中心差分法对极值分布函数的梯度进行估算.最后应用三种算法和Monte Carlo法对受高斯白噪声激励作用的单自由度随机结构进行了计算,结果表明三种方法不但运算简便,而且对比Monte Carlo法计算效率有显著提高. 相似文献
12.
常规电测井一般通过测量供电电极的电流和测量电极之间的电位差来计算围岩的视电阻率并划分地层界面,而忽略了对电流参数的利用。本文基于欧姆定律,通过供电电极电流的变化定性分析其周围岩矿石电阻的变化,进而用电流强度变化曲线识别地层属性。理论分析和试验结果表明,供电电极电流大小取决于电极的接触电阻,而视电阻率大小取决两个测量电极之间或供电与测量电极之间介质的电阻率总和,因此,电流比视电阻率对地层的灵敏度要高一些。绘制电流强度测井曲线可以准确划分地层层位,确定地层厚度,且不增加成本。 相似文献
13.
We present a concept of the hybrid finite volume–integral equation technique for solving Maxwell's equation in a quasi-static form. The divergence correction was incorporated to improve the convergence and stability of the governing linear system equations which pose a challenge on the discretization of the curl–curl Helmholtz equation. A staggered finite volume approach is applied for discretizing the system of equations on a structured mesh and solved in a secondary field technique. The bi-conjugate gradient stabilizer was utilized with block incomplete lower-upper factorization preconditioner to solve the system of equation. To obtain the electric and magnetic fields at the receivers, we use the integral Green tensor scheme. We verify the strength of our hybrid technique with benchmark models relative to other numerical algorithms. Importantly, from the tested models, our scheme was in close agreement with the semi-analytical solution. It also revealed that the use of a quasi-analytical boundary condition helps to minimize the runtime for the linear system equation. Furthermore, the integral Green tensor approach to compute at the receivers demonstrates better accuracy compared with the conventional interpolation method. This adopted technique can be applied efficiently to the inversion procedure. 相似文献
14.
A. ROY 《Geophysical Prospecting》1978,26(2):442-463
For any direct current regime, the theorem holds, where φp is the total measured or calculated potential at any point P, φ is the potential distribution known a priori, r is the distance between P and any volume element dV, the gradients are evaluated at the element, and the current sources and sinks have finite dimensions. Thus, each space element behaves as a dipole of moment (1/4π) ?φdV and contributes its share of signal or potential accordingly. By suitable summation or integration, the contribution from any assigned portion of space to the total measured signal can be determined. Except for the chargeability factor m, the formula also establishes Seigel's initial postulate for the time domain induced polarisation theory. The contribution depends on the potential gradient, not the current density, and the integration extends over the entire space. Although an insulating target carries no current, it contributes a signal that is in general larger than normal by virtue of its higher potential gradient, and thus helps in creating an overall positive anomaly or resistivity high. On the other hand, an infinitely conducting target—even though it supports a larger amount of current than normal—contributes nothing to the measured signal as the potential gradient is zero everywhere inside. Thus, by contributing less than normal, a conducting target promotes the creation of what is usually a resistivity low. In all cases, the contributions from the space elements add up exactly to the measured or total calculated value. Some other consequences of the theorem are also discussed, especially in relation to a simple two-layer earth. For instance, the contribution from the upper half-space (air) turns out to be equal to that from the lower (real ground), for all observation points on the ground surface and for any ground configuration. 相似文献
15.
16.
Resistivity anomaly imaging by probability tomography 总被引：10，自引：0，他引：10
Probability tomography is a new concept reflecting the inherently uncertain nature of any geophysical interpretation. The rationale of the new procedure is based on the fact that a measurable anomalous field, representing the response of a buried feature to a physical stimulation, can be approximated by a set of partial anomaly source contributions. These may be given a multiplicity of configurations to generate cumulative responses, which are all compatible with the observed data within the accuracy of measurement. The purpose of the new imaging procedure is the design of an occurrence probability space of elementary anomaly sources, located anywhere inside an explored underground volume. In geoelectrics, the decomposition is made within a regular resistivity lattice, using the Frechet derivatives of the electric potential weighted by resistivity difference coefficients. The typical tomography is a diffuse image of the resistivity difference probability pattern, that is quite different from the usual modelled geometry derived from standard inversion. 相似文献
17.
The Cagniard-de Hoop method is ideally suited to the analysis of wave propagation problems in stratified media. The method applies to the integral transform representation of the solution in the transform variables (s, p) dual of the time and transverse distance. The objective of the method is to make the p-integral take the form of a forward Laplace transform, so that the cascade of the two integrals can be identified as a forward and inverse transform, thereby making the actual integration unnecessary. That is, the exponent (–sw(p)) is set equal to –sτ, with τ varying from some (real) finite time to infinity. As usually presented, the p-integral is deformed onto a contour on which the exponent is real and decreases to –∞ as p tends to infinity. We have found that it is often easier to introduce a complex variable τ for the exponent and carry out the deformation of contour in the complex τ-domain. In the τ-domain the deformation amounts to ‘closing down’ the contour of integration around the real axis while taking due account of singularities off this axis. Typically, the method is applied to an integral that represents one body wave plus other types of waves. In this approach, the saddle point of w(p) that produces the body wave plays a crucial role because it is always a branch point of the integrand in the τ-domain integral. Furthermore, the paths of steepest ascent from the saddle point are always the tails of the Cagniard path along which w(p) →∞. That is, the image of the pair of steepest ascent paths in the p-domain is a double covering of a segment of the Re τ-axis in the τ-domain. The deformed contour in the p-domain will be the only pair of steepest ascent paths unless the original integrand had other singularities in the p-domain between the imaginary axis and this pair of contours. This motivates the definition of a primary p-domain, i.e. the domain between the imaginary axis and the steepest descent paths, and its image in the τ-domain, the primary τ-domain. In terms of these regions, singularities in the primary p-domain have images in the primary τ-domain and the deformation of contour on to the real axis in the τ-domain must include contributions from these singularities. This approach to the Cagniard-de Hoop method represents a return from de Hoop's modification to Cagniard's original method, but with simplifications that make the original method more tractable and straightforward. This approach is also reminiscent of van der Waerden's approach to the method of steepest descents, which starts exactly the same way. Indeed, after the deformation of contour in the τ-domain, the user has the choice of applying asymptotic analysis to the resulting ‘loop’ integrals (Watson's lemma) or continuing to obtain the exact, although usually implicit, time-domain solution by completing the Cagniard-de Hoop analysis. In developing the method we examine the transformation from a frequency-domain representation of the solution (ω) to a Laplace representation (s). Many users start from the frequency-domain representation of solutions of wave propagation problems. In this case issues arising from the movement of singularities under the transformation from ω to s must be considered. We discuss this extension in the context of the Sommerfeld half-plane problem. 相似文献
18.
Mehdi Eshagh 《Pure and Applied Geophysics》2012,169(12):2201-2215
The Earth’s gravity potential can be determined from its second-order partial derivatives using the spherical gradiometric boundary-value problems which have three integral solutions. The problem of merging these solutions by spectral combination is the main subject of this paper. Integral estimators of biased- and unbiased-types are presented for recovering the disturbing gravity potential from gravity gradients. It is shown that only kernels of the biased-type integral estimators are suitable for simultaneous downward continuation and combination of gravity gradients. Numerical results show insignificant practical difference between the biased and unbiased estimators at sea level and the contribution of far-zone gravity gradients remains significant for integration. These contributions depend on the noise level of the gravity gradients at higher levels than sea. In the cases of combining the gravity gradients, contaminated with Gaussian noise, at sea and 250?km levels the errors of the estimated geoid heights are about 10 and 3 times smaller than those obtained by each integral. 相似文献
19.
Elimination of the water-layer response from multi-component source and receiver marine electromagnetic data 总被引：1，自引：0，他引：1
Janniche Iren Nordskag Lasse Amundsen Lars Løseth Egil Holvik 《Geophysical Prospecting》2009,57(5):897-918
This paper presents the theory to eliminate from the recorded multi‐component source, multi‐component receiver marine electromagnetic measurements the effect of the physical source radiation pattern and the scattering response of the water‐layer. The multi‐component sources are assumed to be orthogonally aligned above the receivers at the seabottom. Other than the position of the sources, no source characteristics are required. The integral equation method, which for short is denoted by Lorentz water‐layer elimination, follows from Lorentz' reciprocity theorem. It requires information only of the electromagnetic parameters at the receiver level to decompose the electromagnetic measurements into upgoing and downgoing constituents. Lorentz water‐layer elimination replaces the water layer with a homogeneous half‐space with properties equal to those of the sea‐bed. The source is redatumed to the receiver depth. When the subsurface is arbitrary anisotropic but horizontally layered, the Lorentz water‐layer elimination scheme greatly simplifies and can be implemented as deterministic multi‐component source, multi‐component receiver multidimensional deconvolution of common source gathers. The Lorentz deconvolved data can be further decomposed into scattering responses that would be recorded from idealized transverse electric and transverse magnetic mode sources and receivers. This combined electromagnetic field decomposition on the source and receiver side gives data equivalent to data from a hypothetical survey with the water‐layer absent, with idealized single component transverse electric and transverse magnetic mode sources and idealized single component transverse electric and transverse magnetic mode receivers. When the subsurface is isotropic or transverse isotropic and horizontally layered, the Lorentz deconvolution decouples into pure transverse electric and transverse magnetic mode data processing problems, where a scalar field formulation of the multidimensional Lorentz deconvolution is sufficient. In this case single‐component source data are sufficient to eliminate the water‐layer effect. We demonstrate the Lorentz deconvolution by using numerically modeled data over a simple isotropic layered model illustrating controlled‐source electromagnetic hydrocarbon exploration. In shallow water there is a decrease in controlled‐source electromagnetic sensitivity to thin resistors at depth. The Lorentz deconvolution scheme is designed to overcome this effect by eliminating the water‐layer scattering, including the field's interaction with air. 相似文献
20.
The phenomenon of acoustic waves inducing electric fields in porous media is called the seismoelectric effect. Earlier investigators proposed the usage of seismoelectric effect for well logging. Soil texture has a strong influence on the coupled wave fields during shallow surface explorations. In this article, we study the borehole pure shear‐horizontal wave and the coupling transverse‐electric field (acoustic–electrical coupling wave fields) in the partially saturated soil. Combined with related theories, we expand the formation parameters to partially saturated forms and discuss the influence of soil texture conditions on the seismoelectric wave fields. The results show that the elastic and electrical properties of porous media are sensitive to water saturation. The compositions of the acoustic and electric fields for different soil textures do not change, but the waveforms differ. We also use the secant integral method to simulate the interface‐converted electromagnetic waves. The results show that interface response strength is greatly influenced by soil texture. In addition, considering the sensitivity of the inducing electric field to fluid salinity, we also simulate the time‐domain waveforms of electric field for different pore fluid salinity levels. The results show that as the salinity increases, the electric field amplitude decreases monotonically. The above conclusions have certain significance for the application of borehole shear wave and its coupled electric fields for resource exploration, saturation assessment and groundwater pollution monitoring. 相似文献