首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
Airborne EM skin depths   总被引:1,自引:0,他引:1  
Skin depth is an electromagnetic (EM) scale length that provides a measure of the degree of attenuation experienced by a particular frequency of an EM system. As has been discussed in the literature, skin depth is not a complete measure of the depth of investigation, but the two may be related. Frequency‐domain airborne EM systems employ pairs of transmitter and receiver coils that use a frequency range from several hundred hertz to over 100 kHz. For elevated dipoles, both geometrical and frequency‐dependent attenuation of the induced fields must be considered. For airborne EM systems it is possible to define a skin depth based only on the electric field induced by the transmitter. A vertical decay scale length, here defined from the at‐surface position of maximum electric field, enables the same skin‐depth estimate to be obtained for both cases of vertical and horizontal dipolar excitation. Such dipolar skin depths associated with towed‐bird and fixed‐wing airborne systems are studied in relation to frequency, conductivity and sensor elevation. Dipolar skin depths are found to be much smaller than their plane‐wave counterparts except at high frequency (>50 kHz) and in combination with high conductivity. For the majority of airborne systems the influence of altitude on skin depth is highly significant. Dipolar skin depths increase with increasing sensor elevation. Low frequencies display the greatest sensitivity. At low elevation (<40 m), geometrical attenuation dominates the behaviour of the skin depth. The study indicates that typical low‐altitude airborne surveys provide vertically compact assessments of subsurface conductivity, well suited to near‐surface, environmental applications.  相似文献   

2.
Offsite movement of waterborne agrochemicals is increasingly targeted as a non-point source of water quality degradation. Our research has indicated that subsurface water movement is variable and site-specific, and that a small soil volume frequently conducts a large volume of flow. This concentrated flow is usually caused by soil morphology, and it often results in water moving rapidly offsite from certain areas of fields; little or no lateral subsurface flow may occur in other areas. Identifying these subsurface regions is difficult using conventional soil survey and vadose zone sampling techniques. In this study, traditional surveying is combined with electromagnetic induction (EMI) and ground-penetrating radar (GPR) mapping to identify areas with high potential for subsurface offsite movement of agrochemicals, optimizing these identification techniques, and expanding the mapping procedures to make them useful at the field-scale for agricultural production practices. Conclusions from this research are: (1) EMI mapping provides rapid identification of areas of soil with a high electrical conductivity and presumably high potential for offsite movement of subsurface water, (2) GPR mapping of areas identified by EMI mapping provides a means to identify features that are known to conduct concentrated lateral flow of water, and (3) combining the capabilities of EMI and GPR instrumentation makes possible the surveys of large areas that would otherwise be impossible or unfeasible to characterize.  相似文献   

3.
This paper presents a computational method for the interpretation of electromagnetic (EM) profile data in the frequency domain using a thin plate model within a two-layer earth. The modelling method is based on an integral equation formulation, where the conductor is represented by a lattice structure composed of two-dimensional surface elements. Several approximations are used to simplify the theoretical basis and to decrease the computation time. The simple parametric model allows efficient use of optimization methods. We employ a linearized inversion scheme based on singular value decomposition and adaptive damping. The new forward computation method and the parameter optimization are combined in the computer program, emplates . The modelling examples demonstrate that the approximate method is capable of describing the characteristic behaviour of the EM response of a thin plate-like conductor in conductive surroundings. The efficacy of the inversion is demonstrated using both synthetic and field data. An optional depth compensation method is used to improve the interpreted values of the depth of burial. The results show that the method is cost effective and suitable for interactive interpretation of EM data.  相似文献   

4.
双平方根单程波动方程叠前τ偏移方法   总被引:1,自引:0,他引:1       下载免费PDF全文
本文将常规双平方根(DSR)单程波动方程从深度域变换到双程垂直走时(τ)域,由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正,可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件,提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难,本文采用稳相近似,在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明,DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.  相似文献   

