Experiments on Exploration of Shallow Fine Structures and the Construction of the 1-D Velocity Model in the Pingtan Island, Fujian

112 short-period seismographs were set up in the 400km² area of Pingtan Island and its surrounding areas in Fujian. The combined observations of the airgun source and ambient noise source were carried out using a dense array to receive the 387 airgun signals excited around the island and one month of continuous ambient noise recording. The 1-D P-wave and S-wave shallow velocity model of Pingtan Island is obtained by the inversion of the airgun body wave''s first arrival time data, and the reliability of the velocity model is verified by using the surface wave phase velocity dispersion curve, which can provide initial model for subsequent 3-D imaging. The experimental results show that this experiment is a successful demonstration of local scale green non-destructive detection, which can provide basic data for shallow surface structure research and strong vibration simulation of the Pingtan Island.     

3D S-wave velocity structure of the Ningdu basin in Jiangxi province inferred from ambient noise tomography with dense array

The Ningdu basin, located in southern Jiangxi province of southwest China, is one of the Mesozoic basin groups which has exploration prospects for geothermal energy. A study on the detailed velocity structure of the Ningdu basin can provide important information for geothermal resource exploration. In this study, we deployed a dense seismic array in the Ningdu basin to investigate the 3D velocity structure and discuss implications for geothermal exploration and geological evolution. Based on the dense seismic array including 35 short-period (5 s-100 ?Hz) seismometers with an average interstation distance of ～5 ?km, Rayleigh surface wave dispersion curves were extracted from the continuous ambient noise data for surface wave tomographic inversion. Group velocity tomography was conducted and the 3D S-wave velocity structure was inverted by the neighborhood algorithm. The results revealed obvious low-velocity anomalies in the center of the basin, consistent with the low-velocity Cretaceous sedimentary rocks. The basement and basin-controlling fault can also be depicted by the S-wave velocity anomalies. The obvious seismic interface is about 2 ?km depth in the basin center and decreases to 700 ?m depth near the basin boundary, suggesting spatial thickness variations of the Cretaceous sediment. The fault features of the S-wave velocity profile coincide with the geological cognition of the western boundary basin-controlling fault, which may provide possible upwelling channels for geothermal fluid. This study suggests that seismic tomography with a dense array is an effective method and can play an important role in the detailed investigations of sedimentary basins.     

A new stacking technique in ambient noise tomography

We have developed a new stacking technique in ambient noise tomography to obtain high-quality dispersion curves of Rayleigh waves.This technique is used to stack the vertical components of the Estimated Green Functions(EGFs) obtained respectively from cross correlation of the ambient noise data recorded by a remote seismic station and one of the short distance seismic stations of a seismic array.It is based on a phase-matched filter and is implemented by a four-step iterative process:signal compression,stacking,signal extraction and signal decompression.The iterative process ends and gives the dispersion curve of Rayleigh wave when the predicted one and the processing result converge.We have tested the method using the vertical components of synthetic Rayleigh wave records.Results show that this new stacking method is stable and it can improve the quality of dispersion curves.In addition,we have applied this method to real data.We see that the results given by our new technique are obviously better than the ones employing the traditional method which is a three-step process:signal compression,signal extraction and signal decompression.In conclusion,the new method proposed in this paper can improve the signal to noise ratio of EGFs,and can therefore potentially improve the resolution of ambient noise tomography.     

龙门山断裂带中上地壳速度结构的短周期环境噪声成像

本文利用在龙门山断裂带周边布设的57个台站自2008年11月至2009年11月为期一年的垂直分量连续地震环境噪声数据,通过短周期地震环境噪声成像方法,获得了龙门山断裂带中北段地壳25km深度范围的S波精细速度结构.结果表明:(1)龙门山断裂带周边区域10km以上的速度结构与地表断裂的分布形态具有良好的一致性,速度结构控制了龙门山主要断层的深部延展特征;在15km及以下深度,S波速度结构呈现沿龙门山和沿岷山隆起走向的交叉构造格局,由此造成的速度结构差异可能影响了汶川地震的破裂过程;(2)速度结构随深度的分布特征为龙门山断裂带主要断层的深部延伸形态给出了良好的约束,结果进一步确认了龙门山断裂中段的高角度铲型断裂构造特征;(3)研究区的南端发现了龙门山断裂下方20km以下深度具有与松潘地块中地壳低速层相关的低速结构的迹象,这可能是汶川地震破裂带南段22km左右深度存在脆韧转换带的一个证据.研究结果显示出密集台阵和短周期环境噪声成像方法在地壳浅部精细结构和断层探测研究中具有巨大潜力.     

