基于最佳小波基的地震面波插值方法

在利用实际地震数据中的面波反演近地表横波速度的过程中,若道间距较大、空间采样率不足,则会产生空间假频现象,从而降低频率速度谱的信噪比,影响频散曲线提取的精度以及反演效果,因此需要针对面波进行插值处理。文中提出了一种基于最佳小波基的地震面波插值方法,通过理论分析和实验误差对比在地震数据处理常用的众多小波基中选出适用于插值处理的最佳小波基bior6.8,提高了插值精度。针对面波同向轴为线性且斜率较大的特点,文中首先采用线性动校正的方法对面波进行拉平处理,再进行小波变换插值,最后进行反线性动校正恢复面波。通过对理论模型与实际资料进行插值处理验证了本文方法的有效性,插值后的面波记录波形恢复较好,显著提高了频率速度谱的信噪比,有效解决了面波数据空间采样率不足引起的假频问题。     

青藏高原东部基于噪声的面波群速度分布特征

通过收集青海、甘肃、四川三省的76个地震台记录的2008年1—12月三分量的连续噪声数据,利用噪声面波层析成像的方法获得了青藏高原东部的面波群速度分布特征。首先采用多重滤波方法提取了1 000多条台站对5～50 s的三分量面波群速度频散曲线,然后将研究区域划分为0.2°×0.2°的网格,利用O ccam方法反演了瑞利波(R-R)和勒夫波(T-T)的群速度分布。反演得到的群速度分布特征与地表地质和构造特征表现出较好的相关性,清晰地揭示了地壳内部的横向速度变化。层析成像的结果显示在短周期(8～20 s)内,拥有较厚的沉积层的四川盆地表现为明显的低速特征,而青藏高原东部则表现为较高的群速度分布特征;随着周期的增加(>20 s),群速度的分布特征呈现出与短周期相反的特性,青藏高原东部下方的速度远远低于四川盆地,这可能与青藏高原东部中、下地壳低速层相关联,同时也意味着研究区域的地壳结构具有明显的横向不均匀性。在群速度分布图上,龙门山不仅是四川盆地与青藏高原的地形和构造分界带,同时也对应着高群速度与低群速度的过渡带。     

Development and application of a GIS-based tool for earthquake-induced hazard prediction

An integrated GIS-based tool (GTIS) was constructed to estimate site effects related to the earthquake hazards in the Gyeongju area of Korea. To build the GTIS for the study area, intensive site investigations and geotechnical data collections were performed and a walk-over site survey was additionally carried out to acquire surface geo-knowledge data in accordance with the procedure developed to build the GTIS. For practical applications of the GTIS used to estimate the site effects associated with the amplification of ground motion, seismic microzoning maps of the characteristic site period and the mean shear wave velocity to a depth of 30 m were created and presented as a regional synthetic strategy addressing earthquake-induced hazards. Additionally, based on one-dimensional site response analyses, various seismic microzoning maps for short- and mid-period amplification potentials were created for the study area. Case studies of seismic microzonations in the Gyeongju area verified the usefulness of the GTIS for predicting seismic hazards in the region.     

Seismic interferometry and ambient noise tomography in the British Isles

Traditional methods of imaging the Earth's subsurface using seismic waves require an identifiable, impulsive source of seismic energy, for example an earthquake or explosive source. Naturally occurring, ambient seismic waves form an ever-present source of energy that is conventionally regarded as unusable since it is not impulsive. As such it is generally removed from seismic data and subsequent analysis. A new method known as seismic interferometry can be used to extract useful information about the Earth's subsurface from the ambient noise wavefield. Consequently, seismic interferometry is an important new tool for exploring areas which are otherwise seismically quiescent, such as the British Isles in which there are relatively few strong earthquakes. One of the possible applications of seismic interferometry is ambient noise tomography (ANT). ANT is a way of using interferometry to image subsurface seismic velocity variations using seismic (surface) waves extracted from the background ambient vibrations of the Earth. To date, ANT has been used successfully to image the Earth's crust and upper-mantle on regional and continental scales in many locations and has the power to resolve major geological features such as sedimentary basins and igneous and metamorphic cores. Here we provide a review of seismic interferometry and ANT, and show that the seismic interferometry method works well within the British Isles. We illustrate the usefulness of the method in seismically quiescent areas by presenting the first surface wave group velocity maps of the Scottish Highlands using only ambient seismic noise. These maps show low velocity anomalies in sedimentary basins such as the Moray Firth, and high velocity anomalies in igneous and metamorphic centres such as the Lewisian complex. They also suggest that the Moho shallows from south to north across Scotland which agrees with previous geophysical studies in the region.     

