On the repeatability and consistency of three-component ambient vibration array measurements

Ambient vibration measurements with small, temporary arrays that produce estimates of surface wave dispersion have become increasingly popular as a low-cost, non-invasive tool for site characterisation. An important requirement for these measurements to be meaningful, however, is the temporal consistency and repeatability of the resulting dispersion and spatial autocorrelation curve estimates. Data acquired within several European research projects (NERIES task JRA4, SESAME, and other multinational experiments) offer the chance to investigate the variability of the derived data products. The dataset analysed here consists of repeated array measurements, with several years of time elapsed between them. The measurements were conducted by different groups in different seasons, using different instrumentations and array layouts, at six sites in Greece and Italy. Ambient vibration amplitude spectra and locations of dominant sources vary between the two measurements at each location. Still, analysis indicates that this does not influence the derived dispersion information, which is stable in time and neither influenced by the instrumentation nor the analyst. The frequency range over which the dispersion curves and spatial autocorrelation curves can be reliably estimated depends on the array dimensions (minimum and maximum aperture) used in the specific deployment, though, and may accordingly vary between the repeated experiments. The relative contribution of Rayleigh and Love waves to the wavefield can likewise change between repeated measurements. The observed relative contribution of Rayleigh waves is generally at or below 50%, with especially low values for the rural sites. Besides, the visibility of higher modes depends on the noise wavefield conditions. The similarity of the dispersion and autocorrelation curves measured at each site indicates that the curves are stable, mainly determined by the sub-surface structure, and can thus be used to derive velocity information with depth. Differences between velocity models for the same site derived from independently determined dispersion and autocorrelation curves—as observed in other studies—are consequently not adequately explained by uncertainties in the measurement part.     

Imaging a shallow salt diapir using ambient seismic vibrations beneath the densely built-up city area of Hamburg,Northern Germany

Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.     

大别-苏鲁及其邻近地区基于背景噪声的瑞雷波和勒夫波相速度层析成像

采用与作者2014年发表的“大别-苏鲁及其邻近地区基于背景噪声的勒夫波群速度成像”文章相同的资料,用频时分析提取5 000余条瑞雷波和4 000余条勒夫波相速度频散曲线,反演得到了8—32 s的瑞雷波和勒夫波相速度分布图像．结果显示,瑞雷波与勒夫波相速度分布具有很好的一致性．8 s的相速度分布与地表构造特征相吻合,造山带与隆起区均表现为高速,盆地因其规模不同而显示不同程度的低速．随着周期的增大,大别 苏鲁的高速带由强变弱,但始终存在．16—24 s的高速可能主要受到中地壳高速的控制,而32 s的高速则可能与上地幔顶部的高速有关．比较大别造山带与苏鲁造山带的平均频散曲线,发现大别造山带和苏鲁造山带的勒夫波频散曲线均高于AK135模型计算的理论频散曲线,而瑞雷波则没有这一现象. 这可能意味着两个地区有比较强烈的径向各向异性．      

Crustal structure beneath Liaoning province and the Bohai Sea and its adjacent region in China based on ambient noise tomography

The velocity structure of the crust beneath Liaoning province and the Bohai sea in China was imaged using ambient seismic noise recorded by 73 regional broadband stations. All available three-component time series from the 12-month span between January and December 2013 were cross-correlated to yield empirical Green's functions for Rayleigh and Love waves. Phasevelocity dispersion curves for the Rayleigh waves and the Love waves were measured by applying the frequencytime analysis method. Dispersion measurements of the Rayleigh wave and the Love wave were then utilized to construct 2D phase-velocity maps for the Rayleigh wave at8–35 s periods and the Love wave at 9–32 s periods,respectively. Both Rayleigh and Love phase-velocity maps show significant lateral variations that are correlated well with known geological features and tectonics units in the study region. Next, phase dispersion curves of the Rayleigh wave and the Love wave extracted from each cell of the 2D Rayleigh wave and Love wave phase-velocity maps,respectively, were inverted simultaneously to determine the3 D shear wave velocity structures. The horizontal shear wave velocity images clearly and intuitively exhibit that the earthquake swarms in the Haicheng region and theTangshan region are mainly clustered in the transition zone between the low-and high-velocity zones in the upper crust, coinciding with fault zones, and their distribution is very closely associated with these faults. The vertical shear wave velocity image reveals that the lower crust downward to the uppermost mantle is featured by distinctly high velocities, with even a high-velocity thinner layer existing at the bottom of the lower crust near Moho in central and northern the Bohai sea along the Tanlu fault, and these phenomena could be caused by the intrusion of mantle material, indicating the Tanlu fault could be just as the uprising channel of deep materials.     