地脉动空间自相关方法反演 浅层S波速度结构

在近地表地球物理领域, 基于地脉动(或称背景噪声)提取的面波频散曲线反演地下S波速度结构是一种简单经济的工程勘察方法. 本文基于地脉动的空间自相关方法对一个微型台阵观测的背景噪声记录进行处理, 介绍了一种简单易行的提取频散曲线的数据处理方法, 获得了6.7—23 Hz频段的可靠频散曲线; 通过对该观测频散曲线与预测模型的频散曲线进行拟合, 反演得到S波速度结构. 结果表明, 该速度结构与钻孔直接测试的结果相吻合.      

新疆和田台阵PSD与PDF分析

和田台阵是我国第一个自主建成并运行至今的小孔径台阵,承担着监测印巴地区核试验以及我国西部地震活动的重要使命.台阵波形数据中充斥着背景噪声,直接影响着数据质量.为了评估台阵噪声水平,本文利用Welch平均周期法对9个子台记录的数十万条噪声样本进行功率谱估计,对出现的谱异常进行了总结归纳,通过绘制概率密度函数图以及单频曲线来研究背景噪声变化范围和规律,最后针对台阵降噪提出了建设性意见.研究结果表明,中心台长周期噪声功率谱密度随季节变化显著,具有周期性;受温度和气压影响,水平分量长周期噪声变化幅度较大,局部频段超出新高噪声模型,建议改善仪器安装条件,或者利用数学方法进行校正.所有子台短周期噪声变化规律与长周期相反,受到采石场影响,谱密度曲线在4~8 Hz之间出现形态规则的高频尖刺,A1、B3、B4子台最为明显,可以通过窄带滤波或者聚束予以压制.本文取得的研究成果为台阵运维提供重要依据,除此之外,总结出的不同地震频谱特征也为地震解释工作提供重要参考.     

Crustal S-wave velocity structure of the Yellowstone region using a seismic ambient noise method

The Yellowstone volcano is one of the largest active volcanoes in the world, and its potential hazards demand detailed seismological and geodetic studies. Previous studies with travel time tomography and receiver functions have revealed a low-velocity layer in the crust beneath the Yellowstone volcano, suggesting the presence of a magma chamber at depth. We use ambient seismic noise from regional seismic stations to retrieve short-period surface waves and then study the shallow shear velocity structure of the Yellowstone region by surface wave dispersion analysis. We first obtained a crustal model of the area outside of the Yellowstone volcano and then constructed an absolute shear wave velocity structure in combination with receiver function results for the crust beneath the Yellowstone volcano. The velocity model shows a low-velocity layer with shear velocity at around 1.3 km/s, suggesting that a large-scale magma chamber exists at shallow levels within the crust of the Yellowstone volcanic region.     

分布式光纤声波传感系统在近地表成像中的应用Ⅱ:背景噪声成像

随着城市化的发展,城市地球物理日益成为地球物理研究的重要方向,地震成像是构建城市地下空间三维/四维图像的重要手段,但面临观测成本高的困难.近年来国际上新发展的分布式光纤声波传感器作为高密度地震观测系统已经在地震层析成像方面得到了应用,在提高成像分辨率的同时,又降低了观测成本.本研究使用国产分布式光纤声波传感器开展了观测实验,利用480m埋地光缆记录了13h背景噪声,计算得到噪声互相关函数,获得了高频Rayleigh面波信号.采用多道面波分析方法提取相速度频散曲线,其结果与传统检波器记录和主动源结果较为一致.采用遗传算法反演得到了研究区内二维S波速度剖面,获得了下方沉积物横向变化特征.通过本次实验,初步验证了国产设备开展地震背景噪声成像研究、构建地下浅层结构模型的可行性.     

四川盐源盆地短周期密集台阵背景噪声分布特征分析