Distinct upper mantle deformation of cratons in response to subduction: Constraints from SKS wave splitting measurements in eastern China

Interaction between the subducting slab, the overriding continental lithosphere and mantle flow are fundamental geodynamic processes of subduction systems. Eastern China is an ideal natural laboratory to investigate the behavior and evolution of cratonic blocks within a subduction system. In this study, we investigate deformation of the upper mantle beneath eastern China. We present seismic shear wave splitting measurements from three networks consisting of over 483 broadband stations, with 157 stations giving a total of 516 results. The splitting parameters exhibit complex regional patterns but are relatively coherent within individual tectonic units. Tectonic blocks exhibited distinctive fast directions relative to regional features. The dominant attitude of fast directions for the North China Craton was subparallel to the direction of subduction, whereas fast directions for Southeastern China were perpendicular to the direction of subduction. The shear wave splitting measurements were interpreted according to a high resolution tomographic body-wave velocity model. Combining these two datasets showed that the predominant geodynamic models for the region (mantle plume, mantle wedge and flat-slab subduction models) are incompatible with the observations presented here. We suggest that the North China Craton, Yangtze Craton and the Cathaysia block have undergone different deformational events due to differing mantle flow patterns, and distinct spatial and temporal subduction histories of the Pacific and Philippine Sea plates.     

基于小波包变换的面波分离技术研究

依据面波视速度的分布范围,提出一种新的叠前面波分离技术,即把地震波场分成不受面波影响的波场和受面波影响的波场,然后利用小波包变换对受面波影响的波场进行时频分析,通过分解与重构,从中提取有效波,再将其与不受面波影响的波场合并,就可得到面波分离后的地震波场。该方法在成功分离面波的同时,可使有效波的能量(尤其是其低频能量)不受损失,既有效地保持信号的频宽不变,又提高了记录的信噪比。实际资料的处理结果表明:此方法在叠前面波分离方面具有有良好的效果。     

Synthetic seismogram and surface wave constraints on crustal models of Spitsbergen

Seismic refraction data collected on Spitsbergen in 1978 are used to obtain a crustal model assuming plane horizontal layering. The observed travel-times and wave forms are compared with those of synthetic seismograms computed for various published crustal models. The more detailed models adequately explain some, but not all, of the features of the synthetics. These models are adjusted, utilizing travel times and wave-form amplitudes until a satisfactory fit is achieved. The best-fitting model consists of a 4-layer crust having thicknesses of 4.1, 10.0 7.4 and 5.8 km and compressional velocities of 4.65, 6.21, 6.30 and 6.65 km/sec with increasing depth. The uppermost mantle has a velocity of 7.90 km/sec. A comparison of observed and synthetic P_n waveforms supports the existence of a thin low-velocity zone beginning at a depth of about 5 km beneath the Moho boundary. An inversion of seismic surface wave group velocity data yields a shear-wave model which is compatible with the compressional wave model.     

地质条件复杂地区瑞雷波勘探正演模拟

为研究复杂地质条件对瑞雷面波频散曲线的影响,采用有限差分地震数值模拟方法,对3种模型进行正演模拟,得出了某一时刻的波场快照与单炮记录。对模拟记录进行频散分析的结果表明:单层均匀介质没有频散现象;起伏地表两层均匀介质的频散曲线相速度只和介质本身的参数有关,不受地形的影响;水平地表含软夹层介质中,瑞雷波频散现象非常明显。      

反射槽波探测技术中速度分析方法研究

