共查询到20条相似文献，搜索用时 936 毫秒
1.
Joint inversion of long‐offset and central‐loop transient electromagnetic data: Application to a mud volcano exploration in Perekishkul,Azerbaijan 
下载免费PDF全文
下载免费PDF全文 J. Adrian R. Bergers M. Gurk B. Tezkan A. L. Mammadov A. G. Novruzov 《Geophysical Prospecting》2015,63(2):478-494
Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly. 相似文献
2.
Effects of the sea floor topography on the 1D inversion of time‐domain marine controlled source electromagnetic data 下载免费PDF全文
下载免费PDF全文 Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets. 相似文献
3.
Sufficient low‐frequency information is essential for full‐waveform inversion to get the global optimal solution. Multi‐scale envelope inversion was proposed using a new Fréchet derivative to invert the long‐wavelength component of the model by directly using the low‐frequency components contained in an envelope of seismic data. Although the new method can recover the main structure of the model, the inversion quality of the model bottom still needs to be improved. Reflection waveform inversion reduces the dependence of inversion on low‐frequency and long‐offset data by using travel‐time information in reflected waves. However, when the underground medium contains strong contrast or the initial model is far away from the true model, it is hard to get reliable reference reflectors for the generation of reflected waves. Here, we propose a combination inversion algorithm, i.e., reflection multi‐scale envelope inversion, to overcome the limitations of multi‐scale envelope inversion and reflection waveform inversion. First, wavefield decomposition was introduced into the multi‐scale envelope inversion to improve the inversion quality of the long‐wavelength components of the model. Then, after the initial model had been established to be accurate enough, migration and de‐migration were introduced to achieve multi‐scale reflection waveform inversion. The numerical results of the salt‐layer model and the SEG/EAGE salt model verified the validity of the proposed approach and its potential. 相似文献
4.
Migration velocity analysis and waveform inversion 总被引：3，自引：0，他引：3
William W. Symes 《Geophysical Prospecting》2008,56(6):765-790
Least‐squares inversion of seismic reflection waveform data can reconstruct remarkably detailed models of subsurface structure and take into account essentially any physics of seismic wave propagation that can be modelled. However, the waveform inversion objective has many spurious local minima, hence convergence of descent methods (mandatory because of problem size) to useful Earth models requires accurate initial estimates of long‐scale velocity structure. Migration velocity analysis, on the other hand, is capable of correcting substantially erroneous initial estimates of velocity at long scales. Migration velocity analysis is based on prestack depth migration, which is in turn based on linearized acoustic modelling (Born or single‐scattering approximation). Two major variants of prestack depth migration, using binning of surface data and Claerbout's survey‐sinking concept respectively, are in widespread use. Each type of prestack migration produces an image volume depending on redundant parameters and supplies a condition on the image volume, which expresses consistency between data and velocity model and is hence a basis for velocity analysis. The survey‐sinking (depth‐oriented) approach to prestack migration is less subject to kinematic artefacts than is the binning‐based (surface‐oriented) approach. Because kinematic artefacts strongly violate the consistency or semblance conditions, this observation suggests that velocity analysis based on depth‐oriented prestack migration may be more appropriate in kinematically complex areas. Appropriate choice of objective (differential semblance) turns either form of migration velocity analysis into an optimization problem, for which Newton‐like methods exhibit little tendency to stagnate at nonglobal minima. The extended modelling concept links migration velocity analysis to the apparently unrelated waveform inversion approach to estimation of Earth structure: from this point of view, migration velocity analysis is a solution method for the linearized waveform inversion problem. Extended modelling also provides a basis for a nonlinear generalization of migration velocity analysis. Preliminary numerical evidence suggests a new approach to nonlinear waveform inversion, which may combine the global convergence of velocity analysis with the physical fidelity of model‐based data fitting. 相似文献
5.
A local plane-wave approach of generalized diffraction tomography in heterogeneous backgrounds, equivalent to Kirchhoff summation techniques when applied in seismic reflection, is re-programmed to act as repeated synthetic aperture radar (SAR) imaging for seismic prestack depth migration. Spotlight-mode SAR imaging quickly provides good images of the electromagnetic reflectivity of the ground via fast Fourier transform (FFT)-based signal processing. By calculating only the Green's functions connecting the aircraft to the centre of the illuminated patch, scattering structures around that centre are also recovered. SAR technology requires us to examine seismic imaging from the local point of view, where the quantity and quality of the available information at each image point are what are important, regardless of the survey geometry. When adapted to seismics, a local image of arbitrary size and sampling is obtained by FFT of seismic energy maps in the scattering wavenumber domain around each node of a pre-calculated grid of Green's functions. These local images can be used to generate a classic prestack depth-migrated section by collecting only their centres. However, the local images also provide valuable information around the centre, as in SAR. They can therefore help to pre-analyse prestack depth migration efficiently, and to perform velocity analysis at a very low cost. The FFT-based signal-processing approach allows local, efficient and automatic control of anti-aliasing, noise and resolution, including optimized Jacobian weights. Repeated local imaging could also be used to speed up migration, with interpolation between local images associated with a coarse grid of Green's functions, as an alternative to interpolation of Green's functions. The local images may, however, show distortions due to the local plane-wave approximation, and the velocity variations across their frame. Such effects, which are not necessarily a problem in SAR, should be controlled and corrected to further enhance seismic imaging. Applications to realistic models and to real data show that, despite the distortion effects, the local images can yield similar information to prestack depth migration, including common-image-point gathers for velocity analyses and AVO/AVA effects, at a much lower cost when a small target is considered. 相似文献
6.
地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具. 相似文献
7.
Scattering theory, a form of perturbation theory, is a framework from within which time‐lapse seismic reflection methods can be derived and understood. It leads to expressions relating baseline and monitoring data and Earth properties, focusing on differences between these quantities as it does so. The baseline medium is, in the language of scattering theory, the reference medium and the monitoring medium is the perturbed medium. The general scattering relationship between monitoring data, baseline data, and time‐lapse Earth property changes is likely too complex to be tractable. However, there are special cases that can be analysed for physical insight. Two of these cases coincide with recognizable areas of applied reflection seismology: amplitude versus offset modelling/inversion, and imaging. The main result of this paper is a demonstration that time‐lapse difference amplitude versus offset modelling, and time‐lapse difference data imaging, emerge from a single theoretical framework. The time‐lapse amplitude versus offset case is considered first. We constrain the general time‐lapse scattering problem to correspond with a single immobile interface that separates a static overburden from a target medium whose properties undergo time‐lapse changes. The scattering solutions contain difference‐amplitude versus offset expressions that (although presently acoustic) resemble the expressions of Landro ( 2001 ). In addition, however, they contain non‐linear corrective terms whose importance becomes significant as the contrasts across the interface grow. The difference‐amplitude versus offset case is exemplified with two parameter acoustic (bulk modulus and density) and anacoustic (P‐wave velocity and quality factor Q) examples. The time‐lapse difference data imaging case is considered next. Instead of constraining the structure of the Earth volume as in the amplitude versus offset case, we instead make a small‐contrast assumption, namely that the time‐lapse variations are small enough that we may disregard contributions from beyond first order. An initial analysis, in which the case of a single mobile boundary is examined in 1D, justifies the use of a particular imaging algorithm applied directly to difference data shot records. This algorithm, a least‐squares, shot‐profile imaging method, is additionally capable of supporting a range of regularization techniques. Synthetic examples verify the applicability of linearized imaging methods of the difference image formation under ideal conditions. 相似文献
8.
In this case study we consider the seismic processing of a challenging land data set from the Arabian Peninsula. It suffers from rough top‐surface topography, a strongly varying weathering layer, and complex near‐surface geology. We aim at establishing a new seismic imaging workflow, well‐suited to these specific problems of land data processing. This workflow is based on the common‐reflection‐surface stack for topography, a generalized high‐density velocity analysis and stacking process. It is applied in a non‐interactive manner and provides an entire set of physically interpretable stacking parameters that include and complement the conventional stacking velocity. The implementation introduced combines two different approaches to topography handling to minimize the computational effort: after initial values of the stacking parameters are determined for a smoothly curved floating datum using conventional elevation statics, the final stack and also the related residual static correction are applied to the original prestack data, considering the true source and receiver elevations without the assumption of nearly vertical rays. Finally, we extrapolate all results to a chosen planar reference level using the stacking parameters. This redatuming procedure removes the influence of the rough measurement surface and provides standardized input for interpretation, tomographic velocity model determination, and post‐stack depth migration. The methodology of the residual static correction employed and the details of its application to this data example are discussed in a separate paper in this issue. In view of the complex near‐surface conditions, the imaging workflow that is conducted, i.e. stack – residual static correction – redatuming – tomographic inversion – prestack and post‐stack depth migration, leads to a significant improvement in resolution, signal‐to‐noise ratio and reflector continuity. 相似文献
9.
10.
Double-square-root one-way wave equation prestack tau migration in heterogeneous media 总被引：1，自引：0，他引：1
In this paper, source‐receiver migration based on the double‐square‐root one‐way wave equation is modified to operate in the two‐way vertical traveltime (τ) domain. This tau migration method includes reasonable treatment for media with lateral inhomogeneity. It is implemented by recursive wavefield extrapolation with a frequency‐wavenumber domain phase shift in a constant background medium, followed by a phase correction in the frequency‐space domain, which accommodates moderate lateral velocity variations. More advanced τ‐domain double‐square‐root wave propagators have been conceptually discussed in this paper for migration in media with stronger lateral velocity variations. To address the problems that the full 3D double‐square‐root equation prestack tau migration could meet in practical applications, we present a method for downward continuing common‐azimuth data, which is based on a stationary‐phase approximation of the full 3D migration operator in the theoretical frame of prestack tau migration of cross‐line constant offset data. Migrations of synthetic data sets show that our tau migration approach has good performance in strong contrast media. The real data example demonstrates that common‐azimuth prestack tau migration has improved the delineation of the geological structures and stratigraphic configurations in a complex fault area. Prestack tau migration has some inherent robust characteristics usually associated with prestack time migration. It follows a velocity‐independent anti‐aliasing criterion that generally leads to reduction of the computation cost for typical vertical velocity variations. Moreover, this τ‐domain source‐receiver migration method has features that could be of help to speed up the convergence of the velocity estimation. 相似文献
11.
Tijmen Jan Moser 《Geophysical Prospecting》2014,62(3):475-490
We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots. 相似文献
12.
Anwar H. Bhuyian Brd P. Thrane Martin Landr Stle E. Johansen 《Geophysical Prospecting》2010,58(3):505-533
A comprehensive controlled source electromagnetic (CSEM) modelling study, based on complex resistivity structures in a deep marine geological setting, is conducted. The study demonstrates the effects of acquisition parameters and multi‐layered resistors on CSEM responses. Three‐dimensional (3D) finite difference time domain (FDTD) grid‐modelling is used for CSEM sensitivity analysis. Interpolation of normalized CSEM responses provides attributes representing relative sensitivity of the modelled structures. Modelling results show that fine grid, 1 × 1 km receiver spacing, provides good correlations between CSEM responses and the modelled structures, irrespective of source orientation. The resolution of CSEM attributes decreases for receiver spacing >2 × 2 km, when using only in‐line data. Broadside data in the grid geometry increase data density by 100 – approximately 200% by filling in in‐line responses and improve the resolution of CSEM attributes. Optimized source orientation (i.e., oblique to the strike of an elongated resistor) improves the structural definition of the CSEM anomalies for coarse‐grid geometries (receiver spacing ≥3 × 3 km). The study also shows that a multi‐resistor anomaly is not simply the summation but a cumulative response with mutual interference between constituent resistors. The combined response of constituent resistors is approximately 50% higher than the cumulative response of the multi‐resistor for 0.5 Hz at 4000 m offset. A gradual inverse variation of offset and frequency allows differentiation of CSEM anomalies for multi‐layered resistors. Similar frequency‐offset variations for laterally persistent high‐resistivity facies show visual continuity with varying geometric expressions. 3D grid‐modelling is an effective and adequate tool for CSEM survey design and sensitivity analysis. 相似文献
13.
Signal detectability of marine electromagnetic methods in the exploration of resistive targets 下载免费PDF全文
下载免费PDF全文 Mark Goldman Vladimir Mogilatov Amir Haroon Eldad Levi Bülent Tezkan 《Geophysical Prospecting》2015,63(1):192-210
We compare selected marine electromagnetic methods for sensitivity to the presence of relatively thin resistive targets (e.g., hydrocarbons, gas hydrates, fresh groundwater, etc.). The study includes the conventional controlled‐source electromagnetic method, the recently introduced transient electromagnetic prospecting with vertical electric lines method, and the novel marine circular electric dipole method, which is still in the stage of theoretical development. The comparison is based on general physical considerations, analytical (mainly asymptotic) analysis, and rigorous one‐dimensional and multidimensional forward modelling. It is shown that transient electromagnetic prospecting with vertical electric lines and marine circular electric dipole methods represent an alternative to the conventional controlled‐source electromagnetic method at shallow sea, where the latter becomes less efficient due to the air‐wave phenomenon. Since both former methods are essentially short‐offset time‐domain techniques, they exhibit a much better lateral resolution than the controlled‐source electromagnetic method in both shallow sea and deep sea. The greatest shortcoming of the transient electromagnetic prospecting with vertical electric lines and marine circular electric dipole methods comes from the difficulties in accurately assembling the transmitter antenna within the marine environment. This makes these methods significantly less practical than the controlled‐source electromagnetic method. Consequently, the controlled‐source electromagnetic method remains the leading marine electromagnetic technique in the exploration of large resistive targets in deep sea. However, exploring laterally small targets in deep sea and both small and large targets in shallow sea might require the use of the less practical transient electromagnetic prospecting with vertical electric lines and/or marine circular electric dipole method as a desirable alternative to the controlled‐source electromagnetic method. 相似文献
14.
2D multiscale non-linear velocity inversion 总被引：2，自引：0，他引：2
An efficient and robust non-linear inversion method for velocity optimization combining a global random search followed by a simplex technique is presented. The background velocity field is estimated at different spatial scales by analysing image gathers after iterative prestack depth migrations. First, the global random search is used to determine the main features/trends of the velocity model (large-scale component). Then, the simplex technique improves the resolution of the velocity field by estimating smaller-scale features. A measure of the quality of the velocity model (objective function) is based on flattening offset events in depth-migrated image gathers. To help constrain the solution, the algorithm can incorporate a priori information about the model and a smoothness condition. This 2D velocity estimation offers the benefit of being semi-automatic (requiring minimal human intervention) as well as providing a global and objective solution (which is a useful approach to an interpretation-derived velocity-estimation technique). The method is applied to a real data set where AVO analysis is carried out after prestack depth migration, as structural effects are non-negligible. It is demonstrated that the method can successfully estimate a laterally inhomogeneous velocity model at a computational cost modest compared with an interpretation-based iterative prestack depth velocity-analysis technique. 相似文献
15.
2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D . 相似文献
16.
拟声波最小二乘逆时偏移是一种极具潜力的地震波成像工具,但该方法遭受各向异性拟声波近似的限制,TTI介质正演模拟不稳定、反偏移记录中遭受伪横波二次扰动及数值频散假象,另外拟声波最小二乘逆时偏移还面临计算效率低、收敛速度慢、对速度等模型参数依赖性高等问题.为了克服各向异性拟声波最小二乘逆时偏移的缺陷,在反演框架下,本文借助Low-rank有限差分算法首次提出并实现了TTI介质纯qP波线性正演模拟及纯qP波最小二乘逆时偏移;为了进一步提升反演成像效率,同时改善反演成像方法对模型参数误差的依赖性及对地震数据噪声的适应性,通过引入叠前平面波优化策略,发展了TTI介质纯qP波叠前平面波最小二乘逆时偏移成像方法.在编程实现方法的基础上,通过开展模型成像测试,展示了本方法的优势和潜力：一方面加快了反演成像效率,另一方面也提升了方法的抗噪性,同时还降低了方法对模型参数的依赖性. 相似文献
17.
18.
An imaging technique is developed which exceeds the resolution limitation prescribed by conventional seismic imaging methods. The high‐resolution imaging is obtained by introducing a sparseness‐constrained least‐squares inversion into the imaging process of prestack depth migration. This is implemented by a proposed interference technique. In contrast to conventional depth migration, a decomposed signal or combined event, instead of the source wavelet, is needed in the proposed scheme. The proposed method aims to image a small local region with a higher resolution using the prestack data set. It should be applied following conventional depth imaging if a higher resolution is needed in a target zone rather than replacing the conventional depth imaging for the entire medium. Synthetic examples demonstrate the significant improvements in the resolution using the proposed scheme. 相似文献
19.
20.
Target detection using ground penetrating radar (GPR) is based on the contrast between the electrical parameters of the target and the background medium, such as dielectric permittivity, conductivity and permeability. The application mainly concentrates on the detection of the medium interface and the target shape. In any theoretical study, a simulation model is built with a homogeneous medium. However, real detection encounters heterogeneous media which might produce scattering and diffraction at electrical interfaces and distort the radar pulse shape and affect the detection resolution. In this paper, we build multi-scale random media model with an ellipsoidal autocorrelation function and use FDTD method to simulate the GPR signal response. We then estimate and analyze the arrival time, layer thickness, permittivity and the physics relation in different scale random models according to the S transform method and the transmission wave method. The results demonstrate that we can use GPR to obtain geophysical information of multi-scale heterogeneous media, and provide a foundation for real media detection and complex media inversion. 相似文献