短周期密集台阵的高频背景噪声互相关函数(NCF)是探查地球浅层精细结构的重要数据.然而高频背景噪声成分复杂且容易分布不均,分析其对NCF信号提取的影响,有助于获取可靠成像结果.本文基于布设于川滇地区盐源盆地的209个短周期台站组成的盐源台阵,利用密集台阵的噪声水平评估以及基于NCF的相干噪声分析两种方法,分析了其记录到的噪声波场特征及其对NCF的影响.结果表明,盐源台阵的整体噪声水平呈现北低南高的不均匀分布,高频噪声水平的强弱受控于当地的人类活动,亦受到浅部松散沉积层的影响.台阵垂直分量NCF中主要信号为基阶Rayleigh波,且产生该信号的相干噪声源的优势方位在不同频带具有较大区别:0.3～0.5Hz的噪声源强度较强且随时间变化较为稳定,主要能量来自台阵的南侧;0.5～1Hz的相干噪声源强度较低,有两个优势方向,其中较强的一个来自于台阵南侧,可能与0.3～0.5 Hz的噪声同源,较弱的一个来自于台阵北偏东方向;1～1.5Hz的背景噪声有四个较弱的优势方向,在台阵的不同区域有不同的优势方向,可能受到不同的局部噪声源的控制.垂向NCF中Rayleigh波的信噪比主要受控于波场的复杂程度,台阵南部受人文活动及沉积层影响,噪声水平较高,且由于盆山边缘复杂的反射、散射作用,其NCF波形复杂,信噪比偏低.受高频噪声源分布不均与及复杂地质结构的共同影响,盐源台阵的高频NCF中的信号复杂,后续对面波频散特征的提取应充分考虑噪声源对NCF的影响以获取可靠结果.     

利用噪声记录估计福建地区中上地壳体波速度结构

In this paper,the dispersion curves of the Rayleigh wave and Love wave were extracted from the seismic noise records of 25 broadband stations of the Fujian Seismic Network, and inverted for the lithosphere velocity structure. Furthermore,the velocity model was verified by the seismic explosion observations. Our results indicate that the resolution of the lithosphere velocity structure obtained by this method is good in the shallow part,but in the deep part,inversion accuracy for the wave velocity structure is low,which is caused mainly by the small inter-station distance chosen in the paper. Thus the wave dispersion curves have high accuracy in the short-period part,but the warp of the wave dispersion curve in long-period part is large. Considering the results from both the noise inversion and the traditional inversion,we finally present a new velocity model,and the theoretical travel time calculated with the new model matches the explosion travel time very well.     

利用噪声记录估计福建地区中上地壳体波速度结构

本文利用福建地震台网25个宽频带台站的噪声记录提取瑞利波和勒夫波频散曲线反演地壳波速结构,并利用爆破观测结果对该速度模型进行检验.结果表明:利用该方法得到的地壳波速结构在浅部分辨率较好,深部波速结构的反演结果精度较低,这主要是由于本文所选台站间距较小,提取得到的面波频散曲线在短周期部分精度较高,而长周期部分面波频散曲线的偏差较大.综合噪声反演结果与传统反演结果,本文最终给出了一个新的速度模型,利用该模型计算得到的理论走时与爆破观测走时具有很好的一致性.     

天然源面波勘探台阵对比试验

为了对比天然源面波勘探不同台阵布局的探测效果, 筛选出探测成果可靠、 效率高和便于野外施工的天然源面波勘探台阵阵形, 在天水市黄土覆盖区的同一场地分别用4种常见的阵形进行数据采集试验, 并对各种阵形数据使用空间自相关法或扩展空间自相关法提取相应的频散曲线, 通过反演得到了试验点地下的浅层速度结构模型. 分析对比试验结果表明: 4种台阵提取的频散曲线数值很相近; 频散谱能量集中度较高的是嵌套式等边三角形和圆形台阵, L形和直线形台阵相对分散; L形台阵低频段(4—8 Hz)比直线形台阵差, 其高频段(8—40 Hz)比直线形台阵好. 针对直线形台阵在高频段信噪比较低的情况, 在确保探测成果可靠性的前提下, 为了提高探测效率, 提出了在同一直线形台阵开展天然源与人工源面波联合勘探的数据采集方法. 实验结果证实, 这种联合勘探方法不仅可弥补直线形台阵高频段的不足, 确保探测精度和结果的可靠性, 而且还能实现“高低”频兼顾, 即“深浅”兼顾.      