在巷道掘进过程中,利用反射槽波探测技术探查巷道侧帮构造发育情况,并准确解释槽波波阻抗反射界面的位置对矿井生产意义重大,但现阶段反射槽波地震解释断层的位置常存在偏差,主要原因是数据处理阶段选取速度存在较大误差。为降低反射槽波速度选取误差,系统分析国内13个煤矿15个煤层的槽波速度特征,总结出反射槽波速度与直达槽波、透射槽波速度的关系。研究发现反射槽波速度高于直达槽波速度3%~12%;低于透射槽波速度1%~6%。将其研究成果引入山西某矿反射槽波探测项目,槽波解释断层位置与探巷揭露位置偏差仅为3 m,取得了较为理想的应用效果。该研究为进一步提高反射槽波探测精度提供了新思路。      

Application of sequential Gaussian simulation to quantify the surface wave magnitude uncertainty within the faulted regions

Geostatistical simulations have been recently widely used in the geological and mining investigations. Variogram, the fundamental tools of geostatistics, can identify the spatial distribution of the regionalized variable within the area. One of the important issues of geostatistical simulation in seismotectonics is producing uncertainty maps, which could be applicable to predict earthquake parameters through the site locations especially for civil structures like bridges. It can help engineers to design the structure of interest better. Earthquake parameters as for example seismic fault and surface wave magnitude (Ms) have significant impact on the feasibility study of the civil structures. In this research, a method is presented to produce uncertainty maps for seismic fault and surface wave magnitude, Ms. For this aim, information related to surface wave magnitude and fault trace in Zagros region (SW of Iran) has been collected. Then, the relationships between them through the site location have been investigated and analyzed by conditional geostatistical simulation. In order to quantify the uncertainty of each parameter, the uncertainty formula after generating the E-type maps has been provided and discussed. Finally, in “Talgah Bridge” site, these uncertainty maps were produced to interpret the impact of the surface wave magnitude and fault trace in this specific civil structure.     

Local modification of the lithosphere beneath the central and western North China Craton: 3-D constraints from Rayleigh wave tomography

We have imaged the lithospheric structure beneath the central and western North China Craton (NCC) with Rayleigh wave tomography. The Rayleigh waveforms of 100 teleseismic events recorded by 208 broadband stations are used to yield high-resolution phase velocity maps at 13 periods from 20 s to 143 s. A 3-D S-wave velocity model is constructed based on the phase velocity maps. Our S-wave velocity model is broadly consistent with the results of previous tomography studies, but shows more detailed variations within the lithosphere. The Trans-North China Orogen (TNCO) is generally characterized by low-velocity anomalies but exhibits great heterogeneities. Two major low-velocity zones (LVZs) are observed in the north and south, respectively. The northern LVZ laterally coincides with sites of Cenozoic magmatism and extends to depths greater than 200 km. We propose that a small-scale mantle upwelling is present, confined to the north of the TNCO. A high-velocity patch in the uppermost mantle is also observed between the two LVZs adjacent to the narrow transtensional zone of the Cenozoic Shanxi–Shaanxi Rift (SSR). We interpret this as the remnant of a cratonic mantle root. The Ordos Block in the western NCC is associated with high-velocity anomalies, similarly reflecting the existence of cratonic mantle root, but a discernible low-velocity layer is observed at depths of 100–150 km in this location. We interpret that this mid-lithospheric structure was probably formed by metasomatic processes during the early formation of the NCC. Based on the observations from our S-wave velocity model, we conclude that the current highly heterogeneous lithospheric structure beneath the TNCO is the result of multiphase reworking of pre-existing mechanically weak zones since the amalgamation of the craton. The latest Cenozoic lithospheric reworking is dominated by the far-field effects of both Pacific plate subduction and the India–Eurasia collision.     

Geospatial contour mapping of shear wave velocity for Mumbai city

Shear wave velocity is one of the most important input parameter in the analysis of geotechnical earthquake engineering problems, particularly to estimate site-specific amplification factor and ground response study. Dynamic in situ tests such as spectral analysis of surface waves (SASW) or multichannel analysis of surface waves (MASW) are very expensive. Also due to lack of specialized personnel, these tests are generally avoided in many soil investigation programs. Worldwide, several researchers have developed correlations between the SPT ‘N’ value and shear wave velocity ‘V _s’, which are useful for determining the dynamic soil properties. In the present study, more than 400 numbers of soil borehole data were collected from various geotechnical investigation agencies, government engineering institutes and geotechnical laboratories from different parts of Mumbai city, which is financial capital of India with highest population density. In this paper, an attempt has been made to develop the correlation between the SPT ‘N’ value and shear wave velocity ‘V _s’ for various soil profile of Mumbai city and compared with other existing correlations for different cities in India. Using Geographical Information System (GIS), a geospatial contour map of shear wave velocity profile for Mumbai city is prepared with contour intervals of 25 and 50 m/s. The scarcity of database or maps of shear wave velocity profile for Mumbai city will make the present geospatial contour maps extremely useful and beneficial to the designer, practitioners for seismic hazard study involved in geotechnical earthquake engineering.     

