首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 921 毫秒
1.
For the exploration of near-surface structures, seismic and geoelectric methods are often applied. Usually, these two types of method give, independently of each other, a sufficiently exact model of the geological structure. However, sometimes the inversion of the seismic or geoelectric data fails. These failures can be avoided by combining various methods in one joint inversion which feads to much better parameter estimations of the model than the independent inversions. A suitable seismic method for exploring near-surface structures is the use of dispersive surface waves: the dispersive characteristics of Rayleigh and Love surface waves depend strongly on the structural and petrophysical (seismic velocities) features of the near-surface Underground. Geoelectric exploration of the structure Underground may be carried out with the well-known methods of DC resistivity sounding, such as the Schlumberger, the radial-dipole and the two-electrode arrays. The joint inversion algorithm is tested by means of synthetic data. It is demonstrated that the geoelectric joint inversion of Schlumberger, radial-dipole and two-electrode sounding data yields more reliable results than the single inversion of a single set of these data. The same holds for the seismic joint inversion of Love and Rayleigh group slowness data. The best inversion result is achieved by performing a joint inversion of both geoelectric and surface-wave data. The effect of noise on the accuracy of the solution for both Gaussian and non-Gaussian (sparsely distributed large) errors is analysed. After a comparison between least-square (LSQ) and least absolute deviation (LAD) inversion results, the LAD joint inversion is found to be an accurate and robust method.  相似文献   

2.
The applicability of three kinds of electrode configurations used to delineate a buried horizontal pipe was studied. A 3D resistivity imaging survey was carried out along eight parallel lines using pole-pole, pole-dipole, and dipole-dipole arrays with 1m minimum electrode spacings. Roll-along measurements were carried out to cover a rectangular grid. The 2D and 3D least squares algorithms based on the robust inversion method were used in the inversion of the apparent resistivity data sets. The 2D inversion of data sets could not delineate the orientation and dimension of the subsurface anomalies clearly. To obtain more accurate results, a 3D joint inversion of the pole-pole and pole-dipole data sets was performed, as well as of pole-pole and dipole-dipole data sets. In this case, both horizontal and vertical dimensions of subsurface structures were resolved. The resulting model obtained from each array was compared to those of joint inversion method. The result showed that the horizontal resolution does not improve so much as that in the vertical direction when joint inversion is applied.  相似文献   

3.
This study applies array methods to measure the relative proportions of Love and Rayleigh waves in the ambient vibration wavefield. Information on these properties is of special relevance for frequencies around the horizontal-to-vertical (H/V) spectral amplitude ratio peak. The analysis of H/V curves, a popular technique in site characterisation, commonly assumes that the curves represent the frequency-dependent Rayleigh wave ellipticity. For the detailed interpretation of amplitudes or the inversion of the curves, it is therefore necessary to estimate and correct for the contribution of other wave types to the ambient vibration wavefield. I use available ambient vibration array measurements to determine the relative amount of Love and Rayleigh waves on the horizontal components by frequency-dependent analysis of the main propagation and polarisation directions, with a special emphasis on the H/V peak frequency as determined from the same recordings. Tests with synthetic data demonstrate the feasibility of this approach, at least in the presence of dominant source regions. Analysis of the data from 12 measurements at nine European sites, which include shallow as well as deep locations that span a wide range of impedance contrasts at the sediment-bedrock interface, indicates that the relative contribution of Rayleigh waves varies widely with frequency, from close to 0% to more than 70%. While most data sets show relative Rayleigh wave contributions between 40% and 50% around the H/V peak, there are also examples where Love waves clearly dominate the wavefield at the H/V peak, even for a site with a low impedance contrast. Longer-term measurements at one site indicate temporal variations in the relative Rayleigh wave content between day- and nighttime. Results calculated with the method introduced herein generally compare well with results of modified spatial autocorrelation analysis. These two methods might be used in a complimentary fashion, as both rely on different properties of the ambient vibration wavefield. This study illustrates that it is possible to measure the relative Rayleigh wave content of the noise wavefield from array data. Furthermore, the examples presented herein indicate it is important to estimate this property, as the assumption that there are an equal proportion of Love and Rayleigh waves is not always correct.  相似文献   