Love wave tomography in Italy from seismic ambient noise

We estimate Love wave empirical Green's functions from cross-correlations of ambient seismic noise to study the crust and uppermost mantle structure in Italy.Transverse-component ambient noise data from October 2005 through March 2007 recorded at 114 seismic stations from the Istituto Nazionale di Geofisica e Vulcanologia (INGV) national broadband network,the Mediterranean Very Broadband Seismographic Network (MedNet) and the Austrian Central Institute for Meteorology and Geodynamics (ZAMG) yield more than 2 000 Love wave group velocity measurements using the multiple-filter analysis technique.In the short period band (5-20 s),the cross-correlations show clearly one-sided asymmetric feature due to non-uniform noise distribution and high local activities,and in the long period band (20 s) this feature becomes weak owing to more diffusive noise distribution.Based on these measurements,Love wave group velocity dispersion maps in the 8-34 s period band are constructed,then the SH wave velocity structures from the Love wave dispersions are inverted.The final results obtained from Love wave data are overall in good agreement with those from Rayleigh waves.Both Love and Rayleigh wave inversions all reveal that the Po plain basin is resolved with low velocity at shallow depth,and the Tyrrhenian sea is characterized with higher velocity below 8 km due to its thin oceanic crust.     

东北地区背景噪声的Rayleigh和Love波相速度层析成像

本文利用中国数字地震台网位于东北地区的122个宽频地震台站的18个月记录的三分量连续地震噪声数据,采用互相关方法提取了Rayleigh和Love波经验格林函数,并利用时频自动分析技术获取了相应的相速度频散曲线.通过反演频散曲线,获得了Rayleigh和Love波周期为8~35s的二维相速度分布.结果表明,东北地区相速度的分布存在横向和垂向的不均匀性.短周期的相速度分布同地表地质构造密切相关,松辽盆地及山间沉积盆地呈现低速异常,而大兴安岭、小兴安岭及东部的一些山岭显示高速异常.随着周期的增加,位于中间的松辽盆地变为高低速相间,两侧的造山带呈现低速异常.这种异常的转变,可能是受构造活动或者莫霍面深度的影响.另外,在周期为20~35s频段内,Rayleigh和Love波同一周期的相速度在松辽盆地和位于吉林地区的郯庐断裂带表现不一致,表明可能存在径向各向异性.     

微动探测方法软件的实现

利用微动进行勘探无需振动源,方便快捷、对环境没有干扰,是一种很有发展前景的勘探方法。针对微动探测数据处理的需求,开发研制了一套可视化微动探测处理分析软件系统MSMSPAC。该软件采用C~#语言开发,主要功能包括台阵设计、台阵方位角分布分析、数据格式转换、微动有效数据段挑选、数据滤波、重采样、参数计算及频散曲线提取、编辑、合并、项目管理等功能。软件系统界面友好,操作简便。     

Dispersion Energy Analysis of Rayleigh and Love Waves in the Presence of Low-Velocity Layers in Near-Surface Seismic Surveys

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.     

利用地脉动台阵观测推算Love波频散曲线

土层剪切波速结构是工程场地条件的重要的基本参数，许多研究者关注利用从地脉动台阵中推得的Rayleigh波或Love波的相速度频散曲线反演实际工程场地剪切波速度结构。探讨从地脉动台阵观测的记录中提取Love波相速度的频散曲线的方法，从厦门某工程场地的地脉动台阵观测三分量记录中提取了Love波相速度的频散曲线，将提取的结果与根据实际场地资料构成的水平成层介质模型计算的Love波相速度的频散曲线进行了对比，为下一步的反演奠定了基础。     

