一种基于层析成像技术提高浅地表面波勘探水平分辨率的方法

在高频面波方法中,水平分辨率是指水平方向上分辨异常体的能力.异常体在水平方向上的长度可用水平方向上横波速度的异常尺度来确定.面波多道分析（MASW）方法被广泛应用于浅地表横波速度结构的探测,然而该方法确定的横波速度是整个检波器排列的平均计算结果,因此水平分辨率较差.另外,采用共中心点（CMP）多次覆盖的方式采集数据亦增加了野外的工作量.我们在MASW方法的基础上,应用面波层析成像方法,提出一套提高面波勘探水平分辨率的完整方法的技术流程.首先,利用波场分离技术获得准确的基阶或高阶模式面波,采用相位扫描的互相关方法测量多道面波记录中任意两道之间的面波走时;然后根据面波层析成像方法,获得高分辨率的各目标网格内的纯路径相速度频散曲线;最后反演所有目标网格内的纯路径相速度频散曲线,得到研究区域的拟二维横波速度结构.这套方法具有一定的抗噪能力,理论上它可以准确地提取相邻两道之间面波的相速度频散曲线;同时由于该方法最少只需要1个排列就可以获得拟二维横波速度结构,因此它显著减小了野外工作量.理论模型和实际资料都证实了这套方法可有效提高面波勘探的水平分辨率.     

Dipping-interface Mapping Using Mode-separated Rayleigh Waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT.     

Research on the middle-of-receiver-spread assumption of the MASW method

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.     

Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure

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.     

S-wave velocity structure in Tangshan earthquake region and its adjacent areas from joint inversion of receiver functions and surface wave dispersion*

Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas (39°N–41.5°N, 115.5°E–119.5°E), we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion. We obtained some detailed information about the Tangshan earthquake region and its adjacent areas, including sedimentary thickness, Moho depth, and crustal and upper mantle S-wave velocity. Meanwhile, we also obtained the v_P/v_S structure along two sections across the Tangshan region. The results show that: (1) the Moho depth ranges from 30 km to 38 km, and it becomes shallower from Yanshan uplift area to North China basin; (2) the thickness of sedimentary layer ranges from 0 km to 3 km, and it thickens from Yanshan uplift region to North China basin; (3) the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure, while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust. Compared with the upper crust, the heterogeneity of the middle and lower crust is more obvious; (4) the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust, and the low v_S and high v_P/v_S beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.     

基于面波频散的三维横波速度方位各向异性层析成像方法

地震各向异性是反映地球内部介质特性的重要指针之一。常用的横波分裂法和二维面波方位各向异性层析成像方法很难准确反映各向异性随深度的变化。将与周期相关的区域化面波方位各向异性转换成与深度相关的一维横波速度方位各向异性可以弥补深度信息不足的缺陷。现有三维横波速度各向异性研究多是通过两步方法来实现的,即逐个周期二维面波方位各向异性层析成像以及逐个格点一维横波速度方位各向异性反演。这种分步反演的方式既不利于三维先验约束的引入,也不利于利用原始观测拟合误差对三维模型进行直接评估。因此本文开发了基于面波频散曲线的三维横波速度方位各向异性层析成像方法,并编制了相关正演和反演程序。为了检测方法和程序的有效性,我们对规律分布的三维检测板模型进行了模拟测试。测试结果显示:该方法可以很好地恢复各向同性波速异常、各向异性相对强度和快波方向等三维结构信息;而且反演模型相对于参考模型明显改善了对观测数据的拟合,降低了对观测数据的均方根误差。但对各向同性理论模型进行各向异性反演时,在波速均匀区可产生小于0.5%的假各向异性幅值,在波速非均匀区该假的各向异性幅值会更大,浅部可达3.5%。因此在实际应用中需要谨慎解释（浅部）非均匀区的各向异性结果。      

A deep-learning-based approach for seismic surface-wave dispersion inversion (SfNet) with application to the Chinese mainland