5.
Pure-mode wave propagation is important for applications ranging from imaging to avoiding parameter tradeoff in waveform inversion. Although seismic anisotropy is an elastic phenomenon, pseudo-acoustic approximations are routinely used to avoid the high computational cost and difficulty in decoupling wave modes to obtain interpretable seismic images. However, such approximations may result in inaccuracies in characterizing anisotropic wave propagation. We propose new pure-mode equations for P- and S-waves resulting in an artefact-free solution in transversely isotropic medium with a vertical symmetry axis. Our approximations are more accurate than other known approximations as they are not based on weak anisotropy assumptions. Therefore, the S-wave approximation can reproduce the group velocity triplications in strongly anisotropic media. The proposed approximations can be used for accurate modelling and imaging of pure P- and S-waves in transversely isotropic media.  相似文献   

6.
New methods for obtaining and quantifying spatially distributed subsurface moisture are a high research priority in process hydrology. We use simple linear regression analyses to compare terrain electrical conductivity measurements (EC) derived from multiple electromagnetic induction (EMI) frequencies to a distributed grid of water‐table depth and soil‐moisture measurements in a highly instrumented 50 by 50 m hillslope in Putnam County, New York. Two null hypotheses were tested: H0(1), there is no relationship between water table depth and EC; H0(2), there is no relationship between soil moisture levels and EC. We reject both these hypotheses. Regression analysis indicates that EC measurements from the low frequency EM31 meter with a vertical dipole orientation could explain over 80% of the variation in water‐table depth across the test hillslope. Despite zeroing and sensitivity problems encountered with the high frequency EM38, EC measurements could explain over 70% of the gravimetrically determined soil‐moisture variance. The use of simple moisture retrieval algorithms, which combined EC measurements from the EM31 and EM38 meters in both their vertical and horizontal orientations, helped increase the r2 coefficients slightly. This first hillslope hydrological analysis of EMI technology in this way suggests that it may be a promising method for the collection of a large number of distributed soilwater and groundwater depth measurements with a reasonable degree of accuracy. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

7.
This study presents the potential of multi-frequency electromagnetic induction (EMI) in archaeology. EMI is currently less employed for archaeological prospection with respect to other geophysical techniques. It is capable of identifying shallow subsurface relics by simultaneously measuring the apparent electrical conductivity (ECa) and apparent magnetic susceptibility (MSa). Moreover, frequency sounding is able to quantify the depths and vertical shapes of buried structures. In this study, EMI surveys with five frequencies were performed at two heritage sites with different geological conditions: Han Hangu Pass characterized by cinnamon soil and Xishan Yang by sandy loams. In the first site, high ECa values were observed with variations in depth correlated to archaeological remains. Moreover, electromagnetic anomalies related to an ancient road and five kiln caves were identified. In the second site, an ancient tomb, indicating extremely low ECa and high MSa, was discovered. Its electromagnetic properties are attributed to the cavity and ferroferric oxides.  相似文献   

8.
Reflection full waveform inversion can update subsurface velocity structure of the deeper part, but tends to get stuck in the local minima associated with the waveform misfit function. These local minima cause cycle skipping if the initial background velocity model is far from the true model. Since conventional reflection full waveform inversion using two‐way wave equation in time domain is computationally expensive and consumes a large amount of memory, we implement a correlation‐based reflection waveform inversion using one‐way wave equations to retrieve the background velocity. In this method, one‐way wave equations are used for the seismic wave forward modelling, migration/de‐migration and the gradient computation of objective function in frequency domain. Compared with the method using two‐way wave equation, the proposed method benefits from the lower computational cost of one‐way wave equations without significant accuracy reduction in the cases without steep dips. It also largely reduces the memory requirement by an order of magnitude than implementation using two‐way wave equation both for two‐ and three‐dimensional situations. Through numerical analysis, we also find that one‐way wave equations can better construct the low wavenumber reflection wavepath without producing high‐amplitude short‐wavelength components near the image points in the reflection full waveform inversion gradient. Synthetic test and real data application show that the proposed method efficiently updates the background velocity model.  相似文献   

