共查询到20条相似文献，搜索用时 31 毫秒
1.
Martin Musil Hansruedi Maurer Klaus Hollinger Alan G. Green 《Geophysical Prospecting》2006,54(3):273-285
Rapid melting of permafrost in many alpine areas has increased the probability of catastrophic rock slides. In an attempt to provide critical structural information needed for the design and implementation of suitable mitigation procedures, we have acquired low frequency (22 MHz) cross‐hole radar data from within a fast‐moving rock glacier, an important form of alpine permafrost. Since the ice, rock and pockets of water and air found in the underground of high alpine areas have very different dielectric permittivities and electrical conductivities, the radar method was well‐suited for investigating the structure and state of the rock glacier. Our interpretation of the radar velocities and attenuations was constrained by geomorphological observations, borehole lithological logs and the results of a surface seismic survey. The radar data revealed the existence of a discontinuous 7–11 m thick ice‐rich zone distinguished by high velocities (0.14–0.17 m/ns) and low attenuations (0.04–0.09 m?1) and a thin underlying ice‐free zone characterized by moderate velocities (0.11–0.12 m/ns) and low attenuations (0.04–0.09 m?1). Beneath these two zones, we observed a prominent band of high velocities (0.14–0.17 m/ns) and moderately high attenuations (0.10–0.20 m?1) associated with unconsolidated glacial sediments and numerous large air‐filled voids, which in the past were probably filled with ice. At greater depths, the variably dry to water‐saturated sediments were represented by generally lower velocities (0.08–0.10 m/ns) and higher attenuations (0.16–0.24 m?1). The bedrock surface was represented by an abrupt ～0.03 m/ns velocity increase. We speculate that the disappearance of ice, both laterally and with depth, occurred during the past one to two decades. 相似文献
2.
The multifold acquisition principle was applied to a borehole radar survey, performed in a granitic site (Grimsel Test Site, Switzerland). Two multifold coverage acquisitions (40-fold and 20-fold) were carried out in a subhorizontal borehole. Instrumental drifts (transmission time and sampling frequency fluctuations) were corrected in order to remove shifts observed on CMP gathers and to optimize velocity analysis and trace stacking. Computation of velocity spectra was adapted in order to take into account the features of the medium investigated (homogeneous velocity, various reflector orientations). The NMO velocities were then interpreted as angles between reflectors and the survey line. The processing, based on the computation of several constant velocity stacked sections performed with different NMO velocities, leads to better results than the standard DMO + NMO processing. The signal-to-noise ratio of the stacked profile is improved in comparison with the single-fold section, which results from a standard acquisition. From a practical point of view, the implementation of a multifold radar survey within a borehole is difficult but a greater investigation range is obtained, more reflectors are detected and the mapping of geological discontinuities is improved. 相似文献
3.
Previous studies of ground ice using moveout type ground-penetrating radar (GPR) surveys indicate that the dielectric permittivity can constrain the type of ground ice present in the subsurface. Due to the high-loss nature of the active layer over permafrost targets, however, the signal strength of GPR signals is often insufficient to resolve the basal boundary required for determining the dielectric permittivity of an underlying unit. We apply a non-conventional antenna orientation and post-processing method to determine the dielectric permittivity of the unit underlying the lowest resolvable boundary. We conduct moveout surveys using a 450 MHz GPR with collinear parallel oriented antennas on two adjacent ground ice formations in the region of Thomas Lee Inlet, Devon Island, Nunavut. We exploit the Brewster angle to calculate the approximate dielectric permittivity of ground ice formations below the active layer. The results agree within 1 dielectric unit with on-ice permittivity measurements made during a complementary study of the site. 相似文献
4.
It is well established that the upper—cooler—part of the crust is brittle, while deeper zones present ductile behaviour. In some cases, this brittle–ductile transition is a single seismic reflector with an associated reflection coefficient. We first develop a stress–strain relation including the effects of crust anisotropy, seismic attenuation and ductility in which deformation takes place by shear plastic flow. Viscoelastic anisotropy is based on the eigenstrain model and the Zener and Burgers mechanical models are used to model the effects of seismic attenuation, velocity dispersion, and steady-state creep flow, respectively. The stiffness components of the brittle and ductile media depend on stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P- and S-wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the reflection and transmission coefficients of a single brittle–ductile interface and of a ductile thin layer. The PP scattering coefficient has a Brewster angle (a sign change) in both cases, and there is substantial PS conversion at intermediate angles. The PP coefficient is sensitive to the layer thickness, unlike the SS coefficient. Thick layers have a well-defined Brewster angle and show higher reflection amplitudes. Finally, we compute synthetic seismograms in a homogeneous medium as a function of temperature. 相似文献
5.
A statistical investigation of the relationship between VHF radar auroral backscatter intensity and Doppler velocity has been undertaken with data collected from 8 years operation of the Wick site of the Sweden And Britain Radar-auroral Experiment (SABRE). The results indicate three different regimes within the statistical data set; firstly, for Doppler velocities <200 m s−1, the backscatter intensity (measured in decibels) remains relatively constant. Secondly, a linear relationship is observed between the backscatter intensity (in decibels) and Doppler velocity for velocities between 200 m s−1 and 700 m s−1. At velocities greater than 700 m s−1 the backscatter intensity saturates at a maximum value as the Doppler velocity increases. There are three possible geophysical mechanisms for the saturation in the backscatter intensity at high phase speeds: a saturation in the irregularity turbulence level, a maximisation of the scattering volume, and a modification of the local ambient electron density. There is also a difference in the dependence of the backscatter intensity on Doppler velocity for the flow towards and away from the radar. The results for flow towards the radar exhibit a consistent relationship between backscatter intensity and measured velocities throughout the solar cycle. For flow away from the radar, however, the relationship between backscatter intensity and Doppler velocity varies during the solar cycle. The geometry of the SABRE system ensures that flow towards the radar is predominantly associated with the eastward electrojet, and flow away is associated with the westward electrojet. The difference in the backscatter intensity variation as a function of Doppler velocity is attributed to asymmetries between the eastward and westward electrojets and the geophysical parameters controlling the backscatter amplitude. 相似文献
6.
沥青混凝土是由骨料、沥青胶浆、空气按照一定的体积百分比混合而成的多相非匀质混合物,其骨料、沥青胶浆和空气的体积不等、形状各异、介电特性不同、空间位置随机分布,具有明显的多相、离散、随机介质特征.本文基于随机介质模型理论,(1)测量与统计了介电常数在典型沥青混凝土芯样空间上的随机分布统计特征;(2)估算了沥青混凝土介质的自相关函数及其特征参数(自相关长度、自相关角度等),确定其随机介质类型;(3)提出了量化约束下的多相离散随机介质建模算法,以混合型椭圆自相关函数为基础,构建了不同粗糙度因子的多相离散随机介质模型;(4)构建了不同空隙率的多相离散随机介质模型,正演模拟与对比分析了探地雷达波在均匀介质、连续型随机介质和多相离散随机介质中的传播特征.结果表明:多相离散随机介质模型不仅描述了沥青混凝土的多相、离散与空间随机分布统计特征,而且进一步描述了其各组成物质体积百分比,能更全面、准确地描述沥青混凝土的介质特征,同时也为描述其他类似材料或介质提供了新的方法和途径;在多相离散随机介质模型中,探地雷达波散射强烈,随机、无序传播的散射波相互叠加干涉,形成了明显的随机扰动和"噪声",致使异常体反射波扭曲变形、不连续,降低了探地雷达回波的信噪比和分辨率.研究探地雷达波的随机扰动特征与多相离散随机介质模型参数之间的关系,将为定量评价多相离散随机介质的属性参数提供参考和帮助. 相似文献
7.
Velocity and absorption tomograms are the two most common forms of presentation of radar tomographic data. However, mining personnel, geophysicists included, are often unfamiliar with radar velocity and absorption. In this paper, general formulae are introduced, relating velocity and attenuation coefficient to conductivity and dielectric constant. The formulae are valid for lossy media as well as high-resistivity materials. The transformation of velocity and absorption to conductivity and dielectric constant is illustrated via application of the formulae to radar tomograms from the Hellyer zinc–lead–silver mine, Tasmania, Australia. The resulting conductivity and dielectric constant tomograms constructed at Hellyer demonstrated the potential of radar tomography to delineate sulphide ore zones. 相似文献
8.
本文应用Chisnell- no方法,求解了在理想介貭中,垂直磁流体冲激波在非均勻磁場中的传播問題。这种方法,把非均勻介貭分解成无限小的弱間断面,根据气体动力学中波与間断面相互作用的原理,算出激波通过弱間断面时的强度变化,然后用积分求得激波通过整个非均勻区时的强度变化。作者引入了激波的特征速度（它是激波在波前后介貭中传播速度的几何平均值）作为輔助参量,得到形式上比較簡单的激波传播方程。然后考虑了磁压力远大于气体压力的强磁介貭中的激波传播問題,并进行了数值积分。采用的介貭密度模型有三种:（1）阿尔芬波速为常数;（2）密度不变;（3）密度与磁場强度成正比。計算結果表明:当激波由弱磁場向强磁場传播时,激波的强度逐漸变弱。其中,在阿尔芬波速为常数的介貭中,激波强度的衰减最为緩慢;在密度不变的介貭中,激波强度的衰減最为迅速;而在密度与磁場成正比的介貭中,激波强度的衰減則介乎上述两种密度分布之間。作者联系磁流体冲激波在地球外层空間的传播問題进行了討論,密度的模型采取大气啃昔的观測結果（卽上述第三种密度分布）,并进行了适当的外推,估計了在十个地球半径处的磁流体冲激波传到地面时的强度,求出了激波在地面引起的磁場变化与激波初始速度之間的关系。根据上述簡化模型,計算結果表明,在十个地球半径处初始速度为108厘米/秒的激波,传到地面引起的磁場变化大約为60伽（亻馬）,这个数值的量級恰好与中低緯度强磁暴的急始变幅相符。 相似文献
9.
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. 相似文献
10.
11.
《Journal of Applied Geophysics》1999,41(1):19-30
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. 相似文献
12.
The EISCAT UHF radar system was used to study the characteristics of E-region coherent backscatter at very large magnetic aspect angles (5–11°). Data taken using 60 s pulses during elevation scans through horizontally uniform backscatter permitted the use of inversion techniques to determine height profiles of the scattering layer. The layer was always singly peaked, with a mean height of 104 km, and mean thickness (full width at half maximum) of 10 km, both independent of aspect angle. Aspect sensitivities were also estimated, with the Sodankylä-Tromsø link observing 5 dB/degree at aspect angles near 5°, decreasing to 3 dB/degree at 10° aspect angle. Observed coherent phase velocities from all three stations were found to be roughly consistent with LOS measurements of a common E-region phase velocity vector. The E-region phase velocity had the same orientation as the F-region ion drift velocity, but was approximately 50% smaller in magnitude. Spectra were narrow with skewness of about –1 (for negative velocities), increasing slightly with aspect angle. 相似文献
13.
Ground-penetrating radar (GPR) is an effective tool for imaging the spatial distribution of water content. An artificial groundwater recharge test was conducted in Nagaoka City in Japan, and time-lapse crosshole GPR data were collected to monitor the infiltration process in the vadose zone. Since electromagnetic wave velocities in the vadose zone are largely controlled by variations in water content, an increase in traveltime is interpreted as an increase in saturation. In the test zone, the infiltrated water penetrated downward with an average velocity of about 2.7 m/h. A finite-difference time-domain method using two-dimensional cylindrical coordinates is applied to simulate radargrams associated with the advancing wetting front and to quantify the effects of critical refraction. Standard zero-offset profiling for which all first-arrivals are assumed to be direct waves results in an underestimation of water content in the transition zone above the wetting front. As a result, correct velocity analysis requires identification of first-arriving critically refracted waves from the traveltime profile to accurately determine a water content profile. 相似文献
14.
Emil Blias 《Geophysical Prospecting》2009,57(3):323-341
Lateral velocity changes (velocity anomalies) in the overburden may cause significant oscillations in normal moveout velocities. Explicit analytical moveout formulas are presented and provide a direct explanation of these lateral fluctuations and other phenomena for a subsurface with gentle deep structures and shallow overburden anomalies. The analytical conditions for this have been derived for a depth-velocity model with gentle structures with dips not exceeding 12°. The influence of lateral interval velocity changes and curvilinear overburden velocity boundaries can be estimated and analysed using these formulas. An analytical approach to normal moveout velocity analysis in a laterally inhomogeneous medium provides an understanding of the connection between lateral interval velocity changes and normal moveout velocities. In the presence of uncorrected shallow velocity anomalies, the difference between root-mean-square and stacking velocity can be arbitrarily large to the extent of reversing the normal moveout function around normal incidence traveltimes. The main reason for anomalous stacking velocity behaviour is non-linear lateral variations in the shallow overburden interval velocities or the velocity boundaries.
A special technique has been developed to determine and remove shallow velocity anomaly effects. This technique includes automatic continuous velocity picking, an inversion method for the determination of shallow velocity anomalies, improving the depth-velocity model by an optimization approach to traveltime inversion (layered reflection tomography) and shallow velocity anomaly replacement. Model and field data examples are used to illustrate this technique. 相似文献
A special technique has been developed to determine and remove shallow velocity anomaly effects. This technique includes automatic continuous velocity picking, an inversion method for the determination of shallow velocity anomalies, improving the depth-velocity model by an optimization approach to traveltime inversion (layered reflection tomography) and shallow velocity anomaly replacement. Model and field data examples are used to illustrate this technique. 相似文献
15.
永久散射体雷达干涉(PSI)技术及其应用于区域地表形变监测已成为雷达遥感领域的研究热点之一.使用单一卫星平台所获取的单侧视SAR影像时间序列进行PSI分析,仅能获取沿雷达视线(LOS)方向的一维地表位移信息.本文提出了基于多平台永久散射体雷达干涉提取三维地表形变速度场的模型与算法,其基本策略是:首先针对每一卫星平台的SAR影像时间序列进行PSI分析,并计算各地面目标沿LOS向的位移速度值,然后联合各平台所对应的LOS向位移速度值进行建模,并基于最小二乘方法解算各地面目标的三维位移速度分量.实验选取天津市西北部作为测试区,使用2007—2010年所获取的39幅TerraSAR-X影像、23幅ENVISAT ASAR影像和16幅ALOS PALSAR影像进行分析,经联合解算得到了该测试区域的垂直位移速度场以及南北向和东西向水平位移速度分量.与地面水准和已有GPS观测结果对比分析表明:多平台PSI的垂直位移速度场精度可达毫米级,而其水平位移速度分量与已有GPS结果基本一致.多平台PSI分析无需引入任何外部形变参考信息,便可以实现形变场的偏差校准和三维形变场的恢复. 相似文献
16.
Sylvie Tillard 《Geophysical Prospecting》1994,42(6):615-636
In order to understand various aspects of radar wave propagation, a survey of electromagnetic wave behaviour relative to the geological characteristics of the formations prospected was undertaken. The sites chosen for the tests were a granite quarry and an underground schist working. By investigating an electrically resistive isotropic site and a conductive anisotropic site, it was demonstrated that non-conventional use of a radar system (antennae raised, various orientations of the transmitter/receiver, etc.) could improve data quality, and could allow information other than reflector depth to be collected (volume scattering intensity, isotropy, etc.). By studying wave propagation velocities, we underlined the difficulties encountered in establishing a velocity versus depth law, despite recourse to seismic data processing, such as NMO corrections. The results of field experiments, complemented by laboratory measurements of dielectric permittivities, clearly showed anisotropy effects: in the case of a path that is perpendicular to the schistosity plane, an electromagnetic wave propagates more slowly and is more attenuated than a wave parallel to the schistosity plane. 相似文献
17.
Roberto de Franco 《Geophysical Prospecting》2001,49(4):395-404
A method to estimate interval velocities and thickness in a horizontal isotropic layered medium from wide-angle reflection traveltime curves is presented. The method is based on a relationship between the squared reflection traveltime differences and the squared offset differences relative to two adjacent reflectors. The envelope of the squared-time versus offset-difference curves, for rays with the same ray parameter, is a straight line, whose slope is the inverse of the square of the interval velocity and whose intercept is the square of the interval time. The method yields velocity and thickness estimates without any knowledge of the overlying stratification. It can be applied to wide-angle reflection data when either information on the upper crust and/or refraction control on the velocity is not available. Application to synthetic and real data shows that the method, used together with other methods, allows us to define a reliable 1D starting model for estimating a depth profile using either ray tracing or another technique. 相似文献
18.
Igor Ravve 《Geophysical Prospecting》2016,64(3):517-542
We suggest a new method to determine the piecewise‐continuous vertical distribution of instantaneous velocities within sediment layers, using different order time‐domain effective velocities on their top and bottom points. We demonstrate our method using a synthetic model that consists of different compacted sediment layers characterized by monotonously increasing velocity, combined with hard rock layers, such as salt or basalt, characterized by constant fast velocities, and low velocity layers, such as gas pockets. We first show that, by using only the root‐mean‐square velocities and the corresponding vertical travel times (computed from the original instantaneous velocity in depth) as input for a Dix‐type inversion, many different vertical distributions of the instantaneous velocities can be obtained (inverted). Some geological constraints, such as limiting the values of the inverted vertical velocity gradients, should be applied in order to obtain more geologically plausible velocity profiles. In order to limit the non‐uniqueness of the inverted velocities, additional information should be added. We have derived three different inversion solutions that yield the correct instantaneous velocity, avoiding any a priori geological constraints. The additional data at the interface points contain either the average velocities (or depths) or the fourth‐order average velocities, or both. Practically, average velocities can be obtained from nearby wells, whereas the fourth‐order average velocity can be estimated from the quartic moveout term during velocity analysis. Along with the three different types of input, we consider two types of vertical velocity models within each interval: distribution with a constant velocity gradient and an exponential asymptotically bounded velocity model, which is in particular important for modelling thick layers. It has been shown that, in the case of thin intervals, both models lead to similar results. The method allows us to establish the instantaneous velocities at the top and bottom interfaces, where the velocity profile inside the intervals is given by either the linear or the exponential asymptotically bounded velocity models. Since the velocity parameters of each interval are independently inverted, discontinuities of the instantaneous velocity at the interfaces occur naturally. The improved accuracy of the inverted instantaneous velocities is particularly important for accurate time‐to‐depth conversion. 相似文献
19.
This study uses a combination of evidence from ground penetrating radar, borehole, video, and wireless probe data to assess temporal changes in englacial water content associated with Briksdalsbreen, a rapidly retreating Norwegian glacier. Over a 13 day period in 2006, ice radar‐wave velocity varied between 0·135 m/ns (± 0·009) and 0·159 m/ns (± 0·003), and water content from 7·8% (+2·6, ?2·8) to 2·5% (+0·9, ?1·1) [derived from the Looyenga (Physica 31 (3): 401–406, 1965) formula]. It is suggested that during warm precipitation free days, void spaces within the glacier become filled with water, resulting in low radar‐wave velocity. This stored water then drained during cold, high precipitation days, allowing the radar‐wave velocity to rise. These changes in englacial storage were caused by the enhanced crevassing generated by the newly floating ice margin, and were associated with accelerated glacier retreat. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
20.
An inclusion model, based on the Kuster–Toksöz effective medium theory along with Gassmann theory, is tested to forward model velocities for fluid-saturated rocks. A simulated annealing algorithm, along with the inclusion model, effectively inverts measured compressional velocity (VP) to achieve an effective pore aspect ratio at each depth in a depth variant manner, continuously along with depth. Early Cretaceous syn-rift clastic sediments at two different depth intervals from two wells [well A (2160–2274 m) and well B (5222–5303 m)], in the Krishna–Godavari basin, India, are used for this study. Shear velocity (VS) estimated using modelled pore aspect ratio offers a high correlation coefficient (>0.95 for both the wells) with measured data. The modelled pore aspect ratio distribution suggests the decrease in pore aspect ratio for the deeper interval, mainly due to increased effective vertical stress. The pore aspect ratio analysis in relation to total porosity and volume of clay reveals that the clay volume has insignificant influence in shaping the pore geometry in the studied intervals. An approach based on multiple linear regression method effectively predicts velocity as a linear function of total porosity, the volume of clay and the modelled pore-space aspect ratio of the rock. We achieved a significant match between measured and predicted velocities. The correlation coefficients between measured and modelled velocities are considerably high (approximately 0.85 and 0.8, for VP and VS, respectively). This process indicates the possible influence of pore geometry along with total porosity and volume of clay on velocity. 相似文献