共查询到20条相似文献，搜索用时 578 毫秒
1.
Computations of the time-domain electromagnetic response of a multi-layered earth have been carried out for different source-receiver coil systems. The primary excitation is a train of half-sinusoidal waveforms of alternating polarity. The conversion into the time-domain involves Fourier series summation of the matched complex mutual coupling ratios of the layered earth models computed by a digital linear filter method. As an example, the response of a perpendicular coil system on the ground surface for two source-receiver separations has been presented for a five-layer earth model. This has been compared with the responses of homogeneous, two-layer, three-layer, and four-layer models. Next, the investigations have been extended to study the problems of equivalence of three-layer models, the intermediate layer of which is either conductive or resistive. For an intermediate conductive layer (H-type), the studies show that in the early portion of the signal the interpretation of a true three-layer earth is possible to some extent, whereas the ambiguity due to equivalence persists in the late samples. On the other hand, for an intermediate resistive layer (K-type), the three-layer earth and its equivalent model cannot be distinguished from each other over the entire sampling period. On the basis of a computational approach, equivalence has been empirically established as √h/ρ=constant for H-type earth-sections, and as h2ρ=constant for K-type earth sections, where h and ρ are respectively the thickness and resistivity of the intermediate layer. 相似文献
2.
3.
Poddar 《Geophysical Prospecting》1999,47(6):945-958
Long-offset transient electromagnetics (LOTEM) is now regarded as a suitable electrical method for deep exploration along with magnetotellurics (MT). In this method, the vertical magnetic-field impulse response and, occasionally, the horizontal electric-field step response of a grounded-wire source on the surface of the earth are measured. Here, these two responses are computed for 3D models of three deep resistivity structures of interest in hydrocarbon exploration: (i) a faulted graben in a resistive basement rock at a depth of 4 km beneath a conductive overburden; (ii) a facies change in a resistive layer buried at a depth of 2 km in the conductive overburden above a resistive basement; and (iii) an anticlinal uplift of a resistive layer at a depth of 1 km in the conductive overburden above a resistive basement. The results show that the sensitivity of the electric-field response to model perturbation is generally greater than that of the magnetic-voltage response. Further, the electric-field sensitivity is confined to early and intermediate times while that of the magnetic-voltage response is confined to intermediate and late times. Hence it is recommended that both electric and magnetic recordings are made in a LOTEM survey so that the final results can be presented as apparent-resistivity curves derived from the two responses jointly as well as separately. 相似文献
4.
T. LEE 《Geophysical Prospecting》1977,25(1):61-75
The transient response of a layered structure to plane wave excitation can be considered to be composed of a series of waves and a ground wave. For the case of a half-space of conductivity σ and permeability μ the maximum in the electric field is found at a depth z and time t when t=z2σμ/2. This formula can be used to estimate the depth to a buried horizontal conductor with an accuracy that depends upon the resistive contrast at the conductor's surface. The above ray type of solution can be converted to a solution composed of a number of modes by the use of a Poisson transform and the transformed solutions yield decay constants that are consistent with the previously reported results. In the case of a finite source, the maximum in the electric field is strongly directed. The direction depends upon the geometry of the source and the air-earth interface. Although the maximum varies with direction it can be shown that in some directions similar laws to that above are valid. The depth to a conductor can be estimated from the early part of the transients when the ground wave is removed. The removal of the ground wave from the transient is facilitated by the use of an apparent conductivity formula. Although these results were obtained under restrictive conditions they do provide some insight into the electrical transients that are encountered by studying more complex models. 相似文献
5.
Velekin and Bulgakov (1967) in an interesting model experiment while studying the transient electromagnetic response of a conductive sphere placed below a thin conductive sheet found that at the earlier stages of the transience, the composite system response corresponded to the response due to the overlying sheet alone and at the later stages, it corresponded to that of the sphere alone. To examine whether such a separation of responses due to individual components can be analytically studied and applied to other source configurations, we have analyzed an idealized model consisting of two spherical shells. We find that in corroboration with the above results, the general nature of the curve consists of two humps representing the responses dominated by the outer shell and the inner shell respectively. In addition, however, we find that the two humps gradually disappear to yield a smooth decay curve for increasing values of the ratio σd1b/σ2d2a (where σ1, σ2 are the conductivities, d1, d2 are the thicknesses of the outer and inner shells respectively, and b and a are their respective distances from the centre) and the effect of inner shell on the composite system response is considerably reduced. 相似文献
6.
T. LEE 《Geophysical Prospecting》1984,32(3):480-496
A thin superparamagnetic layer on the earth's surface greatly affects the transient electromagnetic response of a conducting ground. The effect of the layer is most evident for singleloop transient electromagnetic data where transient voltages decay as 1/t. Even when a separate transmitter and receiver are used, the effect of the superparamagnetic layer is still pronounced. In this case the effect of the 1/t term in the equation is much less. More dominant now is a 1/t2 term. The effect of the superparamagnetism can readily be seen in the analytical expressions for the apparent resistivities. If the presence of the superparamagnetic layer is not recognized, then the apparent resistivities decrease with time rather than approach the true value of the host rock. 相似文献
7.
A. ERCAN 《Geophysical Prospecting》1981,29(1):89-101
The objective of this study is to generate the separation-distance-domain (r-domain) transformation of the theoretically calculated wave number domain (m-domain) electromagnetic induction field component Bz(m, ω) of a stratified medium and to search for interpretive information which has been absent in the previously achieved numerical solutions of the problem. The r-domain kernel R?(r, ω) function defining the induction field appears to adequately reflect the layering and electrical properties of the medium if it is expressed as a function of the frequency if the source-receiver separation r is small with respect to the thickness of the first layer. However, exact values of the conductivity cannot be distinguished from those of the neighboring values unless a resistive basement layer is present. This feature is the result of the truncation in series representation of the kernel function R?(m, ω). However, this truncation is regarded as significant in the case of a conductive first layer. In m-domain static-zone studies, a conductive first layer slightly influences its r-domain correspondent. Although the computational cost of obtaining the kernel B(r, ω) by evaluation of the convolution in a cylindrical coordinate system is high, this semi-analytic solution is still superior to those based on the asymptotic assumptions. 相似文献
8.
The detectability of an intermediate layer in a three-layer earth model in the time domain has been investigated. The calculations were made for the perpendicular loop (designated system II) and vertical-coplanar (designated system III) electromagnetic (EM) sounding systems. The primary excitation employed is a train of half-sinusoidal and square waveforms of alternating polarity. The time-domain response has been determined by Fourier transformation of the matched complex mutual coupling ratios into the time domain and by linear digital filtering. Top and bottom layers have equal resistivity. EM responses have been computed for conductive and resistive intermediate layer with a wide range of thickness and for two values (500 m and 1000 m) of loop-separation. For the detectability analyses, the root mean square (rms) difference between three-layer and homogeneous-earth responses is adapted. The threshold value for detectability is defined as an rms difference of 10% and the measurement error is arbitrarily assumed to be of the order of 3%. It is observed that the perpendicular-loop system is better than the vertical-coplanar system in detecting thin intermediate layers (either conductive or resistive). For a loop separation of 1000 m and half-sinusoidal pulse excitation, the detectable thickness ratio (h2/h1) is 0.10 by system II for the conducting middle layers; for square pulse excitation the corresponding thickness ratios are 0.06 for system II and 0.12 for system III. For a loop separation of 1000 m and half-sinusoidal pulse excitation in detecting the resistive intermediate layers, the corresponding thickness ratios are 0.9 for system II and 2.25 for system III; while for square pulse excitation the thickness ratios are 0.55 for system II and 1.55 for system III. Results in the frequency domain and time domain (for half-sinusoidal and square pulsed field) have also been presented for systems II and III for detecting conducting layers by considering an earth model where p1≠ p3 and p3 > p1 (p is the resistivity). The loop separa- tion used is 1000 m. Direct comparisons between the frequency domain and time-domain results clearly demonstrate the superiority of frequency-domain systems for detecting con- ducting intermediate layers. 相似文献
9.
D. P. GHOSH 《Geophysical Prospecting》1971,19(2):192-217
Koefoed has given practical procedures of obtaining the layer parameters directly from the apparent resistivity sounding measurements by using the raised kernel function H(λ) as the intermediate step. However, it is felt that the first step of his method—namely the derivation of the H curve from the apparent resistivity curve—is relatively lengthy. In this paper a method is proposed of determining the resistivity transform T(λ), a function directly related to H(λ), from the resistivity field curve. It is shown that the apparent resistivity and the resistivity transform functions are linearily related to each other such that the principle of linear electric filter theory could be applied to obtain the latter from the former. Separate sets of filter coefficients have been worked out for the Schlumberger and the Wenner form of field procedures. The practical process of deriving the T curve simply amounts to running a weighted average of the sampled apparent resistivity field data with the pre-determined coefficients. The whole process could be graphically performed within an quarter of an hour with an accuracy of about 2%. 相似文献
10.
瞬变电磁场的直接时域数值分析 总被引：44，自引：9，他引：35
为了深入了解瞬变电磁场的勘探原理,直接在时间域对负阶跃脉冲激发的二维瞬态场进行了数值分析.采用的方法是从反映电磁场基本规律的麦克斯韦方程组出发,导出时域电场的齐次扩散方程,对所研究的空间区域作差分离散,源作为初始条件加入,利用准静态近似处理空中边界,然后进行时间的逐步递推,由此展现瞬变电磁场在地下扩散随时间发展的全过程.通过模拟计算不同时刻瞬态电场在地下的分布形态及地面上感生电动势相应的变化,揭示了低阻异常体对感应涡流的聚集作用,低阻覆盖层对瞬变场扩散的减速作用,及瞬变场的延时效应.因此,瞬变电磁法对低阻体是敏感的,有上覆低阻层时探测同样的深度需要较长的时间,而延时效应瞬变场的晚期时段可反映埋藏较浅的异常体. 相似文献
11.
Abstract Four geoelectrical soundings were measured with a combination of Schlumberger and azimuthal or equatorial dipole electrode arrays on a Carboniferous limestone basin of the Condroz area, Belgium. The measuring technique is briefly outlined as well as the interpretation procedure, which follows a closed-loop scheme with control of calculated model curves. Some general problems of interpretation of geoelectrical sounding curves are tackled, as far as they have a practical bearing on the treatment of Condroz soundings. The problem of determining the very high resistivity of limestone is approached through ARCHIE's formula, an empirical relation between the bulk rock resistivity, the porosity and the electrolyte resistivity. An evaluation of the latter two parameters, combined with electrical horizontal conductance measurements directly made on resistivity sounding curves, offers a possibility for fast determination of the total water storage in a limestone aquifer. Such storage determinations could be applied whenever an aquifer shows up as a conductive layer interbedded between two highly resistant layers (e.g. nonsaturated limestone and compact, non-fractured limestone). 相似文献
12.
R. HAEFELI 《水文科学杂志》2013,58(2):123-125
Abstract The probability-distributed catchment model, as originally proposed by Moore &; Clarke (1981), is re-examined from a maximum statistical entropy viewpoint. The distribution of water within a catchment is treated as a problem of statistical inference and resolved using an entropy maximization technique. A simple runoff generating mechanism is employed, which, together with the catchment mass balance equation, yields a catchment model involving just one dynamic parameter, y, and two constants, k and λ. The parameter y determines the temporal variation of catchment storage V and runoff q. The latter is nonlinearly related to V through q = k(1—λyV), where y provides the nonlinear departure from the simple linear reservoir q = kV. 相似文献
13.
京津唐渤及其周围地区是我国的强烈地震活动区之一。自1976年以来,我们在该区开展大地电磁测深工作,完成了近30个测点。所得结果表明,本区壳内存在高导层,与地震方法确定的壳内低速层一致。平原内上地幔高导层埋深50-80公里,山区大于100公里,与地震方法确定的上地幔低速层基本一致,同时与大地热流测量、居里等温面计算和对新生代玄武岩地球化学研究结果基本吻合。本区绝大多数地震位于壳内高导层之上,强烈地震主要发生在上地幔高导层隆起的边缘。最后讨论了本区强震活动与壳内和上地幔高导层的关系。 相似文献
14.
R. E. WHITE 《Geophysical Prospecting》1973,21(4):660-703
A seismic trace recorded with suitable gain control can be treated as a stationary time series. Each trace, χj(t), from a set of traces, can be broken down into two stationary components: a signal sequence, αj(t) *s(t—τj), which correlates from trace to trace, and an incoherent noise sequence, nj(t), which does not correlate from trace to trace. The model for a seismic trace used in this paper is thus χj(t) =αj(t) * s(t—τj) +nj(t) where the signal wavelet αj(t), the lag (moveout) of the signal τj, and the noise sequence nj(t) can vary in any manner from trace to trace. Given this model, a method for estimating the power spectra of the signal and incoherent noise components on each trace is presented. The method requires the calculation of the multiple coherence function γj(f) of each trace. γj(f) is the fraction of the power on traced at frequency f that can be predicted in a least-square error sense from all other traces. It is related to the signal-to-noise power ratio ρj(f) by where Kj(f) can be computed and is in general close to 1.0. The theory leading to this relation is given in an Appendix. Particular attention is paid to the statistical distributions of all estimated quantities. The statistical behaviour of cross-spectral and coherence estimates is complicated by the presence of bias as well as random deviations. Straightforward methods for removing this bias and setting up confidence limits, based on the principle of maximum likelihood and the Goodman distribution for the sample multiple coherence, are described. Actual field records differ from the assumed model mainly in having more than one correctable component, components other than the required sequence of reflections being lumped together as correlated noise. When more than one correlatable component is present, the estimate for the signal power spectrum obtained by the multiple coherence method is approximately the sum of the power spectra of the correlatable components. A further practical drawback to estimating spectra from seismic data is the limited number of degrees of freedom available. Usually at least one second of stationary data on each trace is needed to estimate the signal spectrum with an accuracy of about 10%. Examples using synthetic data are presented to illustrate the method. 相似文献
15.
Consider the mathematical model of a horizontally layered system subject to an initial downgoing source pulse in the upper layer and to the condition that no upgoing waveforms enter the layered system from below the deepest interface. The downgoing waveform (as measured from its first arrival) in each layer is necessarily minimum-phase. The net downgoing energy in any layer, defined as the difference of the energy spectrum of the downgoing wave minus the energy spectrum of the upgoing wave, is itself in the form of an energy spectrum, that is, it is non-negative for all frequencies. The z-transform of the autocorrelation function corresponding to the net downgoing energy spectrum is called the net downgoing spectral function for the layer in question. The net downgoing spectral functions of any two layers A and B are related as follows: the product of the net downgoing spectral function of layer A times the overall transmission coefficient from A to B equals the product of the net downgoing spectral function of layer B times the overall transmission coefficient from B to A. The net downgoing spectral function for the upper layer is called simply the spectral function of the system. In the case of a marine seismogram, the autocorrelation function corresponding to the spectral function can be used to recursively generate prediction error operators of successively increasing lengths, and at the same time the reflection coefficients at successively increasing depths. This recursive method is mathematically equivalent to that used in solving the normal equations in the case of Toeplitz forms. The upgoing wave-form in any given layer multiplied by the direct transmission coefficient from that layer to the surface is equal to the convolution of the corresponding prediction error operator with the surface seismogram. The downgoing waveform in this given layer multiplied by the direct transmission coefficient from that layer to the surface is equal to the convolution of the corresponding hindsight error operator (i.e., the time reverse of the prediction error operator) with the surface seismogram. 相似文献
16.
The gravitational attraction produced by a geological body of irregular shape can be easily determined by dividing it into cubes of small size. The exact expression of the effect of a cube is very complicated, but it can be calculated by using an electronic computer. 4851 values of the gravitational attraction were determined for different positions of a cube with the side of 2l and the center in M(x0, y0, z0), for x0∈[0;20], y0∈[0;20] and z0∈[0;10]. Using these values, templates in double logarithmic representation were drawn, with x0 and z0 as parameters and y0 as abscissa, and with x0 and y0 as parameters and z0 in abscissa; this double set of templates permits a good interpolation for all cubes in the considered domain of M. The use of templates was tested to approximate the effect produced by a theoretical model of spherical shape and in a real case of a three-dimensional salt body of known shape based on a large number of boreholes. In both cases very good results were obtained. 相似文献
17.
18.
Boniface C. E. Egboka John A. Cherry Robert N. Farvolden 《Pure and Applied Geophysics》1982,120(2):330-347
The bomb tritium (3H) distribution patterns in the aquifer beneath an abandoned landfill at the Canadian Forces Base (CFB) Borden, Ontario, and in a sandy aquifer at Whiteshell Nuclear Research Establishments (WNRE) Pinawa, Manitoba, all in Canada, were delineated in great detail. A sampling and monitoring network of multilevel samplers and bundle piezometers were used. The directions of groundwater flow were established, and the boundary between the tritiated and non-tritiated zones of the two aquifers were closely demarcated. Using a cumulative mass balance method, the3H input mass into the aquifers was compared with the3H mass in groundwater storage to estimate the percentages of annual groundwater recharge from 1953 to 1978. Two recharge calculations for theeffective recharge zone and thetotal recharge area of the aquifers as established from the flownet analysis, and the distributions of dissolved geochemical constitutents showed that theeffective recharge zone calculations gave higher values of 30.6 cm/yr for CFB Borden and 20.1 cm/yr for WNRE while thetotal recharge areas gave lower values of 19.1 and 10.1 cm/yr for the Borden and WNRE aquifers respectively. The two recharge values provide possible minimum and maximum recharge estimates for the two study areas. 相似文献
19.
E. MUNDRY 《Geophysical Prospecting》1967,15(3):468-479
For the computation of the vertical component Hz of the magnetic field of a horizontal A.C. dipole lying on the earth's surface, a recurrence formula is presented for a horizontally stratified half space, to obtain the (n+ 1)-layer case from the w-layer case. By means of several computed diagrams for the two-layer case, Hz can be determined for different ratios of conductivity of the subsoil and that of the overburden. Thereby the distance from the dipole as well as the layer thickness h are expressed in terms of the wave length A of a plain wave in the overburden. Assuming a sufficiently large conductivity difference, the results show that evidence about the subsurface conditions can be obtained if the distance between the measuring coil and the dipole is of the order of A/3, and if the thickness h of the layer varies within the range A/100 < h < A/6. As an example for the 3-layer case, a nonconducting intermediate layer is assumed. 相似文献
20.
An infinitely resistive/conductive horizontal bed is assumed in an otherwise homogeneous and isotropic half space. Schlumberger, three electrode, and unipole profiles are computed at right angles to the strike of the bed. The Schwarz-Christoffel method of conformal transformation and numerical methods of solving non-linear differential equations are used to solve the boundary value problem. It is observed that (i) the three electrode system is the most sensitive gradient electrode configurations for electrical profiling, (ii) the apparent resistivities for Schlumberger, three electrode, and unipole methods become maximum when the depth of the bed is 0.06 L, 0.1 L, and 0.055 L for a resistive bed and minimum when depths are 0.085 L, 0.04 L-0.02 L and indeterminate for conductive beds, respectively, (iii) the limiting depths of detection (defined in the text) by Schlumberger, three electrode, and unipole configurations are respectively 0.9 L, 6.6 L and 2.0 L for resistive beds and 0.58 L, 1.17 L and 1.5 L for conductive beds. The electrode separation L is the distance between the two farthest active electrodes. 相似文献