9.
Alight nonaqueous phase liquid (LNAPL) ground water contaminant plume has been discovered by purely geophysical means at the former Wurtsmith Air Force Base (AFB) near Oscoda, Michigan. It is located near another plume called FT-02, which is a well-studied area undergoing natural bioremediation. The plume was discovered by ground penetrating radar (GPR) profiling while extending a long line from FT-02 to establish background variability around that plume. The new plume was apparent because of a high-conductivity "shadow' or GPR reflection attenuation observed below the conductive zone at the top of the aquifer, identical to the pattern observed at the FT-02 plume. Further GPR surveys were conducted by students of a Western Michigan University geophysics field course to outline the proximal part of the plume. The GPR survey was supplemented by an electromagnetic induction (EM) survey which showed a group of four cables crossing the area. Finally, a magnetometer survey was conducted to search for any buried steel objects which might have been missed by the EM survey. The results of the three geophysical surveys were then used by students of a University of Michigan field course to guide subsurface soil and fluid sampling, which verified the presence of residual LNAPL product and ground water with conductivities 2.5 to 3.3 times above background. The plume source is in the vicinity of a vaulted underground storage tank (UST) formerly used for the collection of waste solvents and fuels for subsequent use in the fire training exercises at FT-02. This newly discovered LNAPL plume, along with other "mature' plumes, fits the electrical model which predicts conductive ground water below the decomposing but electrically resistive LNAPLs. Finally, this is a fine example of the cooperative use of a dedicated research site for training by students of two different universities.  相似文献   

10.
An algorithm for modelling and migrating ground penetrating radar (GPR) data in moderately heterogeneous dispersive media is presented. The method is based on wavefield extrapolation in the frequency–wavenumber (fk) domain, from the solution of the 2D Maxwell's equations. The wavefield is extrapolated by a phase-shift technique using a constant relative permittivity K and a quality factor Q. It is then modified by a correction term to handle the lateral K and Q variations. The spatial distribution of the K and Q-factor values, representing the given model parameters, is introduced into the algorithm by a regular grid parametrization. The radar wave dispersion and attenuation, induced by relaxation processes, are taken into account by a linear frequency-dependent Q model, and expressed by a complex wavenumber in the propagation equation. A synthetic case and a field data set illustrate the potential of the method for frequencies of 300, 500 and 900 MHz. In the first case, a typical civil engineering problem is considered. The frequency dependence of the wave velocity and attenuation is well illustrated. The synthetic data are afterwards migrated using the initial model parameters. The results show the importance of using spatially varying model parameters in the migration processes. The second case concerns an application of the method to a real data set. In order to adjust the model parameters, a forward modelling sequence is performed until the best match between the measured and the synthetic data is achieved. A depth migration is then applied to the data, and the result is compared with the initial model parameters. In conclusion, we assess the contributions of the method to industrial applications, by discussing the performance of the algorithm compared with its limitations.  相似文献   

11.
Abstract

The stratification profile of the Earth's magnetofluid outer core is unknown, but there have been suggestions that its upper part may be stably stratified. Braginsky (1984) suggested that the magnetic analog of Rossby (planetary) waves in this stable layer (the ‘H’ layer) may be responsible for a portion of the short-period secular variation. In this study, we adopt a thin shell model to examine the dynamics of the H layer. The stable stratification justifies the thin-layer approximations, which greatly simplify the analysis. The governing equations are then the Laplace's tidal equations modified by the Lorentz force terms, and the magnetic induction equation. We linearize the Lorentz force in the Laplace's tidal equations and the advection term in the magnetic induction equation, assuming a zeroth order dipole field as representative of the magnetic field near the insulating core-mantle boundary. An analytical β-plane solution shows that a magnetic field can release the equatorial trapping that non-magnetic Rossby waves exhibit. A numerical solution to the full spherical equations confirms that a sufficiently strong magnetic field can break the equatorial waveguide. Both solutions are highly dissipative, which is a consequence of our necessary neglect of the induction term in comparison with the advection and diffusion terms in the magnetic induction equation in the thin-layer limit. However, were one to relax the thin-layer approximations and allow a radial dependence of the solutions, one would find magnetic Rossby waves less damped (through the inclusion of the induction term). For the magnetic field strength appropriate for the H layer, the real parts of the eigenfrequencies do not change appreciably from their non-magnetic values. We estimate a phase velocity of the lowest modes that is rather rapid compared with the core fluid speed typically presumed from the secular variation.  相似文献   