Surface-wave tomography is an important and widely used method for imaging the crust and upper mantle velocity structure of the Earth. In this study, we proposed a deep learning (DL) method based on convolutional neural network (CNN), named SfNet, to derive the v_S model from the Rayleigh wave phase and group velocity dispersion curves. Training a network model usually requires large amount of training datasets, which is labor-intensive and expensive to acquire. Here we relied on synthetics generated automatically from various spline-based v_S models instead of directly using the existing v_S models of an area to build the training dataset, which enhances the generalization of the DL method. In addition, we used a random sampling strategy of the dispersion periods in the training dataset, which alleviates the problem that the real data used must be sampled strictly according to the periods of training dataset. Tests using synthetic data demonstrate that the proposed method is much faster, and the results for the v_S model are more accurate and robust than those of conventional methods. We applied our method to a dataset for the Chinese mainland and obtained a new reference velocity model of the Chinese continent (ChinaVs-DL1.0), which has smaller dispersion misfits than those from the traditional method. The high accuracy and efficiency of our DL approach makes it an important method for v_S model inversions from large amounts of surface-wave dispersion data.     

Velocity and density contrasts across the Moho interface of Guangdong province in south China constrained by receiver functions*

The P receiver function includes P-to-SV converted phases and multiple reverberations of the discontinuities in the crust and mantle. The time of these phases is related to the crustal thickness and v_P/v_S ratio, and the amplitude of these phases is mainly controlled by the velocity and density contrast of interfaces. By using H-κ stacking method, this work estimated the crustal thickness and v_P/v_S ratio beneath the stations in the Guangdong province of South China. The velocity and density contrast (δβ-δρ) scanning stacking algorithm of the receiver function is applied to constrain the velocity and density contrast of the Moho in Guangdong province. This work analyzed the results of the crustal thickness, v_P/v_S ratio, and the velocity and density contrasts of Moho. The results indicate that the velocity contrast is higher beneath Yangjiang area in western Guangdong province and Nanao area in eastern Guangdong, which has a strong correlation with the distribution of geothermal springs in local areas and the characteristics of high heat flow. The velocity contrast of Moho has also a good correlation with the v_P/v_S ratio and the crustal thickness, which indicates that there is a strong material composition contrasts of the Moho in the study area. Velocity and density contrasts of Moho in some local area (such as western Guangdong) are somewhat consistent with the seismic activities.     

3D <Emphasis Type="Italic">v</Emphasis><Subscript>P</Subscript> and <Emphasis Type="Italic">v</Emphasis><Subscript>S</Subscript> models of southeastern margin of the Tibetan plateau from joint inversion of body-wave arrival times and surface-wave dispersion data

A new 3D velocity model of the crust and upper mantle in the southeastern (SE) margin of the Tibetan plateau was obtained by joint inversion of body- and surface-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the v _S model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 M _S 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes (M ≥ 5) in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.     

Three-dimensional velocity structure around the focal area of the 2021 MS6.4 Yangbi earthquake

On May 21, 2021, an M_S6.4 earthquake occurred in Yangbi, Yunnan province, China, which exhibited typical foreshock-mainshock-aftershock characteristics. To better understand the velocity structure of the focal area and adjacent fault zones, Pg/Sg travel times at 12 seismic stations for the local earthquakes with M_L ≥ 1.5 from 2009–2019 and the Yangbi sequence in May of 2021 were used to invert the three-dimensional (3D) structures for both v_P and v_P/v_S. The obtained structure extends deeply to 15 km for area (25°N–26.5°N, 99.5°E–101°E) at a horizontal resolution of 10× 10 km, and the accuracy of the v_P velocity was verified using airgun signals excited by the Binchuan Airgun Transmitting Seismic Station (BATSS). The resulting v_P and v_P/v_S images correlate with existing fault zones and the Yangbi sequence, including: (1) The shallow velocity structure at 0 km agrees with local topography, where the Binchuan basin exhibits low-v_P and high-v_P/v_S values. From 3–15 km, v_P and v_P/v_S show variations, and the boundaries are consistent with the main faults (e.g., the Weixi-Qiaohou-Weishan, Honghe, and Chenghai faults). (2) The largest foreshock (M_S5.6), mainshock (M_S6.4), and largest aftershock (M_S5.2) occurred near the boundaries where both v_P and v_P/v_S have clear contrasts. (3) Small earthquakes are also concentrated in the transition zone between high- and low-v_P and v_P/v_S anomalies, and are biased toward low-v_P/v_S zones. (4) Boundaries in v_P and v_P/v_S are observed at 20 km west of the Weixi-Qiaohou-Weishan fault, indicating that there may exist one hidden fault.     

Joint inversion of high-frequency surface waves with fundamental and higher modes

