共查询到20条相似文献，搜索用时 83 毫秒
1.
The travel time and amplitude of ground-penetrating radar (GPR) waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data. 相似文献
2.
针对Cole-Cole频散介质中的复介电常数是jω的分数次幂函数,传统的时域有限元法难以离散及计算时间域分数阶导数,本文采用Pade逼近算法将含有时间分数阶导数的Cole-Cole频散介质电磁波方程推导为一组整数阶辅助微分方程,提出了一种适用于Cole-Cole频散介质的GPR有限元正演模拟算法.在复数伸展坐标系下,通过在频率域Cole-Cole频散介质电磁波方程中引入2个中间变量,并将其变换到时间域,从而以变分形式将PML边界条件加载到Cole-Cole频散介质GPR有限元方程组中,并给出了详细的求解公式.在此基础上,编制了基于Pade逼近的Cole-Cole频散介质GPR有限元正演程序,利用该程序对均匀模型进行计算,并与解析解进行对比,验证了本文构建的GPR有限元正演算法的正确性和有效性.设计了一个复杂Cole-Cole频散介质GPR模型,利用本文构建的GPR有限元正演算法进行模拟并与非频散介质模型的模拟结果进行对比,分析了电磁波在Cole-Cole频散介质中传播衰减增强、子波延伸,分辨率降低等传播规律,有助于实测雷达资料更可靠、更准确的解释.模拟结果表明,基于Pade逼近的GPR有限元正演算法可用于复杂Cole-Cole频散介质结构模拟,且具有较高的计算精度. 相似文献
3.
沥青混凝土是由骨料、沥青胶浆、空气按照一定的体积百分比混合而成的多相非匀质混合物,其骨料、沥青胶浆和空气的体积不等、形状各异、介电特性不同、空间位置随机分布,具有明显的多相、离散、随机介质特征.本文基于随机介质模型理论,(1)测量与统计了介电常数在典型沥青混凝土芯样空间上的随机分布统计特征;(2)估算了沥青混凝土介质的自相关函数及其特征参数(自相关长度、自相关角度等),确定其随机介质类型;(3)提出了量化约束下的多相离散随机介质建模算法,以混合型椭圆自相关函数为基础,构建了不同粗糙度因子的多相离散随机介质模型;(4)构建了不同空隙率的多相离散随机介质模型,正演模拟与对比分析了探地雷达波在均匀介质、连续型随机介质和多相离散随机介质中的传播特征.结果表明:多相离散随机介质模型不仅描述了沥青混凝土的多相、离散与空间随机分布统计特征,而且进一步描述了其各组成物质体积百分比,能更全面、准确地描述沥青混凝土的介质特征,同时也为描述其他类似材料或介质提供了新的方法和途径;在多相离散随机介质模型中,探地雷达波散射强烈,随机、无序传播的散射波相互叠加干涉,形成了明显的随机扰动和"噪声",致使异常体反射波扭曲变形、不连续,降低了探地雷达回波的信噪比和分辨率.研究探地雷达波的随机扰动特征与多相离散随机介质模型参数之间的关系,将为定量评价多相离散随机介质的属性参数提供参考和帮助. 相似文献
4.
5.
6.
Electromagnetic geophysical methods, such as ground-penetrating radar (GPR), have proved to be optimal tools for detecting and mapping near-surface contaminants. GPR has the capability of mapping the location of hydrocarbon pools on the basis of contrasts in the effective permittivity and conductivity of the subsoil. At radar frequencies (50 MHz to 1 GHz), hydrocarbons have a relative permittivity ranging from 2 to 30, compared with a permittivity for water of 80. Moreover, their conductivity ranges from zero to 10 mS/m, against values of 200 mS/m and more for salt water. These differences indicate that water/hydrocarbon interfaces in a porous medium are electromagnetically 'visible'. In order to quantify the hydrocarbon saturation we developed a model for the electromagnetic properties of a subsoil composed of sand and clay/silt, and partially saturated with air, water and hydrocarbon. A self-similar theory is used for the sandy component and a transversely isotropic constitutive equation for the shaly component, which is assumed to possess a laminated structure. The model is first verified with experimental data and then used to obtain the properties of soils partially saturated with methanol and aviation gasoline. Finally, a GPR forward-modelling method computes the radargrams of a typical hydrocarbon spill, illustrating the sensitivity of the technique to the type of pore-fluid. The model and the simulation algorithm provide an interpretation methodology to distinguish different pore-fluids and to quantify their degree of saturation. 相似文献
7.
The objective of this work is to assess the importance of electrical conductivity and magnetic permeability variations in ground penetrating radar (GPR) reflections commonly interpreted only in terms of permittivity variations. We use the matrix propagator approach to obtain the surface electric field associated with a horizontally layered model whose three electromagnetic properties vary from layer to layer. The solution is based on the plane wave boundary value problem using inverse Fourier transformation to accommodate particular GPR pulses. Our results indicate that while magnetic permeability is unimportant, reflections from electrical conductivity variations can be of the same order as those associated with electrical permittivity boundaries. In particular, we show that a realistic ground model composed of thin conductive layers can produce radargrams similar to those caused by a lossless permittivity contrast. 相似文献
8.
9.
随机弹性介质中地震波散射衰减分析(英文) 总被引：2，自引：0，他引：2
地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。 相似文献
10.
基于可分离小波理论,由一维区间B样条小波尺度函数的张量积构造二维B样条小波基,并将它作为GPR波动方程求解的插值函数,通过引入转换矩阵,实现小波系数空间与雷达电磁场之间的转换.应用Galerkin算法,推导了二维区间B样条小波有限元GPR波动方程离散格式,求出了2阶1尺度与2阶2尺度BSWI尺度函数的积分值及联系系数,给出了该算法的详细求解过程.编制了BSWI的Matlab模拟程序,应用该程序对两个典型实例进行了正演,结果表明:BSWI能采用较少的单元达到与FEM相似的精度,而BSWI算法尺度提升能提高解的精度,但耗时会急剧增加.最后,将BSWI算法应用于双相随机混凝土模型,说明随机介质模型理论能灵活、有效地描述实际混凝土介质的分布,正演剖面与实测剖面特征更相符,能更真实地模拟雷达波的传播过程,可为提高GPR的探测效果和解释准确性提供理论基础. 相似文献
11.
In adapting the prestack migration technique used in seismic imaging to the inversion of ground‐penetrating radar (GPR) from time‐ to depth‐sections, we show that the theoretical integral formulation of the inversion can be applied to electromagnetic problems, albeit with three assumptions. The first two assumptions concern the electromagnetic characteristics of the medium, primarily that the medium must be perfectly resistive and non‐dispersive, and the third concerns the antennae radiation pattern, which is taken to be 2D. The application of this adaptation of the inversion method is confirmed by migrating actual GPR measurements acquired on the test site of the Laboratoire Central des Ponts et Chaussées. The results show good agreement with the geometry of the structures in the medium and confirm that the possible departure from the assumption of a purely resistive medium has no visible effect on the information concerning the geometry of scattering and reflecting structures. The field experiments also show that prestack migration processing is sufficiently robust with regard to the assumption of a non‐dispersive medium. The assumption of a 2D antennae radiation pattern, however, produces artefacts that could be significant for laterally heterogeneous media. Nevertheless, where the medium is not highly laterally heterogeneous, the migration gives a clear image of the scattering potential due to the geometry of structural contrasts in the medium; the scatterers are well focused from diffraction hyperbolae and well localized. Spatial geometry has limited dimensional accuracy and positions are located with a maximum error equal to the minimum wavelength of the signal bandpass. Objects smaller than one wavelength can nevertheless be detected and well focused if their dielectric contrasts are sufficiently high, as in the case of iron or water in gneiss gravels. Furthermore, the suitability of multi‐offset protocols to estimate the electromagnetic propagating velocity and to decrease the non‐coherent noise level of measurements is confirmed. Our velocity estimation is based on the semblance calculation of multi‐offset migrated images, and we confirmed the relevance of this quantification method using numerical data. The signal‐to‐noise ratio is improved by summing multi‐offset results after the addition of random noise on measurements. Thus the adaptation of prestack migration to multi‐offset radar measurements significantly improves the resolution of the scattering potential of the medium. Limitations associated with the methods used here suggest that 3D algorithms should be applied to strongly laterally heterogeneous media and further studies concerning the waveform inversion are necessary to obtain information about the electric nature of the medium. 相似文献
12.
In rough geologic media such as alluvial gravels, glacial tills, talus or colluvium, the grain sizes may span the range of GPR in situ wavelengths. Here we experimentally and numerically modeled the scattering loss from both rough-surface and subsurface dielectric scatterers. The combination of the selected radar frequency and the dimension of the scatterers placed the scattering within the Mie regime. We compared the GPR signal amplitude and waveform reflected from the metal sheet on the bottom of a large box filled with boulders with the numerically computed response from a discrete random medium (DRM) model. The DRM consists of a collection of densely packed ellipsoids. The size and orientation of the ellipsoids are randomized; the size has a Gaussian distribution similar to the physical experiment. The dielectric permittivity of the ellipsoids is constant and their electric conductivity is negligible. The starting in situ dominant pulse wavelength at 900 MHz was about 17 cm, as was about the average rock dimension. Experimentally, the 900-MHz radar pulse underwent most dispersion within the first in situ wavelength of depth, and then, at 500–700 MHz dominant frequency, the pulses underwent a near inverse range dependency loss rate, as if the media were a pure dielectric. The numerical model agrees well with the experimental data. Both experimental and numerical results support a significant scattering loss in Mie regime. Besides the scattering attenuation loss, velocity dispersion has also been observed from both observation and simulation. However, the scattering attenuation and dispersion cannot be fit by the Kramers–Kronig relation that is commonly found in intrinsic attenuation and worth further theoretical investigations. 相似文献
13.
机载探地雷达可以用于人类无法到达的危险地区、植被严重覆盖的地下目标体探测,然而由于机载探地雷达的特殊性,影响机载探地雷达探测效果的因素包括天线的极化方向、天线的飞行高度以及地表粗糙度等.为了研究这些影响因素与探测效果之间的关系,用三维时间域有限差分模拟电磁波的传播过程,以沙漠地区地下空洞掩体的机载探地雷达探测为实例,分别模拟了不同天线极化方向、天线高度及地表粗糙度情况下的机载探地雷达剖面,分析了各因素对机载探地雷达探测地下空洞目标体的影响.天线极化方向与目标体走向垂直更有利于地下目标体探测;天线距离地表越近,可以获得更高分辨率的雷达剖面;沙漠地表起伏越大,雷达剖面中的散射杂波能量越强,浅部地下目标体信号容易被掩盖.为了消除起伏地形造成的散射杂波,提出用逆时偏移成像技术对共炮集机载探地雷达数据进行偏移成像,成像结果优于基尔霍夫偏移成像结果. 相似文献
14.
Maksim Bano 《Geophysical Prospecting》2004,52(1):11-26
The attenuation of ground‐penetrating radar (GPR) energy in the subsurface decreases and shifts the amplitude spectrum of the radar pulse to lower frequencies (absorption) with increasing traveltime and causes also a distortion of wavelet phase (dispersion). The attenuation is often expressed by the quality factor Q. For GPR studies, Q can be estimated from the ratio of the real part to the imaginary part of the dielectric permittivity. We consider a complex power function of frequency for the dielectric permittivity, and show that this dielectric response corresponds to a frequency‐independent‐Q or simply a constant‐Q model. The phase velocity (dispersion relationship) and the absorption coefficient of electromagnetic waves also obey a frequency power law. This approach is easy to use in the frequency domain and the wave propagation can be described by two parameters only, for example Q and the phase velocity at an arbitrary reference frequency. This simplicity makes it practical for any inversion technique. Furthermore, by using the Hilbert transform relating the velocity and the absorption coefficient (which obeys a frequency power law), we find the same dispersion relationship for the phase velocity. Both approaches are valid for a constant value of Q over a restricted frequency‐bandwidth, and are applicable in a material that is assumed to have no instantaneous dielectric response. Many GPR profiles acquired in a dry aeolian environment have shown a strong reflectivity inside dunes. Changes in water content are believed to be the origin of this reflectivity. We model the radar reflections from the bottom of a dry aeolian dune using the 1D wavelet modelling method. We discuss the choice of the reference wavelet in this modelling approach. A trial‐and‐error match of modelled and observed data was performed to estimate the optimum set of parameters characterizing the materials composing the site. Additionally, by combining the complex refractive index method (CRIM) and/or Topp equations for the bulk permittivity (dielectric constant) of moist sandy soils with a frequency power law for the dielectric response, we introduce them into the expression for the reflection coefficient. Using this method, we can estimate the water content and explain its effect on the reflection coefficient and on wavelet modelling. 相似文献
15.
地质雷达是利用电磁波对地下不同电性介质进行探测的地球物理仪器，其探测速率快、分辨率高，可弥补探槽和其他地球物理方法存在探测盲区的缺陷，正在越来越多地应用于活动断层探测领域。本文以乌拉山山前断裂为例开展地质雷达探测工作，使用无人机正射影像技术对测线进行地形校正，获得断层浅部地质雷达图像。研究结果表明，本文研究方法能有效反映探槽揭露的地层单元和断层分布。本次探测中，雷达波形图像特征为:浅地表的土壤层反射波总体较弱；粗粒沉积为主的砾石层反射波总体较强，同相轴连续性好；细粒沉积为主的砂层反射波弱于砾石层，波形以中、高频为主，同相轴具有弱连续性；对于洪冲积地区，地质雷达能分辨具有一定特征的地层单元，这为剖面图像的断层识别提供了标志；通过无人机正射影像技术对地质雷达测线进行地形校正，有利于获得更为准确的探测结果。 相似文献
16.
In this study, we focus on a hydrogeological inverse problem specifically targeting monitoring soil moisture variations using tomographic ground penetrating radar (GPR) travel time data. Technical challenges exist in the inversion of GPR tomographic data for handling non-uniqueness, nonlinearity and high-dimensionality of unknowns. We have developed a new method for estimating soil moisture fields from crosshole GPR data. It uses a pilot-point method to provide a low-dimensional representation of the relative dielectric permittivity field of the soil, which is the primary object of inference: the field can be converted to soil moisture using a petrophysical model. We integrate a multi-chain Markov chain Monte Carlo (MCMC)–Bayesian inversion framework with the pilot point concept, a curved-ray GPR travel time model, and a sequential Gaussian simulation algorithm, for estimating the dielectric permittivity at pilot point locations distributed within the tomogram, as well as the corresponding geostatistical parameters (i.e., spatial correlation range). We infer the dielectric permittivity as a probability density function, thus capturing the uncertainty in the inference. The multi-chain MCMC enables addressing high-dimensional inverse problems as required in the inversion setup. The method is scalable in terms of number of chains and processors, and is useful for computationally demanding Bayesian model calibration in scientific and engineering problems. The proposed inversion approach can successfully approximate the posterior density distributions of the pilot points, and capture the true values. The computational efficiency, accuracy, and convergence behaviors of the inversion approach were also systematically evaluated, by comparing the inversion results obtained with different levels of noises in the observations, increased observational data, as well as increased number of pilot points. 相似文献
17.
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 (f–k) 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. 相似文献
18.
探地雷达图像的正演合成与偏移处理 总被引：8，自引：2，他引：8
本文讨论了射线追踪法实现二维地电断面的探地雷达图像的正演合成问题,其中包括二维地电断面的数字拟合、多层连续界面介质中高频电磁波的反射与透射统一的运动学边界方程;根据Huygen's-Fresnel原理导出了均匀介质中二维管状体模型的衍射雷达波的计算公式;给出了用射线追踪法合成的二维地电模型的探地雷达图像;用衍射波公式合成了二维管状模型的雷达图像,并与物理模拟作了比较.在此基础上,用Kirchhoff积分法对数字模拟与物理模拟记录的雷达图像作了偏移处理,并给出了实例. 相似文献
19.
The dielectric permittivity of fiery spoil tips (Shakhty town, Rostov Region) is studied with the use of a GROT 12E remote-controlled ground-penetrating radar (GPR). An anomalous zone in a combustion source is shown to be clearly pronounced in GPR data due to the temperature dependence of the dielectric permittivity of these spoil tips. To substantiate this statement, the GPR data are compared with direct measurements of soil temperatures at depths from 1.5 to 2.5 m. The experimental results are compared with the variable spectral range of a GPR sounding pulse. GPR is shown to be a promising tool for the mapping of temperature-contrast underground objects. 相似文献
20.
应用小波伽略金方法,对Maxwell方程进行离散化,导出了DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件.在此基础上,开发了探地雷达MRTD(multi-resolution time domain)法正演模拟程序,该程序极大地提高了运算速度,改善了三维探地雷达正演方法,并利用该自制程序,对三角形金属体模型进行了正演模拟,得到了其相应的正演合成三维剖视图及切片图,通过对这些模拟结果进行分析,可以加深对三维雷达反射特征的认识,提高探地雷达探测的可靠性、准确度,同时也说明时域多分辨率法在探地雷达三维正演模拟中的有效性. 相似文献