均匀介质中F-K反偏移法正演模拟

F-K反偏移法正演模拟是在给定地质模型后通过构建人工偏移场,再对人工偏移场进行F-K反偏移来实现。人工偏移场的构建是在地质模型上直接放上经正确拉伸和幅度校正后的地震脉冲子波。F-K反偏移中插值映射采用两点sinc插值,其能避免假同相轴的出现且具有较快的速度。奇异点处理使用泰勒近似展开。最后对凹陷模型进行了F-K反偏移法正演模拟实现,得到了预期的效果。     

The Chelyabinsk meteoroid: A seismologist’s view

This work determines the parameters of the meteoroid that exploded in Chelyabinsk oblast on February 15, 2013. Using video observation and automobile registration cameras, we identified the trajectory and velocity of the meteoric body from the beginning of its glow to the destruction. The precise explosion time (03:20:34UTC) was determined on the basis of records of the broadband seismic stations and arrivals of the surface wave. The estimated explosion energy with assumption of a point source similar to a high thermonuclear explosion was 100 kt in TNT equivalent based on the surface wave amplitude. The archive of records of seismic stations during the meteorite flight is composed.     

Surface wave studies for shear wave velocity and bedrock depth estimation over basalts

Shear wave velocity (V _S) estimation is of paramount importance in earthquake hazard assessment and other geotechnical/geo engineering studies. In our study, the shear wave velocity was estimated from ground roll using multichannel analysis of surface wave (MASW) technique making use of dispersive characteristics of Rayleigh type surface waves followed by imaging the shallow subsurface basaltic layers in an earthquake-prone region near Jabalpur, India. The reliability of MASW depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Inversion of data from surface waves resulted in a shear wave velocity (V _S) in the range of 200–1,200 m/s covering the top soil to weathering and up to bedrock corresponding to a depth of 10–30 m. The P-wave velocity (V _P) obtained from refraction seismic studies at these locations found to be comparable with V _S at an assumed specific Poisson’s ratio. A pair of selected set of V _S profiles over basalt which did not result in a hazardous situation in an earthquake of moderate magnitude are presented here as a case study; in other words, the shear wave velocity range of more than 200 m/s indicate that the area is highly unlikely prone to liquefaction during a moderate or strong earthquake. The estimated depth to basalt is found to be 10–12 m in both the cases which is also supported by refraction studies.     

Lithospheric structure of the Tien Shan region constrained by joint inversion of receiver function and Rayleigh wave dispersion

We obtain a lithospheric shear‐wave velocity model across the Tien Shan orogenic belt by jointly inverting Rayleigh wave group velocities and teleseismic P‐wave receiver functions at 61 broadband seismic stations deployed in this region. Our new model reveals prominent lateral variations of shear‐wave velocity in both the crust and uppermost mantle. This model reveals different structures in the upper and middle crust across the Talas Fergana Fault, which may suggest the presence of a tectonic boundary between the western and central Tien Shan beneath the fault. According to the velocity images, the depth extent of the fault is ~40 km and this is confined to the crust. Pronounced low‐velocity anomalies are imaged in the middle crust and uppermost mantle beneath the southern and middle Tien Shan, implying that the upwelling of the materials from the upper mantle could have played an important role in the mountain building.     

Determination of shear wave velocity and depth to basement using multichannel analysis of surface wave technique

