Feasibility of borehole ambient noise interferometry for permanent reservoir monitoring

The analysis of seismic ambient noise acquired during temporary or permanent microseismic monitoring campaigns (e.g., improved/enhanced oil recovery monitoring, surveillance of induced seismicity) is potentially well suited for time‐lapse studies based on seismic interferometry. No additional data acquisition required, ambient noise processing can be automatized to a high degree, and seismic interferometry is very sensitive to small medium changes. Thus there is an opportunity for detection and monitoring of velocity variations in a reservoir at negligible additional cost and effort. Data and results are presented from an ambient noise interferometry study applied to two wells in a producing oil field in Romania. Borehole microseismic monitoring on three component geophones was performed for four weeks, concurrent with a water‐flooding phase for improved oil recovery from a reservoir in ca. 1 km depth. Both low‐frequency (2 Hz–50 Hz) P‐ and S‐waves propagating through the vertical borehole arrays were reconstructed from ambient noise by the virtual source method. The obtained interferograms clearly indicate an origin of the ambient seismic energy from above the arrays, thus suggesting surface activities as sources. It is shown that ambient noise from time periods as short as 30 seconds is sufficient to obtain robust interferograms. Sonic log data confirm that the vertical and horizontal components comprise first arrivals of P‐wave and S‐waves, respectively. The consistency and high quality of the interferograms throughout the entire observation period further indicate that the high‐frequency part (up to 100 Hz) represents the scattered wave field. The temporal variation of apparent velocities based on first‐arrival times partly correlates with the water injection rate and occurrence of microseismic events. It is concluded that borehole ambient noise interferometry in production settings is a potentially useful method for permanent reservoir monitoring due to its high sensitivity and robustness.     

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

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

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.     

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

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

Retrieving virtual reflection responses at drill‐bit positions using seismic interferometry with drill‐bit noise

In the field of seismic interferometry, researchers have retrieved surface waves and body waves by cross‐correlating recordings of uncorrelated noise sources to extract useful subsurface information. The retrieved wavefields in most applications are between receivers. When the positions of the noise sources are known, inter‐source interferometry can be applied to retrieve the wavefields between sources, thus turning sources into virtual receivers. Previous applications of this form of interferometry assume impulsive point sources or transient sources with similar signatures. We investigate the requirements of applying inter‐source seismic interferometry using non‐transient noise sources with known positions to retrieve reflection responses at those positions and show the results using synthetic drilling noise as source. We show that, if pilot signals (estimates of the drill‐bit signals) are not available, it is required that the drill‐bit signals are the same and that the phases of the virtual reflections at drill‐bit positions can be retrieved by deconvolution interferometry or by cross‐coherence interferometry. Further, for this case, classic interferometry by cross‐correlation can be used if the source power spectrum can be estimated. If pilot signals are available, virtual reflection responses can be obtained by first using standard seismic‐while‐drilling processing techniques such as pilot cross‐correlation and pilot deconvolution to remove the drill‐bit signatures in the data and then applying cross‐correlation interferometry. Therefore, provided that pilot signals are reliable, drill‐bit data can be redatumed from surface to borehole depths using this inter‐source interferometry approach without any velocity information of the medium, and we show that a well‐positioned image below the borehole can be obtained using interferometrically redatumed reflection responses with just a simple velocity model. We discuss some of the practical hurdles that restrict the application of the proposed method offshore.     

3-D surface wave group velocity distribution in Central-Southern Africa

Group velocities estimated from fundamental mode Love and Rayleigh waves are used in a tomography process in central-southern Africa. The waves were generated by eighteen earthquakes, which occurred along the East African Rift and recorded at BOSA, LBTB and SLR seismic stations in southern Africa. The group velocities from Love and Rayleigh waves were isolated using the Multiple Filter Technique (MFT) at the period range of 10 to 50 seconds. The tomography method developed by Ditmar and Yanovskaya (1987) and Yanovskaya and Ditmar (1990), was applied to calculate the lateral distribution of surface wave group velocities in central-southern Africa. The results of the tomographic inversion were plotted as distribution maps. In addition to the maps, I also produced two velocity cross-sections across the area of study. The velocity distribution maps show the regional tectonic units, though with poor resolution. The azimuthal bias of the surface wave paths is reflected in the distribution of the group velocities. The Moho depth appears to correlate with velocities at a period of about 30 s. A low velocity feature observed beneath the Zimbabwe craton implies a thickening upper asthenosphere and lithospheric thinning beneath the Zimbabwe craton. Also estimated was a shear wave velocity model beneath the Zimbabwe craton.     

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.     