4.
The objective of this paper is to investigate the applicability of four-electrode arrays in 3D electrical resistivity imaging survey. A 3D resistivity imaging survey was carried out along fourteen parallel lines using dipole-dipole, Wenner-Schlumberger, and Wenner arrays with 2 m minimum electrode spacings. Roll-along measurements using a line spacing of 1 m were carried out covering a grid of 20 × 14 electrodes. The 3D least squares algorithm, based on the robust inversion method, was used in the inversion of the 3D apparent resistivity data sets. The results show that the 3D electrical resistivity imaging survey using the Wenner-Schlumberger and the dipole-dipole arrays, or the Wenner and the dipole-dipole arrays, in combination with an appropriate 3D inversion method, can be highly useful when the site conditions do not allow using the pole-pole or pole-dipole arrays.  相似文献   

5.
Alfano L. 1980, Dipole-Dipole Deep Geodectric Soundings over Geological Structures, Geophysical Prospecting 28, 283–296. The dipole-dipole array is used for deep vertical geoelectric soundings for practical reasons. However, the data obtained are often too scattered and the quality is lower than that of a Schlumberger sounding. The mathematical transform from these fields data to the corresponding Schlumberger ones is needed. It is demonstrated that for more general underground structures different from plane parallel stratifications, only the continuous polar dipole-dipole array permits a general and reliable transform of this type. It follows that any other dipole-dipole array, (for instance the azimuthal, equatorial, parallel, perpendicular, etc.), are useless over complex structures, or in the presence of irregularities on the ground surface.  相似文献   

6.
In the traditional inversion of the Rayleigh dispersion curve, layer thickness, which is the second most sensitive parameter of modelling the Rayleigh dispersion curve, is usually assumed as correct and is used as fixed a priori information. Because the knowledge of the layer thickness is typically not precise, the use of such a priori information may result in the traditional Rayleigh dispersion curve inversions getting trapped in some local minima and may show results that are far from the real solution. In this study, we try to avoid this issue by using a joint inversion of the Rayleigh dispersion curve data with vertical electric sounding data, where we use the common‐layer thickness to couple the two methods. The key idea of the proposed joint inversion scheme is to combine methods in one joint Jacobian matrix and to invert for layer S‐wave velocity, resistivity, and layer thickness as an additional parameter, in contrast with a traditional Rayleigh dispersion curve inversion. The proposed joint inversion approach is tested with noise‐free and Gaussian noise data on six characteristic, synthetic sub‐surface models: a model with a typical dispersion; a low‐velocity, half‐space model; a model with particularly stiff and soft layers, respectively; and a model reproduced from the stiff and soft layers for different layer‐resistivity propagation. In the joint inversion process, the non‐linear damped least squares method is used together with the singular value decomposition approach to find a proper damping value for each iteration. The proposed joint inversion scheme tests many damping values, and it chooses the one that best approximates the observed data in the current iteration. The quality of the joint inversion is checked with the relative distance measure. In addition, a sensitivity analysis is performed for the typical dispersive sub‐surface model to illustrate the benefits of the proposed joint scheme. The results of synthetic models revealed that the combination of the Rayleigh dispersion curve and vertical electric sounding methods in a joint scheme allows to provide reliable sub‐surface models even in complex and challenging situations and without using any a priori information.  相似文献   