12.
由所建立的三维qP波相速度表示式出发,导出并解析求解各向异性介质中的频散方程,得到三维各向异性介质中的相移算子,进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质,发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化;文中给出的算例虽然是针对二维VTI介质的,但所提出的算法同样适用于三维TI介质.  相似文献   

13.
Complementary geophysical surveys on large landslides help revealgeologic structures and processes, and thus can help devise mitigation strategies. The combined interpretation of these methods enhance the result of each data set interpretation and makes it possible to derive a geological model of the landslide.We chose a test site on the Boup landslide (Wallis, western SwissAlps) to test high-resolution seismic reflection surveyscombined with ground penetratingradar (GPR), electromagnetic (EM) and electrokinetic spontaneous potential (SP) measurements.The results of the high-resolution seismic surveys suggest thatthe sliding is within a gypsum layer at 50 m depth and not as previously believed along a deeper (70 m) gypsum-shale boundary, also mapped seismically. Inversion of electromagnetic profiles (EM-34) with constraints from seismic data provided a model cross-section of conductivity values of the landslide (20–25 mS/m) and of the surrounding stable ground (10–15 mS/m), and it helped outline their boundary at depth.The accurate surface location of the landslide limit could be detected withelectromagnetic measurements of shallower depths of investigation (EM-31). Positive PS anomalies revealed an upward flow of mineralised water interpreted to follow the lateral boundary of the Boup landslide on its east side. Limited success was obtained withGPR profiling. This method can be hampered by conductive shallow layers, and itssuccessful application on landslides is expected to be strongly site dependant.  相似文献   

14.
Migration velocity analysis by depth image-wave remigration: first results   总被引:2,自引:0,他引:2  
The image‐wave equation for depth remigration is a partial differential equation that is similar to the acoustic wave equation. In this work, we study its finite‐difference solution and possible applications. The conditions for stability, dispersion and dissipation exhibit a strong wavenumber dependence. Where higher horizontal than vertical wavenumbers are present in the data to be remigrated, stability may be difficult to achieve. Grid dispersion and dissipation can only be reduced to acceptable levels by the choice of very small grid intervals. Numerical tests demonstrate that, upon reaching the true medium velocity, remigrated images of curved reflectors propagate to the correct depth and those of diffractions collapse to single points. The latter property points towards the method's potential for use as a tool for migration velocity analysis. A first application to inhomogeneous media shows that in a horizontally layered medium, the reflector images reach their true depth when the remigration velocity equals the inverse of the mean medium slowness.  相似文献   

15.
In this work we analyse the applicability of amplitude grid maps to the routine of geological surveys by means of GPR (ground penetrating radar). Although amplitude grid maps have been commonly used in archaeological surveys, their use in geological prospecting (including the detection of voids and determination of the internal geometry of sedimentary bodies) is not widespread. The direct analysis of GPR-profiles permits the analysis of geometrical features and other qualitative aspects that can be related to changes in EM properties. Aspects such as changes in the density of the banded disposition in radargrams, loss of reflector definition or higher scattering in particular zones of the profiles can give useful, though non-quantitative, information. The GPR wave-amplitude is a qualitative measurement of magnetic properties that can be processed as a semi-quantitative layout. The main differences observed in changing wave amplitude are related to the surveyed materials and their geometry. These changes produce variations in the relative wave amplitude or vertical wave-phase changes related to differences in the propagation velocity, the attenuation factor, the reflectivity and the geometry of the materials. Maps based on the lateral correlation between profiles (C-scans) or the lateral correlation of wave amplitudes along the same profile (as a tomography or Am-scans) permit the analysis of these changes. Variations in amplitude grid maps or Am-scans are related to (i) geometrical changes of surveyed materials, (ii) changes in the dielectric constant, and (iii) changes in the potential penetration depth (higher attenuation in particular zones of a profile) than can be correlated to the type, state or clay content of subsoil materials. Direct analysis of exposures helps to constrain interpretation using the geometrical features in radargrams resulting from the geological structure. In the same way, analysis of geometrical features in radargrams, together with time-slices of C-scans, can be used to determine the areal distribution of changes in the subsoil and approach the changes in EM properties. An example with parallel profiles and different devices, constrained by means a broadband multifrequency EM survey is shown.  相似文献   

