共查询到20条相似文献，搜索用时 15 毫秒
1.
《Journal of Applied Geophysics》2007,61(1):56-81
Due to the character of seismic energy generation and propagation, shallow high-resolution seismic-reflection surveys often fail in the identification of the shallowest horizons and, due to the limited offsets, accuracy of velocity analyses is often not very high.In recent years, Rayleigh wave dispersion analysis have proved to have good potential also for near-surface applications but dispersion curve inversion and related uncertainty evaluation pose serious problems to a completely stand-alone application.In order to overcome these problems a joint inversion scheme is proposed, which is based on the identification of the Pareto front, performed in the framework of a Multi-Objective Evolutionary Algorithm (MOEA). Seismic data considered to design the two objectives are the Rayleigh wave dispersion curve and reflection travel times.We initially analyse a set of synthetic cases and evaluate the obtained results. A significant improvement of the retrieved models is observed as long as reflection travel times are added to the dispersion curve alone.Furthermore, the proposed methodology also provides relevant indications about the consistency of the overall inversion process. In fact, the distribution of the models in the objective space, the trend of the objectives over the passing generations and the evolution of the Pareto front can provide useful information to evaluate the provisional tentative interpretation (number of strata and reflector identification) inherently adopted for the data inversion.On the basis of the results obtained from the tests on the synthetic datasets, the analyses of a field dataset are interpreted as possible evidence of lateral heterogeneities. 相似文献
2.
Joint inversion of Rayleigh‐wave dispersion data and vertical electric sounding data: synthetic tests on characteristic sub‐surface models 
下载免费PDF全文
下载免费PDF全文 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. 相似文献
3.
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. 相似文献
4.
Estimating gas saturation in a thin layer by using frequency‐dependent amplitude versus offset modelling 下载免费PDF全文
下载免费PDF全文 Zhaoyu Jin Mark Chapman Xiaoyang Wu Giorgos Papageorgiou 《Geophysical Prospecting》2017,65(3):747-765
Various models have been proposed to link partial gas saturation to seismic attenuation and dispersion, suggesting that the reflection coefficient should be frequency‐dependent in many cases of practical importance. Previous approaches to studying this phenomenon typically have been limited to single‐interface models. Here, we propose a modelling technique that allows us to incorporate frequency‐dependent reflectivity into convolutional modelling. With this modelling framework, seismic data can be synthesised from well logs of velocity, density, porosity, and water saturation. This forward modelling could act as a basis for inversion schemes aimed at recovering gas saturation variations with depth. We present a Bayesian inversion scheme for a simple thin‐layer case and a particular rock physics model and show that, although the method is very sensitive to prior information and constraints, both gas saturation and layer thickness theoretically can be estimated in the case of interfering reflections. 相似文献
5.
Hing-Ho Tsang M. Neaz SheikhNelson T.K. Lam 《Soil Dynamics and Earthquake Engineering》2012,34(1):89-98
Where a distinct soil-rock interface exists, the bedrock medium is commonly treated as elastic half-space and the bedrock surface as the lower boundary of the soil-column model for site response analyses (or the lower boundary of the finite element model for soil-structure interaction analyses). While shear wave velocity in bedrock varies with depth, there has been no consensus amongst scientists and practitioners over the value of “effective depth” into bedrock at which the “half-space” shear wave velocity value should be taken for modeling purposes. This paper reports an interesting and important observation that the effective depth into bedrock is sensitive to the shear wave velocity profile of the overlying soil sediments. A simple and heuristic method, namely Resonant Period Equivalence (RPE) Method, is proposed herein for representing a stratified elastic bedrock of inhomogeneous properties by an equivalent homogeneous elastic half-space medium, which is characterized by a single equivalent shear wave velocity (VR) value. The proposed calculation method has been verified by extensive comparative analyses involving the use of programs SHAKE and NERA and employing the complete shear wave velocity models of both the soil sediments and the underlying stratified bedrock. 相似文献
6.
Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm 总被引：1，自引：0，他引：1
Chong Zeng Jianghai Xia Richard D. MillerGeorgios P. Tsoflias 《Journal of Applied Geophysics》2011,75(4):648-655
Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity. 相似文献
7.
8.
While tilted transverse isotropy (TTI) is a good approximation of the velocity structure for many dipping and fractured strata, it is still challenging to estimate anisotropic depth models even when the tilted angle is known. With the assumption of weak anisotropy, we present a TTI traveltime inversion approach for models consisting of several thickness-varying layers where the anisotropic parameters are constant for each layer. For each model layer the inversion variables consist of the anisotropic parameters ε and δ, the tilted angle φ of its symmetry axis, layer velocity along the symmetry axis, and thickness variation of the layer. Using this method and synthetic data, we evaluate the effects of errors in some of the model parameters on the inverted values of the other parameters in crosswell and Vertical Seismic Profile (VSP) acquisition geometry. The analyses show that the errors in the layer symmetry axes sensitively affect the inverted values of other parameters, especially δ. However, the impact of errors in δ on the inversion of other parameters is much less than the impact on δ from the errors in other parameters. Hence, a practical strategy is first to invert for the most error-tolerant parameter layer velocity, then progressively invert for ε in crosswell geometry or δ in VSP geometry. 相似文献
9.
Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure 下载免费PDF全文
下载免费PDF全文 Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio (i.e., ellipticity), and receiver function data to better resolve 1D crustal shear wave velocity (v S) structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensitivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast (e.g., due to the existence of crustal interfaces) and v P/v S ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the v S model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute v S model and then incorporate receiver function data in the joint inversion to obtain a finer v S model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal v S structures and with little initial model dependency. 相似文献
10.
Ambient noise measurements have been performed at local and small scales in the Neapolitan and surrounding areas (Campania, southern Italy) by employing two broad-band Kinemetrics Q330 stations, equipped with Episensor ES-T three component accelerometers. In both experiments frequency time analysis (FTAN method) has been performed on the vertical and radial components of noise cross correlations to retrieve the Rayleigh wave dispersion (Green??s function). At local scale, over an interstation distance of about 26?km, the group velocity dispersion values have been compared with those obtained from FTAN analysis on recordings of two earthquakes with similar path. At small scale, measurements have been carried out over an interstation distance of about 440?m in the public gardens of Scampia, the northernmost quarter of Naples. The Rayleigh wave group velocity dispersion data obtained from noise cross correlation, have been combined with those from active seismic experiment along the same alignment, but shorter (120?m offset). The non linear inversion of such a dispersion curve has allowed the definition of V S models to depths of 100?m, in agreement with nearby stratigraphy. Moreover, a good agreement has resulted for the resonance frequency among the H/V ratio, the ellipticity of the fundamental mode computed for the chosen V S model, and the average two-dimensional (2D) spectral amplification computed along a cross section representative of the Scampia quarter. 相似文献
11.
面波多道分析方法(MASW)是获取垂向剪切波速度剖面的一种有效方法。频散曲线反演是MASW中关键的一步。由于瑞雷波频散曲线反演具有非线性、多参数和多极值的特征,这对于常规的局部线性化反演方法是极大的挑战。为此,本文采取确定性的全局优化算法,广义模式识别算法(GPS)对瑞雷波频散曲线进行反演。其原理可以简述为:算法首先通过模式以确定性的方式对目标函数进行采样来搜索一个点序列;然后使序列中每一个点到下一个点的目标函数值逐渐减少,从而使点序列逐渐逼近全局最优解,最后的解便为待求的最优模型参数。为验证GPS的有效性,首先利用设计的3种典型的6层地质模型通过快速矢量传递算法正演模拟产生基模式频散曲线(频率范围为5~101Hz,频率间隔为2Hz,频点数为49),并对理论频散曲线进行反演。反演结果表明,模型的真实值已经被高度精确地重建。说明GPS可以用于实际勘探中的基模式频散曲线反演。为进一步验证GPS的有效性,在吉林大学校园采集瑞雷波实测数据,并提取基模式频散曲线,应用GPS进行反演。反演重建的横波速度剖面与先验的地质信息吻合得很好。理论模型和真实数据的反演结果表明,GPS可以应用在瑞雷波频散曲线非线性反演中。 相似文献
12.
目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型... 相似文献
13.
地球物理反演浅地表土层波速剖面通常基于初始速度结构模型,为避免对地勘资料产生较多依赖,充分利用地震和背景噪声等被动源观测记录,快速简便地构建出土层反演的初始模型,本文基于场地土层的频散曲线和水平与竖向谱比提出了一种改进半波长法。该方法首先通过水平与竖向谱比确定场地土层的卓越频率,获取场地覆盖层厚度,确定反演所需的频带范围;其次,利用瑞雷波相速度对剪切波速的偏导数随土层深度的变化规律反演出场地的初始速度结构,并结合工程场地勘探中的半波长法,对常见的三类典型土层模型和日本Kushiro场地实测模型进行了实例分析;最后,将反演得到的速度结构与理论或实测速度结构进行误差对比分析。结果表明:改进半波长法获得的初始速度结构相对于理论或实测速度结构的最大误差不超过35%,可为利用地球物理方法反演工程场地波速剖面构造一个较小的搜索模型空间,进而提高反演计算的速度和结果的可靠性。 相似文献
14.
基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27. 相似文献
15.
The reliability of inversion of apparent resistivity pseudosection data to determine accurately the true resistivity distribution over 2D structures has been investigated, using a common inversion scheme based on a smoothness‐constrained non‐linear least‐squares optimization, for the Wenner array. This involved calculation of synthetic apparent resistivity pseudosection data, which were then inverted and the model estimated from the inversion was compared with the original 2D model. The models examined include (i) horizontal layering, (ii) a vertical fault, (iii) a low‐resistivity fill within a high‐resistivity basement, and (iv) an upfaulted basement block beneath a conductive overburden. Over vertical structures, the resistivity models obtained from inversion are usually much sharper than the measured data. However, the inverted resistivities can be smaller than the lowest, or greater than the highest, true model resistivity. The substantial reduction generally recorded in the data misfit during the least‐squares inversion of 2D apparent resistivity data is not always accompanied by any noticeable reduction in the model misfit. Conversely, the model misfit may, for all practical purposes, remain invariant for successive iterations. It can also increase with the iteration number, especially where the resistivity contrast at the bedrock interface exceeds a factor of about 10; in such instances, the optimum model estimated from inversion is attained at a very low iteration number. The largest model misfit is encountered in the zone adjacent to a contact where there is a large change in the resistivity contrast. It is concluded that smooth inversion can provide only an approximate guide to the true geometry and true formation resistivity. 相似文献
16.
在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想. 相似文献
17.
—An algorithm has been developed to compute the dispersive and dissipative seismic response using FUTTERMAN’S (1962) third attenuation-dispersion relationship. In the computation, frequency-dependent velocity and quality factor Q have been used but in the case of the nondispersive synthetic seismogram, frequency-independent velocity has been used. The model’s parameters are density, phase velocity, quality factors and thicknesses of the layers. Dispersive and nondispersive synthetic seismograms have been computed with and without absorption for a layered earth geological model. Fast Fourier transform (FFT) technique has been adopted for converting the frequency domain response into the time domain. The frequency spacing, Δf = 0.976?Hz, has been considered to avoid the aliasing effect. The results have revealed changes in the reflected waveforms in the frequency domain as well as in the time domain for absorption and dispersion cases. It is also concluded that dispersion reduces the arrival time and this effect is increasing with the travel time. The effect of constant Q on the seismic response has also been studied. 相似文献
18.
Ozlem Karagoz Kosuke Chimoto Hiroaki Yamanaka Oguz Ozel Seckin Citak 《Bulletin of Earthquake Engineering》2018,16(1):23-43
On 24 May 2014, a Mw 6.9 earthquake occurred in the west of Gokceada Island, northern Aegean Sea. The earthquake was close to Canakkale, Enez, Tekirdag cities, and damaged 300 buildings in the Marmara Region, NW Turkey. We simulated its broadband (0.1–10 Hz) ground motions including 1D deep and shallow structures soil amplification effects at the 12 strong ground motion stations in the western Marmara Region. The 1D deep velocity structures from the focal layer to the engineering bedrock with an S-wave velocity of 0.78 km/s in different azimuthal directions were tuned by comparing the observed group-velocity dispersion curves of Rayleigh and Love waves from the mainshock with theoretical ones. We also added the shallow parts from previous surveys into the 1D models. Synthetic seismograms on the engineering bedrock were generated using the discrete wave number method with a source model and the 1D deep velocity structures. Then the surface motion was generated considering shallow soil amplification. The synthetic seismograms are generally in good agreement with the observed low and high-frequency parts at most of the stations indicating an appropriateness of the source model and the 1D structural model. 相似文献
19.
This study proposes the use of multi-layer perceptron neural networks (MLPNN) to invert dispersion curves obtained via multi-channel analysis of surface waves (MASW) for shear S-wave velocity profile. The dispersion curve used in inversion includes the fundamental-mode dispersion data. In order to investigate the applicability and performance of the proposed MLPNN algorithm, test studies were performed using both synthetic and field examples. Gaussian random noise with a standard deviation of 4 and 8% was added to the noise-free test data to make the synthetic test more realistic. The model parameters, such as S-wave velocities and thicknesses of the synthetic layered-earth model, were obtained for different S/N ratios and noise-free data. The field survey was performed over the natural gas pipeline, located in the Germencik district of Ayd?n city, western Turkey. The results show that depth, velocity, and location of the embedded natural gas pipe are successfully estimated with reasonably good approximation. 相似文献
20.
瑞利面波垂直-水平振幅比(或ZH振幅比)是一个随频率变化的函数,对于台站下方浅层地壳结构非常敏感,且具有和频散资料不同的深度敏感核,是传统频散反演方法的一个很好的补充,从而可以将基阶瑞利面波的ZH振幅比和面波频散数据联合起来更好地反演获得观测台站下方的速度结构.本文提出了基于邻域算法的面波频散曲线与ZH振幅比联合反演方法,我们进行了基于理论模型的模拟测试,证明了联合反演是一种更为可靠的反演方法,且能更好地约束浅层地壳结构.相比于频散曲线单独反演,联合反演不仅可以精确反演获得地壳的Vs结构,对分层地壳的Vp/Vs也能很好地约束.然后我们将联合反演算法应用于实际测量数据,获得了中国西南昆明台(KMI)下方更为准确的地壳横波速度结构及Vp/Vs模型. 相似文献