The dispersive characteristics of Rayleigh type surface waves were utilized to estimate shear wave velocity (Vs) profile followed by imaging the shallow subsurface granitic layers in the heart of Hyderabad. The reliability of Multichannel Analysis of Surface Waves (MASW) depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Multichannel recording leads to effective identification and isolation of various factors of noise. Calculating the 1-D shear wave velocity (Vs) field from surface waves ensures high degree of accuracy irrespective of cultural noise. The main advantage of mapping the bed rock surface with shear wave velocity (Vs) is the insensitivity of MASW to velocity inversion besides being free from many constraints such as contrast in physical properties etc. Modeling of surface waves data results a shear wave velocity (Vs) of 250?C750 m/s covering the top soil to weathering and up to bedrock corresponding to a depth range of 10?C30 m. Further, the computed N values (which is an indicator of site characteristic) based on the harmonic shear wave velocity up to a depth of 5 m is found to be quite high (> 25?C30) well above 5 indicating the site to be safe and strong enough and not prone to liquefaction. A pair of selected set of results over granites are presented here as a case study highlighting the salient features of MASW.     

密集线性台阵地震背景噪声速度成像及其在湖南沃溪金钨锑矿勘探中的应用

多道面波分析技术在近地表勘探领域有着广泛应用,准确的提取频散曲线成为面波勘探成像的关键。文章介绍了一种新的地震背景噪声互相关面波频散成像方法——拓距相移法。该方法在传统相移法的基础上,利用阵内相移对小孔径范围的面波中高频信号进行提取,并利用阵外相移对大孔径范围的面波中低频信号进行提取,然后将两部分频散曲线融合从而得到更宽频带的面波频散曲线用于地下速度结构的反演。该方法在保证对近地表结构进行较高分辨率成像的同时,大大增加对深部结构的有效约束。2019年9月到10月期间,作者在湖南沃溪布设了8条密集测线,进行了1个月的地震背景噪声数据采集,并利用上述拓距相移法提取了0.1~2 s的瑞利面波宽频带相速度频散曲线。通过初步反演其中3条测线的背景噪声数据,获得了该矿区深度2.5 km以浅的地震横波速度结构。经与已知地质资料比对,160测线的地震横波速度反演结果与断层、岩性分界面及矿脉有着较好的对应关系,表明获得的沃溪矿区地震横波速度结构较好地反映了控矿构造、岩性分界以及矿体的分布位置等信息,为该区中—深部找矿提供了重要依据。该研究利用实际数据检验了拓距相移法的有效性,为今后深部找矿提供了一个有效的高精度成像方法。     

夹层厚度对瑞利面波频散曲线的影响

文章以弹性波横波速度为参数,利用瑞利面波频散曲线的快速标量传递算法进行正演模拟计算,系统研究软弱层厚度对瑞利面波基阶模频散曲线的影响,分析和总结了频率与相速度（f-c）、频率与波数（f-k）、深度与相速度（h-c）等瑞利面波基阶模曲线随软弱层厚度变化的规律。     

Crustal and upper mantle S-wave velocity structure beneath the Bransfield Strait (West Antarctica) from regional surface wave tomography

Regional surface wave tomography in the sub-Antarctic Scotia Sea is helpful in revealing the nature of the crust and the S-wave seismic velocity profile beneath the Bransfield Strait. The joint use of our regional network, global seismographic network stations and local temporary arrays provide better lateral resolution than that obtained in our previous studies concerning the Scotia Sea region.Tomographic analysis of data obtained using 10 broad band seismic stations and more than 300 regional events, shows that the Bransfield Basin is characterised by a strong group velocity reduction of 8% with respect to the surrounding areas, in the period range from 15 s to 50 s.The crustal and upper mantle models of the eastern, central and western Bransfield Basin are obtained by joint inversion of Rayleigh and Love local dispersion curves from 15 s to 50 s. In addition our data set is expanded to a broader period interval (1–80 s), in central Bransfield Strait in order to better constrain the upper mantle and shallow crust.The main results can be summarized as follows: (a) the crust thins distinctly from W toward E; the variation is consistent with the type of volcanism, earthquake distribution and bathymetric observations, (b) low upper mantle velocities (soft lid) extend down to depths exceeding 70 km as a consequence of elevated temperatures, (c) the crust beneath the central Bransfield Basin displays continental characteristics with a gradually increasing S-wave velocity distribution versus depth analogous to the East African Rift structure of Kenya, (d) negative velocity gradients are present in the lower crust beneath the eastern Bransfield Basin; these could be interpreted as magmatic bodies originating from decompression melting of the mantle.