16.
三维VTI介质中波动方程深度偏移的最优分裂Fourier方法   总被引:17,自引:5,他引:12       下载免费PDF全文
从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点.  相似文献   

17.
Optimal electromagnetic wave propagation velocities and subsurface images for ground-penetrating radar (GPR) data can be specified by using an imaging scanning method. In addition to time-migrating the unmigrated GPR section, we remigrate the already time-migrated section by a one-step remigration operator using different velocities. This creates many time-migrated images for different constant migration velocities. In this way, the computation time for time-migration is very much reduced. Time-migrated reflector images `propagate' when the constant migration velocity is continuously changed. For this `propagation' there exists a wave-equation-type partial differential equation. Each time-migrated section can thus be viewed as a snapshot for a certain migration velocity. The time-migrated reflector images behave like `waves', called image waves. This is applied to real GPR data acquired over a concrete body within which a steel cable frame is buried. The method produces a quick velocity scan to find a reliable migration velocity leading to the best time-migrated image.  相似文献   

18.
Electrical, seismic, and electromagnetic methods can be used for noninvasive determination of subsurface physical and chemical properties. In particular, we consider the evaluation of water salinity and the detection of surface contaminants. Most of the relevant properties are represented by electric conductivity, P-wave velocity, and dielectric permittivity. Hence, it is important to obtain relationships between these measurable physical quantities and soil composition, saturation, and frequency. Conductivity in the geoelectric frequency range is obtained with Pride's model for a porous rock. (The model considers salinity and permeability.) White's model of patchy saturation is used to calculate the P-wave velocity and attenuation. Four cases are considered: light nonaqueous phase liquid (LNAPL) pockets in water, dense nonaqueous phase liquid (DNAPL) pockets in water, LNAPL pockets in air, and DNAPL pockets in air. The size of the pockets (or pools), with respect to the signal wavelength, is modeled by the theory. The electromagnetic properties in the GPR frequency range are obtained by using the Hanai–Bruggeman equation for two solids (sand and clay grains) and two fluids (LNAPL or DNAPL in water or air). The Hanai–Bruggeman exponent (1/3 for spherical particles) is used as a fitting parameter and evaluated for a sand/clay mixture saturated with water.Pride's model predicts increasing conductivity for increasing salinity and decreasing permeability. The best-fit exponent of the Hanai–Bruggeman equation for a sand/clay mixture saturated with water is 0.61, indicating that the shape of the grains has a significant influence on the electromagnetic properties. At radar frequencies, it is possible to distinguish between a water-saturated medium and a NAPL-saturated medium, but LNAPL- and DNAPL-saturated media have very similar electromagnetic properties. The type of contaminant can be better distinguished from the acoustic properties. P-wave velocity increases with frequency, and has dissimilar behaviour for wet and dry soils.  相似文献   

19.
接地导线源电磁场全域有效趋肤深度   总被引:3,自引:2,他引:1       下载免费PDF全文
为了更好地推进广域电磁法的发展和应用,本文以接地长导线源为例,研究了可控源电磁场全场域的有效趋肤深度.利用频域电偶极源在均匀半空间产生电磁场的闭合表达式,计算了不同电磁场分量定义的有效趋肤深度,并讨论了在不同频率、不同偏移距、不同电导率情况下,有效趋肤深度的变化特性.根据有效趋肤深度随偏移距的变化特征以及与平面波趋肤深度之间的关系,利用多项式拟合的办法在五个不同频率范围内给出了适用于全场域的有效趋肤深度快速估算公式.研究结果表明:不同分量定义的有效趋肤深度是不同的,但是它们随偏移距、频率、电导率等参数的变化趋势是类似的,而且在一定范围内都趋近于平面波趋肤深度.理论模型的研究表明,有效趋肤深度可以作为测量参数选取和数据解释工作的参考依据.  相似文献   

20.
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号