7.
Horizontally layered (1D) earth models are often assumed as a model estimate for the interpretation of geophysical data measured along 2D geological structures. In this process, the individual data sets are usually inverted independently, and it is considered only in a later phase of interpretation that these local (1D) models have common characteristic features. Taking account of these common attributes, instead of the successive independent interpretations, the lateral variations of geometrical and petrophysical parameters can be efficiently determined for the whole 2D structure by applying a series expansion. Using global basis functions, two advantages can be achieved: (i) choosing an appropriate number of basis functions helps us to restrict the complexity of the model; (ii) the integration of all the data sets measured along the profile gives rise to the application of simultaneous or joint inversion methods. This results in a decrease of the number of independent unknowns, a higher stability during the inversion and a more accurate and reliable parameter estimation.In this paper, a joint inversion algorithm is presented using DC geoelectric apparent resistivities and refraction seismic travel times measured along various layouts above a 2D geological model. To describe lateral variations series, expansions are used, and furthermore, to improve the often used approximation of a (locally) 1D forward modelling, the integral mean value of the horizontally changing model parameters (calculated along an appropriately defined interval) is introduced. We call the inversion procedure that combines series expansions and the concept of integral mean Generalised Series Expansion (GSE) inversion. The method was developed and tested for both the simultaneous (integrating data sets of one method or methods on the same physical basis) and the joint inversion (where data sets of methods on different physical bases are joined together), using synthetic and field data sets. It is also demonstrated that the equivalence problem inherent in the independent inversion of DC geoelectric data can efficiently be resolved by the use of the joint GSE inversion method in the cases of conductive and resistive equivalent geological models.  相似文献   

8.
The methods of in-mine seam-sounding and transillumination (geoelectric tomography) for the detection of tectonic disturbances of coal seams were developed by the Department of Geophysics of the University of Miskolc in the 1970–80’s with the effective support of the former “Borsod” Coal Mines Ltd. The paper gives an overview about the theory of seam-sounding and a special geoelectric tomographic inversion, and introduces the in-mine geoelectric seam-sounding and transillumination measurement systems using vertical electrode dipoles. In the second part the paper, the results of an in-mine geoelectric measurement are presented, which was carried out in order to detect tectonic disturbances of the Miocene aged coal seams situated in Slovakia. As results of the geophysical investigation, the authors forecasted the tectonic features in the coal seam. The company confirmed the results by independent information about seam disturbances and tectonic features arising from the excavation of the investigated area.  相似文献   

9.
The local geology and shallow S-wave velocity structure of a site are recognized to be key factors for the increase in the damaging potential of seismic waves. Indeed, seismic amplitudes may be amplified in frequency ranges unfavorable for building stock by the presence of soft sedimentary covers over lying hard bedrock. Hence, microzonation activities, which aim at assessing the site response as accurately as possible, have become a fundamental task for the seismic risk reduction of urbanized areas. Methods based on the measurement of seismic noise, which typically are fast, non-invasive, and low cost, have become a very attractive option in microzonation studies.Using observations derived from seismic noise recordings collected by two-dimensional arrays of seismic stations, we present a novel joint inversion scheme for surface wave curves. In particular, the Love wave, the Rayleigh wave dispersion and the HVSR curves are innovatively combined in a joint inversion procedure carried out following a global search approach (i.e., the Genetic Algorithm).The procedure is tested using a data set of seismic noise recordings collected at the Bevagna (Italy) test-site. The results of the novel inversion scheme are compared with the inversion scheme proposed by Parolai et al. (2005), where only Rayleigh wave dispersion and HVSR curves are used, and with a cross-hole survey.  相似文献   