Estimation of rocking and torsion associated with surface waves extracted from recorded motions

By exploiting the capability of identifying and extracting surface waves existing in a seismic signal, we can proceed to estimate the angular displacement (rotation about the horizontal axis normal to the direction of propagation of the wave; rocking) associated with Rayleigh waves as well as the angular displacement (rotation about the vertical axis; torsion) associated with Love waves.For a harmonic Rayleigh (Love) wave, rocking (torsion) would be proportional to the harmonic vertical (transverse horizontal) velocity component and inversely proportional to the phase velocity corresponding to the particular frequency of the harmonic wave (a fact that was originally exploited by Newmark (1969) [15] to estimate torsional excitation). Evidently, a reliable estimate of the phase velocity (as a function of frequency) is necessary. As pointed out by Stockwell (2007) [17], because of its absolutely referenced phase information, the S-Transform can be employed in a cross-spectrum analysis in a local manner. Following this suggestion a very reliable estimate of the phase velocity may be obtained from the recordings at two nearby stations, after the dispersed waves have been identified and extracted. Synthesis of the abovementioned harmonic components can provide a reliable estimate of the rocking (torsional) motion induced by an (extracted) Rayleigh (Love) wave.We apply the proposed angular displacement estimation procedure for two well recorded data sets: (1) the strong motion data generated by an aftershock of the 1999 Chi-Chi, Taiwan earthquake and recorded over the Western Coastal Plain (WCP) of Taiwan, and (2) the strong motion data generated by the 2010 Darfield, New Zealand earthquake and recorded over the Canterbury basin. The former data set is dominated by basin-induced Rayleigh waves while the latter contains primarily Love waves.     

Lunar noise correlation,imaging and monitoring

Passive seismic techniques have revolutionarised seismology,leading for example to increased resolution in surface wave tomography,to the possibility to monitor changes in the propagation medium,and to many new processing strategies in seismic exploration.Here we review applications of the new techniques to a very particular dataset,namely data from the Apollo 17 lunar network.The special conditions of the lunar noise environment are investigated,illustrating the interplay between the properties of the noise and the ability to reconstruct Green's functions.With a dispersion analysis of reconstructed Rayleigh waves new information about the shallow shear velocity structure of the Moon are obtained.Passive image interferometry is used to study the effect of temperature changes in the subsurface on the seismic velocities providing direct observation of a dynamic process in the lunar environment.These applications highlight the potential of passive techniques for terrestrial and planetary seismology.     

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

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

高频面波方法的若干新进展

面波多道分析方法(MASW)通过分析高频瑞雷波确定浅地表剪切波速度.在过去的20年中,由于该方法具有非侵入性、无损、高效及价格低的特点,越来越受到浅地表地球物理和地质工程学界的重视,视为未来最有希望的技术之一.这篇综述论文将介绍中国地质大学(武汉)浅地表地球物理团队近年来在研究高频面波的传播理论和应用中取得的部分成果.非几何波是一种仅存在于浅地表介质,尤其是未固结的沉积物中的独特的地震波.它的存在对快速而准确地获得表层S波速度有一定价值.我们的研究表明非几何波是一种具有频散特性的泄漏波.泄漏波的存在可能导致将其误认为瑞雷波的基阶或高阶能量,从而造成模式误判.这种模式误判会导致错误的反演结果.我们通过求取高基阶分离后的瑞雷波格林函数证明虚震源法瑞雷波勘探的可行性.这个结果将极大地降低野外瑞雷波勘探成本.勒夫波多道分析方法(MALW)中未知参数比瑞雷波的少,这使得勒夫波的频散曲线比瑞雷波的简单.因此,勒夫波反演更稳定,非唯一性更低.勒夫波数据生成的能量图像通常比瑞雷波的清晰,并具有更高的分辨率,从而可以更容易地拾取精确的勒夫波的相速度.利用雅克比矩阵分析波长与探测深度的关系表明对相同波长的基阶模式而言,瑞雷波的探测深度是勒夫波的1.3~1.4倍;而两种波的相同波长的高阶模式波的探测深度相同.我们也尝试了时间域勒夫波反演.按照勒夫波分辨率将地球模型剖分成了不同尺寸的块体,利用反卷积消除了地震子波对勒夫波波形的影响,通过更新每个块体的S波速度来拟合勒夫波波形,从而获得地下S波速度模型.该方法不基于水平层状模型假设,适用于任意二维介质模型.     

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

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

