共查询到20条相似文献,搜索用时 156 毫秒
1.
由于瑞雷波存在多个导波模式,仅用基阶模式波频散数据反演而忽略高模式波,会影响反演精度.本文利用高阶模式瑞雷波频散数据与基阶模式频散数据互相独立的特性,构建瑞雷波基阶波、一阶高模式波和二阶高模式波的联合反演目标函数.通过增加反演变量,利用权系数调节各变量对反演结果的影响,实现同步联合反演,从而保证反演稳定性.使用人工蜂群智能算法求解目标函数,有效解决了传统反演方法将非线性问题线性化处理产生的细节信息丢失问题,提高了反演精度.经过对较为复杂的四层含低速夹层和四层含高速夹层理论模型和实例数据的反演验证,表明基于人工蜂群算法的瑞雷波多阶模式波非线性联合反演能够提高反演的稳定性和精度. 相似文献
2.
反演高频率(≥2Hz)的瑞雷波的相速度可得到分层地球模型从地表到地下30m以内的剪切(S)波速度。如果已知S波速度VS、压缩(P)波速度VP和瑞雷波速度,通过反演瑞雷波衰减系数就可以获得分层地球模型的P波品质因子QP和S波品质因子QS。模拟结果证明,根据瑞雷波衰减系数反演品质因子QS是可行的。当VS/VP达到0.45时,不能忽略QP对瑞雷波衰减系数的贡献,这种情况在近地表构造中并非罕见。从某些地质构造中的瑞雷波衰减系数反演得到QP是可能的,这是一种不同于一般认为的观点,即相对于P波的品质因子QP,瑞雷波的衰减系数对S波的品质因子QS更为敏感。在亚利桑那沙漠,采集了60道的面波数据。对一个层厚度超过20m的10层模型,首先利用多道面波分析(MASW)方法反演数据得到S波速度,然后通过反演衰减系数确定品质因子。 相似文献
3.
为充分利用微动信号中基阶和高阶模式瑞雷波,本文研究了基于多阶瑞雷波SPAC系数直接反演的方法.该方法首先基于地层介质响应计算多阶瑞雷波的能量占比,考虑实际观测台阵有限台站个数对SPAC系数影响,正演计算多阶瑞雷波SPAC系数,再采用快速模拟退火算法对其反演以获得地下介质横波速度结构.在此基础上,本文通过数值模拟验证该方法的可靠性,分别选取三种典型地质模型,基于模式叠加算法合成理论微动信号,采用本文方法计算其理论多阶瑞雷波SPAC系数并反演,给出反演结果与真实模型对比.我们将该方法应用于上海中心城区的地质调查中,通过与钻探结果对比,进一步验证该方法的有效性.本文理论与实际应用研究表明,基于多阶瑞雷波SPAC系数直接反演的微动探测方法有助于提高反演结果的可靠性,尤其对含软硬夹层的复杂地层介质,可提高探测精度. 相似文献
4.
研究了六层层状介质模型瑞利波基阶和高模式波相速度对横波速度、深度的敏感性,结果表明:基阶波较高模式波对7 m以内浅部地层的横波速度更敏感,敏感性频带在10~25 Hz范围内,峰值频带集中在18 Hz左右;高模式波较基阶波对深部地层的横波速度更敏感,敏感性频带宽,峰值分散.基阶波对浅层的敏感性和高模式波穿透深度更深的特点为近地表岩土层二维横波速度结构的联合反演提供了前提条件.利用阻尼最小二乘SVD(Singular Value Decomposition)算法联合基阶与高模式波对理论模型和实例数据进行横波速度反演,反演结果表明联合反演增强了反演的稳定性,提高了反演的精度. 相似文献
5.
为充分利用微动信号中基阶和高阶模式瑞雷波,本文研究了基于多阶瑞雷波SPAC系数直接反演的方法.该方法首先基于地层介质响应计算多阶瑞雷波的能量占比,考虑实际观测台阵有限台站个数对SPAC系数影响,正演计算多阶瑞雷波SPAC系数,再采用快速模拟退火算法对其反演以获得地下介质横波速度结构.在此基础上,本文通过数值模拟验证该方法的可靠性,分别选取三种典型地质模型,基于模式叠加算法合成理论微动信号,采用本文方法计算其理论多阶瑞雷波SPAC系数并反演,给出反演结果与真实模型对比.我们将该方法应用于上海中心城区的地质调查中,通过与钻探结果对比,进一步验证该方法的有效性.本文理论与实际应用研究表明,基于多阶瑞雷波SPAC系数直接反演的微动探测方法有助于提高反演结果的可靠性,尤其对含软硬夹层的复杂地层介质,可提高探测精度.
相似文献6.
《地球物理学进展》2016,(6)
高频面波方法是以瑞雷波和勒夫波为研究对象,在水文、工程和环境研究中具有广泛的应用,其主要思想是利用面波的频散特性来反映浅层地质问题,提取高质量的频散曲线成了面波勘探中的一个重要步骤.本文通过对理论模型研究说明,利用τ-p变换与F-K变换两种方法均能有效地进行频散能量成像,其效果在高频范围差别不大,在低频范围F-K变换提取频散曲线更为精确;与瑞雷波记录相比,勒夫波记录信噪比较高,通过对同一地区采集的瑞雷波和勒夫波进行频散分析可知,瑞雷波记录频散能量成像效果差,频散能量出现间断、跳跃的现象,而勒夫波记录频散能量连续、集中,相速度分辨率较高,提取频散曲线简单清晰,可为反演地表横波精细速度结构提供更好的频散曲线数据. 相似文献
7.
将一种新的方法——频率-贝塞尔变换法(F-Jmethod)应用于日本NIED在关东盆地布设的MeSO-net台网的背景噪声数据中,证明频率-贝塞尔变换法可以有效地从背景噪声中提取高阶频散曲线.利用提取的基阶和高阶频散曲线反演关东盆地区域的沉积层和地震基岩层的S波速度结构,并将我们反演得到的S波速度结构与Koketsu等提出的日本综合速度结构模型进行对比讨论.我们的例子证明,在基阶面波的基础上,高阶面波能减少在反演中的非唯一性,得到更为准确的S波速度结构.
相似文献8.
9.
微动探测方法是在背景噪声成像理论的基础上发展起来的一种新型地球物理探测手段,其根据获取野外数据的观测台站方式分为微动阵列方法和微动单台方法.本文综述了前人对微动探测方法的研究进展,并对目前该方法存在的问题及研究的发展趋势进行了展望.微动阵列方法采用频率波数法或空间自相关法从微动信号中提取频散曲线,进而反演获得地下横波速度结构.微动单台方法即H/V谱比法以三分量台站为观测基础,对水平分量信号与垂直分量信号进行傅里叶谱比,以此估算沉积层厚度、获得共振频率以及场地放大因子等.目前微动探测发展还有诸多问题需要进行深入研究,例如瑞雷波与勒夫波浅层探测理论方法、同时考虑基阶和高阶面波的影响进行多阶面波反演方法研究以及H/V谱比曲线与频散曲线联合反演方法研究等,从而提高微动探测方法浅层探测精度以满足城市地下空间精细探测的需求,加深探测深度为透明地壳、地下第三空间的精细探测提供一种新型地球物理手段. 相似文献
10.
采用与作者2014年发表的“大别-苏鲁及其邻近地区基于背景噪声的勒夫波群速度成像”文章相同的资料,用频时分析提取5 000余条瑞雷波和4 000余条勒夫波相速度频散曲线,反演得到了8—32 s的瑞雷波和勒夫波相速度分布图像.结果显示,瑞雷波与勒夫波相速度分布具有很好的一致性.8 s的相速度分布与地表构造特征相吻合,造山带与隆起区均表现为高速,盆地因其规模不同而显示不同程度的低速.随着周期的增大,大别 苏鲁的高速带由强变弱,但始终存在.16—24 s的高速可能主要受到中地壳高速的控制,而32 s的高速则可能与上地幔顶部的高速有关.比较大别造山带与苏鲁造山带的平均频散曲线,发现大别造山带和苏鲁造山带的勒夫波频散曲线均高于AK135模型计算的理论频散曲线,而瑞雷波则没有这一现象. 这可能意味着两个地区有比较强烈的径向各向异性. 相似文献
11.
Advantages of Using Multichannel Analysis of Love Waves (MALW) to Estimate Near-Surface Shear-Wave Velocity 总被引:2,自引:1,他引:1
Jianghai Xia Yixian Xu Yinhe Luo Richard D. Miller Recep Cakir Chong Zeng 《Surveys in Geophysics》2012,33(5):841-860
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. 相似文献
12.
Rayleigh-wave phase velocities have been utilized to determine shear (S)-wave velocities in near-surface geophysics since early 1980s. One of the key steps is to calculate theoretical dispersion curves of an earth model. When the S-wave velocity of the surface layer is higher than some of the layers below, however, the Rayleigh-wave phase velocity in a high-frequency range calculated by existing algorithms approaches the lowest S-wave velocity among the layers above the half-space, rather than a value related to the S-wave velocity of the surface layer. According to our numerical modeling results based on wave equation, trends of the Rayleigh-wave dispersive energy approach about a 91% of the S-wave velocity of the surface layer at a high-frequency range when its wavelength is much shorter than the thickness of the surface layer, which cannot be fitted by a dispersion curve calculated by existing algorithms. We propose a method to calculate Rayleigh-wave phase velocities of models with a high-velocity surface layer by considering its penetration depth. We build a substituted model that only contains the layer with the lowest S-wave velocity among the layers above the half-space and the layers above it. We use the substituted model to replace the original model to calculate phase velocities when the Rayleigh-wave wavelength is not long enough to penetrate the lowest S-wave velocity layer. Several synthetic models are used to verify fitness between the dispersion curve calculated by our proposed method and the trend of the highest dispersive energy. Examples of inversion also demonstrate high accuracy of using our method as the forward calculation method during the inversions. 相似文献
13.
High-frequency (≥ 2 Hz) Multi-channel Analysis of Love Waves (MALW) provides a practical way to determine velocity of horizontally polarized shear (SH) waves for a layered earth model up to 30 m below the ground surface in many geological settings. The information used in the MALW method is phase of Love waves. Information on amplitude of Love waves is not utilized in the MALW method. In this paper we present a method that uses information on amplitude of high-frequency Love waves to estimate quality factors (Qs) of near-surface materials. Unlike Rayleigh waves, attenuation coefficients (amplitude) of Love waves are independent of quality factors for P waves and are function of quality factors of Love waves. In theory, a fewer parameters make the inversion of attenuation coefficients of Love waves more stable and reduce the degree of nonuniqueness. We discussed sensitivity of an inversion system based on a linear relationship between attenuation coefficients and dissipation factors (1/Qs). The sensitivity analysis suggested that damping and constraints to an inversion system are necessary to obtain a smooth and meaningful quality factor model when no other information is available. We used synthetic and real-world data to demonstrate feasibility of inversion of attenuation coefficients of high-frequency Love-wave data acquired with the MALW method for quality factors with a linear, damped and constrained system. 相似文献
14.
Finite-Difference Modeling and Dispersion Analysis of High-Frequency Love Waves for Near-Surface Applications 总被引:3,自引:0,他引:3
Yinhe Luo Jianghai Xia Yixian Xu Chong Zeng Jiangping Liu 《Pure and Applied Geophysics》2010,167(12):1525-1536
Love-wave propagation has been a topic of interest to crustal, earthquake, and engineering seismologists for many years because it is independent of Poisson’s ratio and more sensitive to shear (S)-wave velocity changes and layer thickness changes than are Rayleigh waves. It is well known that Love-wave generation requires the existence of a low S-wave velocity layer in a multilayered earth model. In order to study numerically the propagation of Love waves in a layered earth model and dispersion characteristics for near-surface applications, we simulate high-frequency (>5 Hz) Love waves by the staggered-grid finite-difference (FD) method. The air–earth boundary (the shear stress above the free surface) is treated using the stress-imaging technique. We use a two-layer model to demonstrate the accuracy of the staggered-grid modeling scheme. We also simulate four-layer models including a low-velocity layer (LVL) or a high-velocity layer (HVL) to analyze dispersive energy characteristics for near-surface applications. Results demonstrate that: (1) the staggered-grid FD code and stress-imaging technique are suitable for treating the free-surface boundary conditions for Love-wave modeling, (2) Love-wave inversion should be treated with extra care when a LVL exists because of a lack of LVL information in dispersions aggravating uncertainties in the inversion procedure, and (3) energy of high modes in a low-frequency range is very weak, so that it is difficult to estimate the cutoff frequency accurately, and “mode-crossing” occurs between the second higher and third higher modes when a HVL exists. 相似文献
15.
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. 相似文献
16.
Data-resolution Matrix and Model-resolution Matrix for Rayleigh-wave Inversion Using a Damped Least-squares Method 总被引:1,自引:0,他引:1
Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface
geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might
match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and
a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (≥2 Hz) Rayleigh-wave phase
velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated
and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there
are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides
insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure.
Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities
within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because
of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that
selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution
matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices
in these examples to show the potential of increasing model resolution with selected surface-wave data. 相似文献
17.
Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media 总被引:1,自引:0,他引:1
Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity–stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10–16% when the frequency is above 10 Hz due to the velocity dispersion of P and S waves. 相似文献
18.
Yinhe Luo Jianghai Xia Jiangping Liu Yixian Xu Qingsheng Liu 《Soil Dynamics and Earthquake Engineering》2009
The multichannel analysis of surface wave (MASW) method has been effectively used to determine near-surface shear- (S-) wave velocity. Estimating the S-wave velocity profile from Rayleigh-wave measurements is straightforward. A three-step process is required to obtain S-wave velocity profiles: acquisition of a multiple number of multichannel records along a linear survey line by use of the roll-along mode, extraction of dispersion curves of Rayleigh waves, and inversion of dispersion curves for an S-wave velocity profile for each shot gather. A pseudo-2D S-wave velocity section can be generated by aligning 1D S-wave velocity models. In this process, it is very important to understand where the inverted 1D S-wave velocity profile should be located: the midpoint of each spread (a middle-of-receiver-spread assumption) or somewhere between the source and the last receiver. In other words, the extracted dispersion curve is determined by the geophysical structure within the geophone spread or strongly affected by the source geophysical structure. In this paper, dispersion curves of synthetic datasets and a real-world example are calculated by fixing the receiver spread and changing the source location. Results demonstrate that the dispersion curves are mainly determined by structures within a receiver spread. 相似文献
19.
High-frequency surface-wave analysis methods have been effectively and widely used to determine near-surface shear (S) wave velocity. To image the dispersion energy and identify different dispersive modes of surface waves accurately is one of key steps of using surface-wave methods. We analyzed the dispersion energy characteristics of Rayleigh and Love waves in near-surface layered models based on numerical simulations. It has been found that if there is a low-velocity layer (LVL) in the half-space, the dispersion energy of Rayleigh or Love waves is discontinuous and ‘‘jumping’’ appears from the fundamental mode to higher modes on dispersive images. We introduce the guided waves generated in an LVL (LVL-guided waves, a trapped wave mode) to clarify the complexity of the dispersion energy. We confirm the LVL-guided waves by analyzing the snapshots of SH and P–SV wavefield and comparing the dispersive energy with theoretical values of phase velocities. Results demonstrate that LVL-guided waves possess energy on dispersive images, which can interfere with the normal dispersion energy of Rayleigh or Love waves. Each mode of LVL-guided waves having lack of energy at the free surface in some high frequency range causes the discontinuity of dispersive energy on dispersive images, which is because shorter wavelengths (generally with lower phase velocities and higher frequencies) of LVL-guided waves cannot penetrate to the free surface. If the S wave velocity of the LVL is higher than that of the surface layer, the energy of LVL-guided waves only contaminates higher mode energy of surface waves and there is no interlacement with the fundamental mode of surface waves, while if the S wave velocity of the LVL is lower than that of the surface layer, the energy of LVL-guided waves may interlace with the fundamental mode of surface waves. Both of the interlacements with the fundamental mode or higher mode energy may cause misidentification for the dispersion curves of surface waves. 相似文献
20.
Love and Rayleigh Wave Tomography of the Qinghai-Tibet Plateau and Surrounding Areas 总被引:8,自引:0,他引:8
Surface wave data were initially collected from events of magnitude Ms ≥ 5.0 and shallow or moderate focal depth occurred between 1980 and 2002: 713 of them generated Rayleigh waves and 660 Love waves, which were recorded by 13 broadband digital stations in Eurasia and India. Up to 1,525 source-station Rayleigh waveforms and 1,464 Love wave trains have been processed by frequency-time analysis to obtain group velocities. After inverting the path-averaged group times by means of a damped least-squares approach, we have retrieved location-dependent group velocities on a 2° × 2°-sized grid and constructed Rayleigh- and Love-wave group velocity maps at periods 10.4–105.0 s. Resolution and covariance matrices and the rms group velocity misfit have been computed in order to check the quality of the results. Afterwards, depth-dependent SV- and SH-wave velocity models of the crust and upper mantle are obtained by inversion of local Rayleigh- and Love-wave group velocities using a differential damped least-squares method. The results provide: (a) Rayleigh- and Love-wave group velocities at various periods; (b) SV- and SH-wave differential velocity maps at different depths; (c) sharp images of the subducted lithosphere by velocity cross sections along prefixed profiles; (d) regionalized dispersion curves and velocity-depth models related to the main geological formations. The lithospheric root presents a depth that can be substantiated at ~140 km (Qiangtang Block) and exceptionally at ~180 km in some places (Lhasa Block), and which exhibits laterally varying fast velocity very close to that of some shields that even reaches ~4.8 km/s under the northern Lhasa Block and the Qiangtang Block. Slow-velocity anomalies of 7–10% or more beneath southern Tibet and the eastern edge of the Plateau support the idea of a mechanically weak middle-to-lower crust and the existence of crustal flow in Tibet. 相似文献