Joint inversion of multimode surface waves for estimating the shear (S)-wave velocity has received much attention in recent years. In this paper, we first analyze sensitivity of phase velocities of multimodes of surface waves for a six-layer earth model, and then we invert surface-wave dispersion curves of the theoretical model and a real-world example. Sensitivity analysis shows that fundamental mode data are more sensitive to the S-wave velocities of shallow layers and are concentrated on a very narrow frequency band, while higher mode data are more sensitive to the parameters of relatively deeper layers and are distributed over a wider frequency band. These properties provide a foundation of using a multimode joint inversion to define S-wave velocities. Inversion results of both synthetic data and a real-world example demonstrate that joint inversion with the damped least-square method and the singular-value decomposition technique to invert high-frequency surface waves with fundamental and higher mode data simultaneously can effectively reduce the ambiguity and improve the accuracy of S-wave velocities.     

A physical model for shear-wave velocity prediction1

The clay-sand mixture model of Xu and White is shown to simulate observed relationships between S-wave velocity (or transit time), porosity and clay content. In general, neither S-wave velocity nor S-wave transit time is a linear function of porosity and clay content. For practical purposes, clay content is approximated by shale volume in well-log applications. In principle, the model can predict S-wave velocity from lithology and any pair of P-wave velocity, porosity and shale volume. Although the predictions should be the same if all measurements are error free, comparison of predictions with laboratory and logging measurements show that predictions using P-wave velocity are the most reliable. The robust relationship between S- and P-wave velocities is due to the fact that both are similarly affected by porosity, clay content and lithology. Moreover, errors in the measured P-wave velocity are normally smaller than those in porosity and shale volume, both of which are subject to errors introduced by imperfect models and imperfect parameters when estimated from logs. Because the model evaluates the bulk and shear moduli of the dry rock frame by a combination of Kuster and Toksöz’ theory and differential effective medium theory, using pore aspect ratios to characterize the compliances of the sand and clay components, the relationship between P- and S-wave velocities is explicit and consistent. Consequently the model sidesteps problems and assumptions that arise from the lack of knowledge of these moduli when applying Gassmann's theory to this relationship, making it a very flexible tool for investigating how the v_P-v_s relationship is affected by lithology, porosity, clay content and water saturation. Numerical results from the model are confirmed by laboratory and logging data and demonstrate, for example, how the presence of gas has a more pronounced effect on P-wave velocity in shaly sands than in less compliant cleaner sandstones.     

面波频散谱多模式高分辨率成像的多道信号比较法

高分辨率的面波频散谱成像是浅层地震勘探领域基于频散性质反演横波速度结构中的一个关键步骤.在天然地震探测领域,仅利用两个台站记录的线性信号比较法(LSC),被广泛用来计算面波的频散谱,并用于大尺度的面波层析成像.然而互相关的成像方式会造成频散谱在低频端较低的分辨率.非线性信号比较法(NLSC)利用指数函数克服了这个问题,同时极大地提高了频散谱的成像分辨率.然而,在研究中我们发现,仅利用两个台站的地震记录,并不能将面波的频散特性完整地考虑在内,导致LSC和NLSC方法对高阶模式的成像存在较大的误差.由于主动源面波勘探多道采集的方式,基于信号比较理论的多道信号比较法(MSC)充分利用多道地震信号,可以获得准确的多模式成像,然而该方法需要计算任意两道的频散谱并叠加,存在冗余的计算,导致计算效率较低.因此,本研究对MSC方法进行了相应的改进,通过追踪炮集记录上的面波波组提高了原方法的计算效率,同时,利用理论频散曲线进行叠加分析,验证了改进的MSC方法的正确性和有效性.通过与相移法、LSC和NLSC方法的对比分析,展示了MSC方法是一种准确的、高分辨率的面波多模式频散谱成像方法.实际地震资料的应用...     

利用瑞雷面波反演华北克拉通东北部边界的岩石圈结构

运用瑞雷面波相速度频散曲线分析和反演得到了华北克拉通东北部边界及其邻近区域岩石圈的精细S波速度结构.利用11个地震事件、60个台站的瑞雷面波波形资料,得到了周期从25-150 S的相速度频散曲线,并且通过线性反演方法得到了深度从40-300 km的S波速度结构.结果表明,该研究区域S波速度存在强烈的不均匀性.从东南部的...     

Joint Inversion of Body-Wave Arrival Times and Surface-Wave Dispersion for Three-Dimensional Seismic Structure Around SAFOD