Characterization of site conditions (soil class,V<Subscript>S30</Subscript>, velocity profiles) for 33 stations from the French permanent accelerometric network (RAP) using surface-wave methods

Data provided by accelerometric networks are important for seismic hazard assessment. The correct use of accelerometric signals is conditioned by the station site metadata quality (i.e., soil class, V_S30, velocity profiles, and other relevant information that can help to quantify site effects). In France, the permanent accelerometric network consists of about 150 stations. Thirty-three of these stations in the southern half of France have been characterized, using surface-wave-based methods that allow derivation of velocity profiles from dispersion curves of surface waves. The computation of dispersion curves and their subsequent inversion in terms of shear-wave velocity profiles has allowed estimation of V_S30 values and designation of soil classes, which include the corresponding uncertainties. From a methodological point of view, this survey leads to the following recommendations: (1) perform both active (multi-analysis surface waves) and passive (ambient vibration arrays) measurements to derive dispersion curves in a broadband frequency range; (2) perform active acquisitions for both vertical (Rayleigh wave) and horizontal (Love wave) polarities. Even when the logistic contexts are sometimes difficult, the use of surface-wave-based methods is suitable for station-site characterization, even on rock sites. In comparison with previous studies that have mainly estimated V_S30 indirectly, the new values here are globally lower, but the EC8-A class sites remain numerous. However, even on rock sites, high frequency amplifications may affect accelerometric records, due to the shallow relatively softer layers.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

Phase and Correlation in `Random' Seismic Fields and the Reconstruction of the Green Function

We first present a summary of recent results on coda interpretation. We emphasize the observation of the stabilization of P to S energy ratio indicating the modal equipartition of the wavefield. This property clearly shows that the coda waves are in the regime of multiple scattering. Numerical solutions of the elastic radiative transfer equation are used to illustrate the evolution of the wave-field towards P-to-S energy stabilization, and asymptotically to complete isotropy. The energy properties of the coda have been widely studied but the phase properties have often been neglected. The recently observed coherent backscattering enhancement, an expression of the so-called `weak localization', demonstrates that interference effects still persist for multiple diffracted waves. Another manifestation of the persistence of the phase is the possibility to reconstruct the Green function between two stations by averaging the cross correlation of coda waves produced by distant earthquakes and recorded at those two stations. This reconstruction is directly related to the properties of reciprocity and time reversal of any wavefield. Using broadband seismic coda waves, we show that the dominant phases of the Green function in the band 2 s–10 s, namely fundamental mode Rayleigh and Love waves, are reconstructed. We analyze the time symmetry of the cross correlation and show how the level of symmetry evolves with the isotropization of the diffuse field with lapse time. Similarly we investigate the correlation in continuous ambient noise records. Whereas the randomness of the coda results from multiple scattering by randomly distributed scatterers, we assume that the seismic noise is random mostly because of the distribution of sources at the surface of the Earth. Surface waves can be extracted from long time series. The dispersion curves of Rayleigh waves are deduced from the correlations. On paths where measurements from earthquake data are also available, we show that they are in good agreement with those deduced from noise correlation. The measurement of velocities from correlation of noise along paths crossing different crustal structures opens the way for a `passive imaging' of the Earth's structure.     

Ocean-generated microseismic noise located with the Gräfenberg array