Array performances for ambient vibrations on a shallow structure and consequences over <Emphasis Type="Italic">V</Emphasis><Subscript><Emphasis Type="Italic">s</Emphasis></Subscript> inversion

Valuable information about one-dimensional soil structures can be obtained by recording ambient vibrations at the surface, in which the energy contribution of surface waves predominates over the one of other types of waves. The dispersion characteristics of surface waves allow the retrieval of the shear-wave velocity as a function of depth. Microtremor studies are usually divided in two stages: deriving the dispersion (or auto-correlation) curve from the recorded signals and inverting it to obtain the site velocity profile. The possibility to determine the dispersion curve over the adequate frequency range at one site depends on the array aperture and on the wavefield spectra amplitude that can be altered by filtering effects due to the ground structure. Microtremors are usually recorded with several arrays of various apertures to get the spectral curves over a wide frequency band, and different methods also exist for processing the raw signals. With the objective of defining a strategy to achieve reliable results for microtremor on a shallow structure, we analyse synthetic ambient vibrations (vertical component) simulated with 333 broadband sources for a 25-m deep soil layer overlying a bedrock. The first part of our study is focused on the determination of the reliable frequency range of the spectral curves (dispersion or auto-correlation) for a given array geometry. We find that the wavenumber limits deduced from the theoretical array re sponse are good estimates of the valid spectral curve range. In the second part, the spectral curves are calculated with the three most popular noise-processing techniques (frequency–wavenumber, high-resolution frequency–wavenumber and spa tial auto-correlation methods) and inverted indi vidually in each case. The inversions are performed with a tool based on the neighbour hood algorithm that offers a better estimation of the global uncertainties than classical linearised methods, especially if the solution is not unique. Several array apertures are necessary to construct the dispersion (auto-correlation) curves in the appropriate frequency range. Considering the final velocity profiles, the three tested methods are almost equivalent, and no significant advantage was found for one particular method. With the chosen model, all methods exhibit a penetration limited to the bedrock depth, as a consequence of the filtering effect of the ground structure on the vertical component, which was observed in numerous shallow sites.     

1D velocity structure of the Po River plain (Northern Italy) assessed by combining strong motion and ambient noise data

Strong ground motions recorded on the sedimentary deposits of the Po River alluvial plain during the Emilia (Northern Italy) Mw 5.7 earthquake of May 29, 2012 are used to assess the vertical profile of shear-wave velocity above the limestone basement. Data were collected by a linear array installed for site effect studies after the Mw 5.9 mainshock of May 20, 2012. The array stations, equipped with both strong and weak motion sensors, are aligned in the South–North direction, at distances ranging from 1 to 26 km from the epicenter. The vertical components of ground motion show very distinctive, large-amplitude, low-frequency dispersive wave trains. Wavelet analysis yields group-velocity dispersion curve in the 0.2–0.7 Hz frequency band. The availability of a long ambient noise record allows estimates of the site resonance frequency along with its stability among stations. The joint inversion of dispersion of surface waves and ellipticity curves derived from ambient noise H/V allows extending investigations down to the sediment-limestone interface, at a depth of about 5,000 m. Our results add new information about the velocity structure at a scale that is intermediate between the local scale already investigated by other authors with small-aperture arrays using ambient noise and the regional scale inferred from modeling of seismogram waveforms recorded at hundreds of kilometers from the source.     

Seismic noise array measurements using broadband stations and vertical geophones: preliminary outcomes for the suitability on f-k analysis

In recent years, seismic noise array measurements have emerged as a promising way for estimating Vs profiles and characterizing the corresponding soil structure. Although methods based on analysis of seismic noise are simpler, cheaper and faster than conventional soil penetrating techniques, array deployment requires the use of several stations which are not always available for all the research groups. Broadband stations are expensive, but even short-period sensors can require a significant economical effort as it is necessary to use several. In this paper, the results from broadband stations and vertical geophones using the f-k method have been compared. Circular arrays with both types of instruments were installed north and south of the Himalayan Frontal Thrust (HFT) also known as Mohand Thrust, with different soil characteristics. The comparison carried out in terms of dispersion curves shows an excellent agreement between broadband stations and geophones recordings in the frequency range of applicability of the arrays used. These results indicate the possibility of using seismic refraction equipments for the f-k analysis in the frequency band of geotechnical engineering interest. This has potentially great importance since the same equipment is indicated to be applicable for both seismic noise measurements as well as for active reflection seismic surveys.     