10.
As theory dictates, for a series of horizontal layers, a pure, plane, horizontally polarized shear (SH) wave refracts and reflects only SH waves and does not undergo wave-type conversion as do incident P or Sv waves. This is one reason the shallow SH-wave refraction method is popular. SH-wave refraction method usually works well defining near-surface shear-wave velocities. Only first arrival information is used in the SH-wave refraction method. Most SH-wave data contain a strong component of Love-wave energy. Love waves are surface waves that are formed from the constructive interference of multiple reflections of SH waves in the shallow subsurface. Unlike Rayleigh waves, the dispersive nature of Love waves is independent of P-wave velocity. Love-wave phase velocities of a layered earth model are a function of frequency and three groups of earth properties: SH-wave velocity, density, and thickness of layers. In theory, a fewer parameters make the inversion of Love waves more stable and reduce the degree of nonuniqueness. Approximating SH-wave velocity using Love-wave inversion for near-surface applications may become more appealing than Rayleigh-wave inversion because it possesses the following three advantages. (1) Numerical modeling results suggest the independence of P-wave velocity makes Love-wave dispersion curves simpler than Rayleigh waves. A complication of “Mode kissing” is an undesired and frequently occurring phenomenon in Rayleigh-wave analysis that causes mode misidentification. This phenomenon is less common in dispersion images of Love-wave energy. (2) Real-world examples demonstrated that dispersion images of Love-wave energy have a higher signal-to-noise ratio and more focus than those generated from Rayleigh waves. This advantage is related to the long geophone spreads commonly used for SH-wave refraction surveys, images of Love-wave energy from longer offsets are much cleaner and sharper than for closer offsets, which makes picking phase velocities of Love waves easier and more accurate. (3) Real-world examples demonstrated that inversion of Love-wave dispersion curves is less dependent on initial models and more stable than Rayleigh waves. This is due to Love-wave’s independence of P-wave velocity, which results in fewer unknowns in the MALW method compared to inversion methods of Rayleigh waves. This characteristic not only makes Love-wave dispersion curves simpler but also reduces the degree of nonuniqueness leading to more stable inversion of Love-wave dispersion curves.  相似文献   

11.
Love and Rayleigh wave phase velocities are analyzed with the goal of retrieving information about the anisotropic structure of the Iberian lithosphere. The cross-correlation method is used to measure the interstation phase velocities between diverse stations of the ILIHA network at periods between 20 and 120 s. Despite the 2-D structure of the network, the Love wave data are too few to enable an analysis of phase velocity azimuthal variations. Azimuthal averages of Love and Rayleigh wave phase velocities are calculated and inverted both in terms of isotropic and anisotropic structures. Realistic isotropic models explain the Rayleigh wave and short-period Love wave phase velocities. Therefore no significant anisotropy needs to be introduced in the crust and down to 100 km depth in the upper mantle to explain our data. A discrepancy is observed only at long periods, where the data are less reliable. Love wave data at periods between 80 and 120 s remain 0.15 km/s faster than predicted by isotropic models explaining the long-period Rayleigh wave data. Possibilities of biases in the measurements due to interferences with higher modes are examined but seem unlikely. A transversely isotropic model with 8% of S-wave velocity anisotropy in the upper mantle at depths larger than 100 km can explain the whole set of data. In terms of a classical model of mantle anisotropy, this corresponds to 100% of the crystals perfectly oriented in the horizontal plane in a pyrolitic mantle. This is a rather extreme model, which predicts at time delay between 0 and 2 seconds for split SKS.  相似文献   

12.
利用高密度电阻率法确定滑坡面的正演研究   总被引：2，自引：0，他引：2  
高密度电法探测是滑坡面探测的有效手段,适宜 的现场工作方案和相应的资料解译模型,是获取 可靠探测成果的前提。本文针对高密度电法勘探 中常用的四类电极排列方式,即维纳、施伦贝谢、 偶极、二极,利用有限元的方法二维正演模拟建 立了均质土坡滑面、软弱夹层滑面、堆积层滑面、 构造破碎带滑面的标准地电模型剖面;明确了进 行不同类型滑坡工程地质调查时不同电测装置的 探测能力和视电阻率异常模式;提出复杂地质条 件下斜坡问题探测数据有效处理的设计比值参数 方法。最后给出若干利用高密度电阻率法进行滑 坡工程地质调查的典型实例。  相似文献   