We incorporate body-wave arrival time and surface-wave dispersion data into a joint inversion for three-dimensional P-wave and S-wave velocity structure of the crust surrounding the site of the San Andreas Fault Observatory at Depth. The contributions of the two data types to the inversion are controlled by the relative weighting of the respective equations. We find that the trade-off between fitting the two data types, controlled by the weighting, defines a clear optimal solution. Varying the weighting away from the optimal point leads to sharp increases in misfit for one data type with only modest reduction in misfit for the other data type. All the acceptable solutions yield structures with similar primary features, but the smaller-scale features change substantially. When there is a lower relative weight on the surface-wave data, it appears that the solution over-fits the body-wave data, leading to a relatively rough V _s model, whereas for the optimal weighting, we obtain a relatively smooth model that is able to fit both the body-wave and surface-wave observations adequately.     

Crustal and upper-mantle structure of the Baltic Shield investigated by a combined interpretation of deep-seismic-sounding data and surface-wave analysis

P- and S-wave velocity distributions obtained from DSS measurement have been used as a constraint in the inversion of surface-wave dispersion data. The combined interpretation was made as an attempt to enhance the resulting velocity models and to test the possibility to draw conclusions about the density distribution. The result indicates a potential value of a combined interpretation but it is obvious that very accurate velocity distributions are needed. The achieved density distribution is in good agreement with reported densities derived from gravimetric studies.     

Calculation of Rayleigh-wave phase velocities due to models with a high-velocity surface layer

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.     

Effect of lateral heterogeneity on surface wave testing: Numerical simulations and a countermeasure

The surface wave (S-wave) method has gained popularity in engineering practice for determining S-wave velocity depth profiles. A growing trend is towards the application of S-wave testing for spatially 2-D S-wave velocity tomography, ignoring the assumption of horizontally layered medium. A fourth-order velocity-stress finite difference method is used to perform numerical simulations of S-wave testing in earth models with lateral variation. Results show that the lateral heterogeneity induces a non-stationary property in the space domain, resulting in false depth-related dispersion or higher modes if conventional approach based on stationary assumption is used for the dispersion analysis. Artifacts maybe introduced in spatially 2-D S-wave velocity imaging if the effect of lateral heterogeneity is not accounted for. As a potential countermeasure, a high-lateral-resolution S-wave method is proposed to reduce the effect of lateral heterogeneity while maintaining the resolution and depth range of dispersion analysis. It consists of a walk-away survey and a phase-seaming procedure when synthesizing seismograms with different nearest source-to-receiver offset, allowing wide-wavelength dispersion analysis within a small spatial range. The proof of concept is given with several numerical examples. They show that the high-resolution S-wave method can greatly alleviate the effect of lateral heterogeneity and increase spatial resolution.     

利用微动勘察方法探测煤矿陷落柱

陷落柱是危害煤矿安全生产的主要因素之一,如何有效探测陷落柱是实现煤矿安全、高效生产的当务之急.本文首次将微动勘察方法应用于煤矿采区勘探,结果表明,该方法对陷落柱异常区反映敏感.在山西潞安漳村煤矿2002工作面测区已知陷落柱上,无论是单点反演获得的S波速度结构,还是微动视S波速度剖面,均能清晰显示陷落柱.微动剖面确定的陷落柱位置与巷道揭露位置一致,边界误差在10 m左右.微动勘察方法是探测陷落柱的一种行之有效的物探方法,由于其分辨率高、无需人工源、野外施工方法灵活便捷、不受施工场地限制等特点,对探测村庄覆盖区之下的煤层构造、圈定陷落柱等,具有得天独厚的技术优势和应用前景.     

The Selection of Field Acquisition Parameters for Dispersion Images from Multichannel Surface Wave Data

The accuracy and resolution of surface wave dispersion results depend on the parameters used for acquiring data in the field. The optimized field parameters for acquiring multichannel analysis of surface wave (MASW) dispersion images can be determined if preliminary information on the phase velocity range and interface depth is available. In a case study on a fill slope in Hong Kong, the optimal acquisition parameters were first determined from a preliminary seismic survey prior to a MASW survey. Field tests using different sets of receiver distances and array lengths showed that the most consistent and useful dispersion images were obtained from the optimal acquisition parameters predicted. The inverted S-wave velocities from the dispersion curve obtained at the optimal offset distance range also agreed with those obtained by using direct refraction survey.