The main cause for mid-period seismic ground distortions are ocean waves generated by atmospheric disturbances. These act upon the earth through different mechanisms. The microseismic wavefield can be divided into primary (T =12–18 s) and secondary (T = 6–9 s) noise. Classical theory tells that the origin of these induced ground distortions depends on the location and the intensity of the low pressure region. A considerable part of the microseismic wave field reaches the GRF-array in southern Germany with high coherency and almost constant amplitudes. Thus it is possible to locate the generating areas using frequency-wavenumber analysis. Five discrete generating areas for secondary microseisms and three generating areas for primary microseisms could be determined in the Atlantic Ocean, the Arctic Sea and the Mediterranean Sea by investigating broadband continuous recordings over four months in winter 1995/96. An essential result is the long-time constancy of the backazimuths of the coherent part of the microseismic wavefield with respect to the origin areas, independent of the location of the moving low pressure zone. Results from a triangulation using additionally broadband data from the NORSAR-array and an independent estimation of the distance of the source region with water wave dispersion data indicate an origin of the secondary microseismic wavefield near the north-Norwegian coast for the strongest source. The array analysis of a temporary network of ten three-component broadband stations in south-east Germany shows that the ratio of energy between coherent Love and Rayleigh waves is much higher for the primary than for the secondary microseismic noise wavefield. This indicates differences in the source mechanisms.     

Inversion of local S-wave velocity structures from average H/V ratios, and their use for the estimation of site-effects

H/V spectral ratios from microtremors areused to retrieve the S-velocity structurefrom a single ambient vibration record, byusing its relation to the ellipticity ofthe fundamental mode Rayleigh wave and theamplitude of observed H/V ratio.Constraints are needed in order to restrictthe possible range of solutions, and theinversion is applied to sites where thethickness of the unconsolidated sedimentsis approximately known from boreholeinformation. Within the uncertainty, theinverted structures agree well with theresults from other S-wave measuringtechniques such as downhole and cross-holemeasurements, and the analysis of ambientvibrations measured on an array.The influence of the inversion uncertaintyon site-amplification estimates forearthquakes is then investigated. For allinverted models, site response is computedfor a large number of events, which allowsto define the uncertainty by the aprioriunknown source position and mechanism of afuture earthquake. In most cases thevariability between the results obtainedfor the different models is much smallerthan the variability introduced by theunknown source position. The accuracy withwhich S-wave velocity structures can beretrieved from observed H/V ratios istherefore sufficient for an application ofthe method in seismic hazard analysis for aspecific site.     

Time-domain Pure-state Polarization Analysis of Surface Waves Traversing California

-- A time-domain pure-state polarization analysis method is used to characterize surface waves traversing California parallel to the plate boundary. The method is applied to data recorded at four broadband stations in California from twenty-six large, shallow earthquakes which occurred since 1988, yielding polarization parameters such as the ellipticity, Euler angles, instantaneous periods, and wave incident azimuths. The earthquakes are located along the circum-Pacific margin and the ray paths cluster into two groups, with great-circle paths connecting stations MHC and PAS or CMB and GSC. The first path (MHC-PAS) is in the vicinity of the San Andreas Fault System (SAFS), and the second (CMB-GSC) traverses the Sierra Nevada Batholith parallel to and east of the SAFS. Both Rayleigh and Love wave data show refractions due to lateral velocity heterogeneities under the path, indicating that accurate phase velocity and attenuation analysis requires array measurements. T he Rayleigh waves are strongly affected by low velocity anomalies beneath Central California, with ray paths bending eastward as waves travel toward the south, while Love waves are less affected, providing observables to constrain the depth extent of anomalies. Strong lateral gradients in the lithospheric structure between the continent and the ocean are the likely cause of the path deflections.     

利用地震背景噪声提取西山村滑坡高频面波信号

滑坡体三维速度结构为滑坡治理、灾害防治、风险防范以及理解滑坡的动力学过程提供了关键信息,而高频面波成像是研究浅层速度结构的重要手段.本文详细介绍了利用2016年12月1日在四川理县西山村滑坡上布设的38个仪器记录的三分量连续地震噪声数据,提取1～5Hz的基阶Love波和Rayleigh波经验格林函数,分析了不同参数和处理方法对背景噪声互相关计算结果的影响.结果表明,1天左右的连续波形记录裁剪至1200s的时间长度,进行互相关叠加就可以得到较为稳定的经验格林函数.西山村滑坡体上Love波和Rayleigh波的群速度分别为约400m·s-1和700m·s-1,并且Love波信噪比高于Rayleigh波.此外,我们还利用聚束分析方法对噪声源的位置进行了分析,发现1～5Hz的背景噪声主要来自滑坡东南侧附近杂谷脑河水的搬运作用.这些高频面波数据和噪声源位置为获取滑坡浅层三维速度结构提供了重要输入,同时也为研究滑坡体速度结构随时间的变化提供了基础.     