13.
傅淑芳  程宁亚 《地震学报》1988,10(4):352-362
本文提出了一个利用平面上长方形区域內的面波频散资料,求区域的三维速度结构的方法。将地震波慢度表示成二重Fourier级数,反演其系数为深度的函数,最后合成速度值,并以Love面波为例作了数值计算试验,绘制出了深度为150km的速度平面分布图。   相似文献   

14.
An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results.  相似文献   

15.
高频面波方法的若干新进展   总被引：12，自引：5，他引：7       下载免费PDF全文
面波多道分析方法(MASW)通过分析高频瑞雷波确定浅地表剪切波速度.在过去的20年中,由于该方法具有非侵入性、无损、高效及价格低的特点,越来越受到浅地表地球物理和地质工程学界的重视,视为未来最有希望的技术之一.这篇综述论文将介绍中国地质大学(武汉)浅地表地球物理团队近年来在研究高频面波的传播理论和应用中取得的部分成果.非几何波是一种仅存在于浅地表介质,尤其是未固结的沉积物中的独特的地震波.它的存在对快速而准确地获得表层S波速度有一定价值.我们的研究表明非几何波是一种具有频散特性的泄漏波.泄漏波的存在可能导致将其误认为瑞雷波的基阶或高阶能量,从而造成模式误判.这种模式误判会导致错误的反演结果.我们通过求取高基阶分离后的瑞雷波格林函数证明虚震源法瑞雷波勘探的可行性.这个结果将极大地降低野外瑞雷波勘探成本.勒夫波多道分析方法(MALW)中未知参数比瑞雷波的少,这使得勒夫波的频散曲线比瑞雷波的简单.因此,勒夫波反演更稳定,非唯一性更低.勒夫波数据生成的能量图像通常比瑞雷波的清晰,并具有更高的分辨率,从而可以更容易地拾取精确的勒夫波的相速度.利用雅克比矩阵分析波长与探测深度的关系表明对相同波长的基阶模式而言,瑞雷波的探测深度是勒夫波的1.3~1.4倍;而两种波的相同波长的高阶模式波的探测深度相同.我们也尝试了时间域勒夫波反演.按照勒夫波分辨率将地球模型剖分成了不同尺寸的块体,利用反卷积消除了地震子波对勒夫波波形的影响,通过更新每个块体的S波速度来拟合勒夫波波形,从而获得地下S波速度模型.该方法不基于水平层状模型假设,适用于任意二维介质模型.  相似文献   

16.
Joint analysis of Rayleigh- and Love-wave dispersion is performed with the aim of evaluating how their joint use can improve retrieved vertical VS profiles. In fact, non-uniqueness of the solution and complex energy distribution among different modes represent problems which, if not properly considered, can eventually lead to ambiguous or erroneous subsurface models.Some tests performed on synthetic datasets show that for the deepest layers the improvements obtained by the joint inversion cannot be considered as fully decisive in terms of ultimate solution of non-uniqueness. Nevertheless joint analysis of dispersive properties of Rayleigh and Love waves reveals as a highly valuable tool able to clarify possible interpretation issues of the single components. Under some stratigraphical circumstances, velocity spectra of Rayleigh waves can in fact be extremely complex in terms of energy distribution among different modes and erroneous interpretations of dispersion curves can thus occur. Beneficial aspects of the joint analysis is shown in the light of possible inconsistencies of the Pareto front, since major interpretative errors can be revealed in the outcomes of the proposed inversion procedure. Two field datasets are analysed also suggesting some improvements in the field acquisition procedures aimed at the acquisition of both Rayleigh and Love waves.  相似文献   