Crustal Structure of the Korean Peninsula Using Surface Wave Dispersion and Numerical Modeling

Surface wave dispersion is studied to obtain the 1-D average velocity structure of the crust in the Korean Peninsula by inverting group- and phase-velocities jointly. Group velocities of short-period Rayleigh and Love waves are obtained from cross-correlations of seismic noise. Multiple-filter analysis is used to extract the group velocities at periods between 0.5 and 20 s. Phase velocities of Rayleigh waves in 10- and 50-s periods are obtained by applying the two-station method to teleseismic data. Dispersion curves of all group and phase velocities are jointly inverted for the 1-D average model of the Korean Peninsula. The resultant model from surface wave analysis can be used as an initial model for numerical modeling of observations of North Korean events for a velocity model appropriated to the Korean Peninsula. The iterative process is focused especially on the surface sedimentary layer in the numerical modeling. The final model, modified by numerical modeling from the initial model, indicates that the crust shear wave velocity increases with depth from 2.16 km/s for a 2-km-thick surface sedimentary layer to 3.79 km/s at a Moho depth of 33 km, and the upper mantle has a velocity of 4.70 km/s.     

Seismic crust structure beneath the Aegean region in southwest Turkey from radial anisotropic inversion of Rayleigh and Love surface waves

The Turkish plate is covered by hundreds of accelerometer and broadband seismic stations with less than 50 km inter-station distance providing high-quality earthquake recordings within the last decade. We utilize part of these stations to extract the fundamental mode Rayleigh and Love surface wave phase and group velocity data in the period range 5–20 s to determine the crust structure beneath the Aegean region in southwest Turkey. The observed surface wave signals are interpreted using both single-station and two-station techniques. A tomographic inversion technique is employed to obtain the two-dimensional group velocity maps from the single-station group velocities. One-dimensional velocity–depth profiles under each two-dimensional mesh point, which are jointly interpreted to acquire the three-dimensional image of the shear-wave velocities underneath the study area, are attained by utilizing the least-squares inversion technique, which is repeated for both Rayleigh and Love surface waves. The isotropic crust structure cannot jointly invert the observed Rayleigh and Love surface waves where the radial anisotropic crust better describes the observed surface wave data. The intrusive magmatic activity related to the northward subducting African plate under the Turkish plate results the crust structure deformations, which we think, causing the observed radial anisotropy throughout complex pattern of dykes and sills. The magma flow resulting in the mineral alignment within dykes and sills contributes to the observed anisotropy. Due to the existence of dykes, the radial anisotropy in the upper crust is generally negative, i.e., vertically polarized S-waves (Vsv) are faster than horizontally polarized S-waves (Vsh). Due to the existence of sills, the radial anisotropy in the middle-to-lower crust is generally positive, i.e., horizontally polarized S-waves (Vsh) are faster than vertically polarized S-waves (Vsv). Similar radial anisotropic results to those of the single-station analyses are obtained by the two-station analyses utilizing the cross-correlograms. The widespread volcanic and plutonic rocks in the region are consistent with the current seismic interpretations of the crustal deformations.     

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

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

Shallow seismicity, triggered seismicity, and ambient noise tomography at the long-dormant Uturuncu Volcano, Bolivia