长江中下游成矿带及邻区地壳剪切波速度结构和径向各向异性

收集了安徽、江西、浙江、江苏、湖北和河南6个省的区域地震台网138个宽频地震台站以及中国地质大学(北京)在长江中下游成矿带布设的19个流动宽频地震台站的三分量背景噪声数据,利用背景噪声面波层析成像方法,获得了长江中下游成矿带及其邻区地壳三维剪切波速度结构和径向各向异性特征.首先获得了5~38s周期的瑞利波和勒夫波相速度,结果显示短周期(16s)的瑞利波和勒夫波相速度与研究区内的主要地质构造单元具有良好的相关性,但在中长周期(20~30s)瑞利波相速度显示大别造山带东部为明显低速特征,而勒夫波相速度并未表现出异常特征.研究区域地壳三维有效剪切波速度和径向各向异性结果显示:苏北盆地和江汉盆地上地壳都表现为低速和正径向各向异性特征,华北克拉通东南部也表现为正径向各向异性,这可能与盆地浅部沉积层的水平层理结构相关.大别造山带中地壳显示为弱的正径向各向异性,同时其东部下地壳显示为低剪切波速度和强的正径向各向特征,可能是由于其在造山后发生了中下地壳的流变变形,引起各向异性矿物近水平排列所导致的.长江中下游成矿带内的鄂东南和安庆—贵池矿集区中地壳弱的负径向各向异性可能是由于深部岩浆向上渗透时所产生的有限应力导致结晶各向异性矿物的垂直排列所引起的.整个长江中下游成矿带下地壳都表现出正径向各向异性特征,可能是由于在伸展拉张的构造作用力下,下地壳矿物的晶格优势水平排列所引起的.     

2015年8月12日天津化学爆炸产生的多模式面波分析及其应用研究

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的浅地表速度梯度较大,符合典型的盆地结构特征.本文的研究结果为类似爆炸等突发事件快速定位提供了新的思路,有助于灾后救援的迅速展开;同时得到天津滨海新区及周边浅层精细的速度结构,对于地震灾害评估有很大帮助.     

地震环境噪声成像中的地形影响校正

以传统地震环境噪声面波成像方法研究地壳速度结构时,在一些极端的地形条件下,结果与真实结构会存在较大偏差.我们以地震波场三维正演模拟为基础,提出了一种地形校正方法.我们保留了传统噪声面波成像简单的两步反演法,在面波层析成像和一维速度结构反演的基础上,通过地震波场三维模拟近似估计地形和散射波场的影响,并据此校正瑞利波频散曲线,最终反演得到校正地形影响的S波速度结构.理论测试与在实际观测数据上的应用都证明了校正方法的有效性,同时也显示了地形校正的必要性.     

Soil characterization in urban areas of the Bajo Segura Basin (Southeast Spain) using H/V, F-K and ESAC methods

Bajo Segura Basin, placed in the south of Alicante province (southeast Spain), is a region of moderate seismicity, which presents frequent seismic episodes. This region, as the rest of south-eastern part of Spain, is growing both economically and demographically. Hence, the presence of soft soils in a seismic region where many towns and villages are placed shows the importance of an adequate characterization of the site effects study in this area.The present work investigates the use of noise measurements for soil characterization by the application of H/V method and array techniques: extended spatial autocorrelation (ESAC) and frequency-wave number (F-K) in three urban areas of the Bajo Segura Basin. The application of these methods will help us to estimate the resonance frequencies, as well as shear-wave velocity profiles, having a better understanding of the sediment depth and site effects phenomena at the studied places.The estimated velocity profiles are used to model the fundamental mode Rayleigh wave ellipticity curves and the S-wave site transfer functions. Finally, the correlation between this curves and the experimental H/V analysis has allowed us a better assessment of the site response in the studied areas.