17.
Refracted arrivals are analysed to estimate the near‐surface anisotropy of marine sediments using a vertical‐cable (VC) configuration. In the presence of dip, the horizontal and vertical ray‐slownesses are obtained from the observed apparent slownesses in the up‐ and downdip directions using a sum or difference at each azimuth. The multiple azimuths generated by a VC geometry permit the ray‐slowness distribution of the marine sediments to be determined. An inversion procedure is developed to provide dip and anisotropy parameters for refractive layers from the measured refraction traveltimes in multilayered azimuthally isotropic and anisotropic media. Two sets of transversely isotropic models are used to analyse the azimuthal variations of apparent and ray slownesses. In the first set, we fix the anisotropic parameters of the models but vary the dip (0°, 5° and 10°) to test the effects of the presence of dip. In the second set, we vary the P‐wave anisotropy strength (5.2%, 10.3%, 15.8% and 22.0%) to examine the sensitivity and accuracy of ray‐slowness approximations which are independent of dip. We test this inversion procedure on synthetic P‐wave VC data calculated for six different models by a finite‐difference method. The results of applications to real VC data acquired from the North Sea are also presented.  相似文献   

18.
冀战波  王宝善 《地球物理学报》1954,63(11):4097-4113
2015年8月12日天津滨海新区发生的强烈化学品爆炸造成了巨大的经济损失和社会影响.天津爆炸产生了清晰的大振幅面波信号,分析结果表明这组信号由基阶和高阶面波组成,可以追踪到约135 km外的远处台站.利用这组面波信号分别开展了以下研究：（1）利用附近三个台站记录的四个单频基阶Rayleigh波信号对爆破事件的绝对位置进行了网格搜索,结果与利用GPS测量的位置相差仅0.498 km;（2）分别利用网格搜索和主事件定位法,对两次子事件的相对位置进行了确定,距离约75 m左右,与前人研究结果吻合;（3）从面波记录中测量到36条基阶Rayleigh波、49条第一高阶Rayleigh波、9条基阶Love波和29条第一高阶Love波的频散曲线,并进一步反演获得研究区域地下4 km内的S波速度结构.反演结果显示地表处S波速度低至0.375 km·s-1,在小于1 km的浅地表速度梯度较大,符合典型的盆地结构特征.本文的研究结果为类似爆炸等突发事件快速定位提供了新的思路,有助于灾后救援的迅速展开;同时得到天津滨海新区及周边浅层精细的速度结构,对于地震灾害评估有很大帮助.  相似文献   

19.
By exploiting the capability of identifying and extracting surface waves existing in a seismic signal, we can proceed to estimate the angular displacement (rotation about the horizontal axis normal to the direction of propagation of the wave; rocking) associated with Rayleigh waves as well as the angular displacement (rotation about the vertical axis; torsion) associated with Love waves.For a harmonic Rayleigh (Love) wave, rocking (torsion) would be proportional to the harmonic vertical (transverse horizontal) velocity component and inversely proportional to the phase velocity corresponding to the particular frequency of the harmonic wave (a fact that was originally exploited by Newmark (1969) [15] to estimate torsional excitation). Evidently, a reliable estimate of the phase velocity (as a function of frequency) is necessary. As pointed out by Stockwell (2007) [17], because of its absolutely referenced phase information, the S-Transform can be employed in a cross-spectrum analysis in a local manner. Following this suggestion a very reliable estimate of the phase velocity may be obtained from the recordings at two nearby stations, after the dispersed waves have been identified and extracted. Synthesis of the abovementioned harmonic components can provide a reliable estimate of the rocking (torsional) motion induced by an (extracted) Rayleigh (Love) wave.We apply the proposed angular displacement estimation procedure for two well recorded data sets: (1) the strong motion data generated by an aftershock of the 1999 Chi-Chi, Taiwan earthquake and recorded over the Western Coastal Plain (WCP) of Taiwan, and (2) the strong motion data generated by the 2010 Darfield, New Zealand earthquake and recorded over the Canterbury basin. The former data set is dominated by basin-induced Rayleigh waves while the latter contains primarily Love waves.  相似文献   

20.
Data collected in a dipole-dipole sounding along a crooked line can be transformed to form an approximately equivalent Schlumberger sounding, using a simple matrix inversion technique. The equivalent curve can be interpreted using rapid interpretation methods.  相似文献   

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

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