基于频率-贝塞尔变换法的关东盆地S波速度成像

将一种新的方法——频率-贝塞尔变换法（F-Jmethod）应用于日本NIED在关东盆地布设的MeSO-net台网的背景噪声数据中,证明频率-贝塞尔变换法可以有效地从背景噪声中提取高阶频散曲线.利用提取的基阶和高阶频散曲线反演关东盆地区域的沉积层和地震基岩层的S波速度结构,并将我们反演得到的S波速度结构与Koketsu等提出的日本综合速度结构模型进行对比讨论.我们的例子证明,在基阶面波的基础上,高阶面波能减少在反演中的非唯一性,得到更为准确的S波速度结构.     

Effect of earthquakes on ambient noise cross-correlation function

Surface wave tomography method based on analysis of ambient noise is widely used during the last decade. It is assumed that correlated component of noise is composed of surface waves generated by sources distributed over the Earth’s surface more or less uniformly. In such a case the cross-correlation function (CCF) at two stations may be considered as the Green’s function of surface wave. This function should be symmetric relatively to zero time. However analysis of CCF at the stations located at the East-European Platform shows that as a rule CCF is characterized with a strong asymmetry. Since “purered noise cannot be extracted from seismic records due to superposition of earthquake signals, the method for calculation of CCF includes amplitude normalization for suppression of earthquakes that reduces signals from earthquakes to a noise level. The parts of records containing waves from earthquakes are neglected because of their short duration. Present study shows that this contribution turns out to be dominant at periods larger than 20–40 s. In other words, what is assumed as a “noisered in reality is a superposition of signals from earthquakes. This fact results in distortion of the Green’s function and of surface wave dispersion curve used in surface wave tomography if in the time interval used for calculation of CCF many earthquakes occur within a small area apart of an extension of the interstation path (clustering). Numerical modeling shows how clusters of sources affect CCF and dispersion curve correspondingly. Means for reducing this effect are outlined.     

利用背景噪声数据提取地震台站间 面波的可靠性分析 以中国大陆中东部地区的宽频带台站为例

利用中国大陆中东部地区以国家台网为主的100个分布均匀的宽频带地震台记录到的21个月的连续波形数据, 经过单台数据处理和互相关叠加计算后, 由时频分析法提取了研究区各台站对间瑞雷波的格林函数. 为了检验经验格林函数的可靠性和稳定性, 对沿部分路径的经验格林函数和频散曲线进行了质量评估. 检测结果表明, 自21个月叠加的台站对间背景噪声中提取的经验格林函数与实际的地震面波一致, 提取的格林函数可靠. 此外, 统计了使用从3—21个月不同长度数据叠加后, 经验格林函数信噪比大于10的频散曲线数目. 结果表明, 至少要使用12个月的数据才能提取到信噪比足够大, 数目足够多, 可用于反演面波速度结构的经验格林函数; 12个月的叠加时长, 可以保证30 s以下周期的频散曲线在时间上稳定.      

基于贝叶斯理论的接收函数与环境噪声联合反演

基于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.     

利用噪声研究地下介质速度结构及其变化在中国大陆的应用

背景噪声层析成像技术已被广泛应用于地壳和上地幔速度结构的研究,这种方法不依靠地震的发生和人工源爆破,只需记录连续的噪声信号而无需产生信号,因为噪声穿过地下介质时会携带信息,然后通过利用台站记录到的连续背景噪声数据进行互相关计算和叠加,即可得到台站间的经验格林函数,从而获取对地下结构的认识。这种方法已经很好地应用于中国的东北地区、华北克拉通、青藏高原以及华南地区,并成功地揭示了这些地区地壳与上地幔顶部的速度结构。此外近年来,一些学者还利用噪声互相关技术研究地下介质地震波速度随时间的变化,通过对汶川大地震前后连续噪声记录的研究发现,大震发生后呈现同震波速降低和震后波速逐渐恢复的特点,这表明可以通过观测地震波特性的变化来监测地下应力的变化,从而为大震的预测预防工作提供科学依据。本文主要综述了近些年来背景噪声技术及其在中国大陆地区的应用。