共查询到20条相似文献，搜索用时 22 毫秒
1.
A new formulation is proposed for the electrical potential developed inside a horizontally‐layered half‐space for a direct current point‐source at the surface. The recursion formula for the kernel coefficient in the potential integral is simpler than the generally used two‐coefficient recursion. The numerical difficulties that may occur during the computation of the integrals and near the source axis are examined and solutions are proposed. The set of equations permits a stable and accurate computation of the tabular potential everywhere in the medium. 相似文献
2.
Nagendra Pratap Singh Mitsuru Utsugi Tsuneomi Kagiyama 《Pure and Applied Geophysics》2009,166(12):2037-2058
In the present research, an attempt is made to derive the generalized expressions for the transient electromagnetic response of a large loop source over the surface of a homogeneous earth model for arbitrary receiver positions inside or outside the source loop. Expressions are derived for the impulse as well as step excitations of the source loop. As a cross check for validity of expressions, the step response expression is obtained from the impulse response expression and vice versa. Computations are performed for the TEM response over a homogeneous earth model for source-receiver offset (r = 0) pertinent to the central loop configuration and the results are compared with the published results for the central loop TEM responses. The results are in well coincidence with each other and thus provide the check for the authenticity of the expressions. To exemplify the nature of TEM response at various source receiver offsets, results are presented for the TEM response at source receiver offsets r = 0, r = a/2, r = a and r = 2a, 6a, 12a pertaining to the central loop, in-loop, on-loop and offset loop configurations, respectively. The results depict their characteristic variations. At receiver positions inside the loop source, both the impulse as well as step responses are of same sign, whereas at the receiver positions outside the loop source, both the curves exhibit a change of sign that shifts towards the later times with increase in the offset distances. The change of sign in impulse response occurs at a relatively later time than that in the step response. This is the initial presentation of TEM response expressions for the large loop source over a homogeneous earth model for arbitrary receiver position inside and/ or outside the loop source except for the case of receiver at the center of the loop and at the coincident loop point. This research would be of immense use in the development and use of the large loop TEM method in its various configurations and thus would enhance the applicability and cost effectiveness of the large loop source TEM method. 相似文献
3.
本文给出一种用积分方程法配合有限差分法计算位于导电介质中三维导体的电磁响应的新途径,在二层大地条件下给出了具体计算方法和结果。详细讨论了围岩介质的导电性对导体异常的影响。这对于开展低阻覆盖层地区的电磁法工作具有一定的指导意义。 相似文献
4.
5.
M. PODDAR 《Geophysical Prospecting》1983,31(3):508-523
An electromagnetic sounding experiment with a large square loop as source was carried out on the dried-up bed of a water reservoir near the town of Dharmavaram in Andhra Pradesh. The sounding was performed in both geometric and parametric modes, and involved measurements of phase as well as amplitude of the vertical magnetic field inside and outside the loop. The six-frequency EM system used for the experiment was found to be a workable system for electromagnetic sounding. The results of the experiment more or less confirm the conclusions from theoretical modeling. They also show that even though the earth is not always electrically horizontally layered over a sufficiently large lateral extent for the one-dimensional model to be strictly valid, it is still possible to apply such modeling to sounding curves taken one part at a time and obtain layer parameters which check qualitatively with the layer parameters obtained from direct current resistivity sounding. 相似文献
6.
In present‐day land and marine controlled‐source electromagnetic (CSEM) surveys, electromagnetic fields are commonly generated using wires that are hundreds of metres long. Nevertheless, simulations of CSEM data often approximate these sources as point dipoles. Although this is justified for sufficiently large source‐receiver distances, many real surveys include frequencies and distances at which the dipole approximation is inaccurate. For 1D layered media, electromagnetic (EM) fields for point dipole sources can be computed using well‐known quasi‐analytical solutions and fields for sources of finite length can be synthesized by superposing point dipole fields. However, the calculation of numerous point dipole fields is computationally expensive, requiring a large number of numerical integral evaluations. We combine a more efficient representation of finite‐length sources in terms of components related to the wire and its end points with very general expressions for EM fields in 1D layered media. We thus obtain a formulation that requires fewer numerical integrations than the superposition of dipole fields, permits source and receiver placement at any depth within the layer stack and can also easily be integrated into 3D modelling algorithms. Complex source geometries, such as wires bent due to surface obstructions, can be simulated by segmenting the wire and computing the responses for each segment separately. We first describe our finite‐length wire expressions and then present 1D and 3D examples of EM fields due to finite‐length sources for typical land and marine survey geometries and discuss differences to point dipole fields. 相似文献
7.
Summary The basic theory of airborne EM surveying, in the time domain, is considered. Rather than resorting to tediousdouble numerical integration, a more direct approach is adopted. This method, valid in the quasi-static regime, is illustrated for a homogeneous flat earth. The results exhibit a number of clear-cut features that are relevant to remote sensing. For example, a vertical co-axial loop system has a desirable transient response from the standpoint of yielding conductivity data without requiring accurate information on the height of the device above the ground.The research reported here was supported by AFCRL Contract No. PRO-Y-71-872. 相似文献
8.
S.H. Hall 《Geophysical Prospecting》1996,44(2):251-277
The rotating current EM method has been applied to the delineation of two conductive orebodies, Elura near Cobar, NSW, and Thalanga near Charter's Towers, Queensland. The field data were collected in the form of observations of the vertical magnetic field strength ratio and phase difference using a Turam-style receiver with twin vertical coils. By reconstituting this data back to the ring source field and phase, i.e. the observed Hz, phasor, it is possible to present contoured maps of the EM field. Anomaly phasors are obtained by subtracting theoretical phasors from the observed phasors in the complex plane of the Hz phasor. The theoretical phasors for the finite source are based on horizontally layered, half-space earth models, computed at each point of the survey grids, then normalized to selected points of the observed fields. Use is made of the intrinsic circular symmetry of the method in X–Y plots of field versus source-receiver distance to ascertain geoelectric parameters for the earth models. A steel picket fence at Thalanga is modelled by a line source grounded at each end and its Hz, phasor is removed by the same process. A considerable improvement in anomaly delineation is gained over previous Turam-style anomalies and the two survey examples illustrate the limitations of the method in the presence of a conductive overburden (Elura) and its abilities in the absence of a conductive overburden (Thalanga). 相似文献
9.
The electric and magnetic fields generated by an individual horizontal current ring induced inside a homogeneous conductive half‐space, originating from an external large circular loop source of current in the presence of a flat half‐space, are deduced. A check of self‐consistency for these expressions led to the known general functions for these fields due to the same external source in the presence of that medium. The current rings’ mutual coupling related to the magnetic field's radial component is thoroughly analysed and its specific members are presented. The existence of a relatively small zone inside the half‐space responsible for the main contribution for the signal measured at the observation point, with the source and receiver on the ground surface, is made evident. For increasing values of frequency, at a given transmitter‐receiver (T–R) configuration, this zone shrinks and its central point moves away from a maximum depth of about 30% and horizontal distance of nearly 85%, of the T–R separation, to a point very close to the receiver position. The coordinates of the central point of this zone of main contribution are provided as approximated functions in terms of the induction number . 相似文献
10.
The elastic Kirchhoff-Helmholtz integral expresses the components of the monochromatic displacement vector at any point A in terms of the displacement field and the stress field at any closed surface surrounding A. By introducing Green's functions for P- and S-waves, the elastic Kirchhoff-Helmholtz integral is modified such that it expresses either the P-wave or the S-wave at A in terms of the elastic wavefield at the closed surface. This modified elastic Kirchhoff-Helmholtz integral is transformed into one-way elastic Rayleigh-type integrals for forward extrapolation of downgoing and upgoing P- and S-waves. We also derive one-way elastic Rayleigh-type integrals for inverse extrapolation of downgoing and upgoing P- and S-waves. The one-way elastic extrapolation operators derived in this paper are the basis for a new prestack migration scheme for elastic data. 相似文献
11.
M. PODDAR 《Geophysical Prospecting》1982,30(1):101-114
Using the concept of reciprocity and the known solution for the electric field of a vertical oscillating magnetic dipole source placed over a two-layered halfspace, an integral expression for the vertical magnetic field produced by a horizontal rectangular loop, carrying an oscillating current and placed on the surface of the same halfspace, is deduced. This expression is such that it can be evaluated by a combination of straightforward numerical integration and digital linear filter techniques. Displacement currents everywhere in space are neglected. Fields both inside and outside the loop are presented as sounding curves. A little over a decade of frequencies (100–3000 Hz) appears to be the minimum requirement and two decades (100–10000 Hz) are sufficient for most situations. Parametric and geometric sounding curves both show good resolution of subsurface layering. Phase shows better resolution than amplitude. Similarly, resolution is higher when conductivity increases with depth than when it decreases. Besides being useful for the electromagnetic depth sounding for layered earth structures, such computations can be applied to determine normal corrections to Turam observations, whenever a rectangular loop of finite size is used as a transmitter. 相似文献
12.
A. T. BAOKUR 《Geophysical Prospecting》1984,32(1):132-138
Different sets of filter coefficients for the linear filter technique for the computations of resistivity and EM sounding curves are evaluated for several electrode and coil configurations. Instead of this procedure, the two-electrode filter can be used for computations of Wenner, Schlumberger, and dipole—dipole apparent resistivity model curves by defining convolutional expressions which contain the new input functions in terms of the resistivity transform function. Similarly, the Schlumberger filter performs the computations of dipole—dipole apparent resistivity model curves. The Wenner, Schlumberger, and dipole—dipole filter functions are defined in terms of the two-electrode filter using the new convolutional expressions. A relationship between the Schlumberger and dipole—dipole filter functions is given. The above arguments are adopted for the computations of EM sounding curves. It is shown that the EM filter for the horizontal coplanar loop system (which is identical to the two-electrode filter) performs the computations of the mutual coupling ratios for perpendicular, vertical coplanar, and vertical coaxial loop systems. In the same way, the Schlumberger filter can be used to compute vertical coaxial sounding curves. The corresponding input functions are defined in terms of the EM kernel for all convolutional expressions presented. After these considerations, integral expressions of the mutual coupling ratios involving zero-order Bessel function are derived. The mutual coupling ratio for the vertical coaxial loop system is given in the same form as the mutual coupling ratio for the vertical coplanar loop system. 相似文献
13.
14.
The two-dimensional problem of a long displacement dislocation in a multilayered half-space is studied. Both plane strain and antiplane strain cases are considered. The Thomson-Haskell matrix method is used to obtain the static field caused by the dislocation. The source is represented in terms of the jumps at the source level in the depth-dependent factors in the Fourier integral representation of the displacements and stresses due to the source in an infinite medium. Explicit expressions for the surface displacements due to dip-slip and strike-slip faults of arbitrary dip are obtained. The formulation developed is quite convenient for numerical computation. 相似文献
15.
In this paper, closed-form analytic expressions for the frequency-wave number domain Fourier amplitudes of the displacement field at the free surface of a layered, anelastic half-space are established. The displacement field is caused by a seismic source described by a shear dislocation propagating with constant velocity over a rectangular fault (Haskell's model). Three-dimensional plane wave propagation is considered in the layered half-space using a propagator-based formalism. The wave radiation from the source is decoupled into P-SV and SH motions and the two problems are treated separately. First, analytic expressions are calculated for the displacement field at the free surface due to unidirectional unit impulses. Then, these expressions are used to compute solutions for the displacement field due to effective point sources associated with a pure strike slip and a pure dip slip. Finally, these solutions are combined and integrated over the rectangular fault area to establish closed-form analytic expressions of the total displacement field at the free surface. 相似文献
16.
D. PATELLA 《Geophysical Prospecting》1972,20(3):561-579
In this paper it is shown how one may obtain a generalized Ohm's law which relates the induced polarization electric field to the steady-state current density through the introduction of a fictitious resistivity defined as the product of the chargeability and the resistivity of a given medium. The potential generated by the induced polarization is calculated at any point in a layered earth by the same procedure as used for calculating the potential due to a point source of direct current. On the basis of the definition of the apparent chargeability ma, the expressions of ma for different stratigraphie situations are obtained, provided the IP measurements are carried out on surface with an appropriate AMNB array. These expressions may be used to plot master curves for IP vertical soundings. Finally some field experiments over sedimentary formations and the quantitative interpretation procedure are reported. 相似文献
17.
本文研究声电效应测井波场的有限差分模拟算法.忽略井外地层中诱导电磁场对孔隙弹性波的影响,将求解动电耦合波方程组的问题解耦,先计算孔隙弹性波,再计算其诱导电磁场.基于轴对称柱坐标系下的速度-应力交错网格,采用时域有限差分计算井孔流体声波和井外地层孔隙弹性波.将电磁场近似看作似稳场,基于轴对称柱坐标系下的5点式有限差分网格,求解不同时刻的电位Poisson方程,计算诱导电场.结果表明:本文算法可准确模拟频率6.0 kHz的声电效应测井全波;在声波测井频率范围内,电导率、动电耦合系数和动态渗透率的低频近似对伴随电磁场的计算影响不大;地层水平界面导致伴随反射斯通利波的电场和显著的界面电磁波,后者对于探测地层界面具有潜在的应用价值. 相似文献
18.
19.
Complex image method for calculating electric and magnetic fields produced by an auroral electrojet of finite length 总被引：3，自引：0，他引：3
The electromagnetic field due to ionospheric currents has to be known when evaluating space weather effects at the earths surface. Forecasting methods of these effects, which include geomagnetically induced currents in technological systems, are being developed. Such applications are time-critical, so the calculation techniques of the electromagnetic field have to be fast but still accurate. The contribution of secondary sources induced within the earth leads to complicated integral formulas for the field at the earths surface with a time-consuming computation. An approximate method of calculation based on replacing the earth contribution by an image source having mathematically a complex location results in closed-form expressions and in a much faster computation. In this paper we extend the complex image method (CIM) to the case of a more realistic electrojet system consisting of a horizontal line current filament with vertical currents at its ends above a layered earth. To be able to utilize previous CIM results, we prove that the current system can be replaced by a purely horizontal current distribution which is equivalent regarding the total (= primary + induced) magnetic field and the total horizontal electric field at the earths surface. The latter result is new. Numerical calculations demonstrate that CIM is very accurate and several magnitudes faster than the exact conventional approach. 相似文献
20.
An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed. 相似文献