Using a network of 15 seismometers around the inflating Uturuncu Volcano from April 2009 to 2010, we find an average rate of about three local volcano-tectonic earthquakes per day, and swarms of 5–60 events a few times per month with local magnitudes ranging from −1.2 to 3.7. The earthquake depths are near sea level, more than 10 km above the geodetically inferred inflation source and the Altiplano Puna Magma Body. The Mw 8.8 Maule earthquake on 27 February 2010 triggered hundreds of earthquakes at Uturuncu with the onset of the Love and Rayleigh waves and again with the passage of the X2/X3 overtone phases of Rayleigh waves. This is one of the first incidences in which triggering has been observed from multiple surface wave trains. The earthquakes are oriented NW–SE similar to the regional faults and lineaments. The b value of the catalog is 0.49, consistent with a tectonic origin of the earthquakes. We perform ambient noise tomography using Love wave cross-correlations to image a low-velocity zone at 1.9 to 3.9 km depth below the surface centered slightly north of the summit. The low velocities are perhaps related to the hydrothermal system and the low-velocity zone is spatially correlated with earthquake locations. The earthquake rate appears to vary with time—a seismic deployment from 1996 to 1997 reveals 1–5 earthquakes per day, whereas 60 events/day were seen during 5 days using one seismometer in 2003. However, differences in analysis methods and magnitudes of completeness do not allow direct comparison of these seismicity rates. The rate of seismic activity at Uturuncu is higher than at other well-monitored inflating volcanoes during periods of repose. The frequent swarms and triggered earthquakes suggest the hydrothermal system is metastable.     

青藏高原中部Quasi-Love波的识别及其转换点揭示的东西向方位各向异性变化

目前人们利用4种基本的地震波现象研究地震各向异性,如横波双折射、面波散射、与传播方向有关的走时异常和PS转换波震相.本文利用面波散射产生的Quasi-Love（QL）波研究青藏高原上地幔顶部的各向异性结构特征.首先利用中国地震台网昌都（CAD）台记录的地震波形资料识别出产生QL波的路径,并利用合成地震记录和垂直偏振极性分析证实所观测到的为QL波,而不是高阶振型的Rayleigh波或其他体波震相;然后由Rayleigh波、Love波和QL波的群速度估算了各向异性结构横向变化的转换点;不同周期时,转换点的位置不同,这种频率依赖性还需要进一步的模拟研究.Love波向Rayleigh波耦合（产生QL波）的转换点位置揭示了青藏高原面波方位各向异性变化特征,并以南北向构造带的东西分段性、上地幔流引起的地球内力诱导岩石形变解释了青藏高原各向异性的东西向差异性.     

Numerical simulations of full-wave fields and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling twodimensional wave fields and therefore cannot accurately simulate three-dimensional (3-D) full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity–stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity V_x, V_y, and V_z for the same node in 3-D staggered-grid finite difference models by calculating the average value of V_y, and V_z of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.     

Joint inversion of teleseismic P waveforms and surface-wave group velocities from ambient seismic noise in the Bohemian Massif

Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle. We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly felsic.     

Surface Wave Group Velocity Tomography Imaging from Ambient Noise for Fujian Province and Its Adjacent Areas

Two-month continuous waveforms of 108 broadband seismic stations in Fujian Province and its adjacent areas are used to compute noise cross-correlation function (NCF). The signal quality of NCF is improved via the application of time-frequency phase weighted stacking. The Rayleigh and Love waves group velocities between 1s-20s are measured on the symmetrical component of the NCF with the multiple filter method. More than 5,000 Rayleigh wave dispersion curves and about 4,000 Love wave dispersion curves are obtained and used to invert for group velocity maps. This data set provides about 50km resolution that is demonstrated with checkerboard tests. Considering the off great circle effect in inhomogeneous medium, the ray path is traced based on the travel time field computed with a finite difference method. The inverted group velocity maps show good correlation with the geological features in the upper and middle crust. The Fuzhou basin and Zhangzhou basin showed low velocity on the short period group velocity maps. On the long period group velocity maps, the low velocity anomaly in the high heat flow region near Zhangzhou and clear velocity contrast across the Zhenghe-Dapu faults, which suggests that the Zhenghe-Dapu fault might be a deep fault.