用能量累积法检测地震波雷达信号

将具有高度重复性的地震波雷达长时间向地壳内发射线性调频信号,经过地下介质的传播后到达地面,用地震仪器在监测点和检测点记录下来,通过分析数据来了解地壳速度结构及波速变化.线性调频信号是一种非平稳信号,它的频率随时间线性变化,有很好的能量聚集性,非常适合做时间-频率分析.本文用短时傅里叶变换对监测点的信号进行时间-频率分析,以检验地震波雷达发射信号的时间和频率是否和控制系统一致.通过WignerVille分布将地震波雷达发射的信号能量聚集在线性调频直线上,再用Hough变换累积聚集的能量形成波峰,按照线性调频直线的倾角提取波峰所在行,计算到时后构成地震波走时曲线图.用靠近本次实验地点的H-21剖面得到的地壳速度结构正演该测线的Pg、Sg、PmP和SmS的折合走时曲线,并与用能量累积法提取出的地震波走时曲线进行对比,分析结果表明:地震波雷达发射线性调频信号的时间和频率都符合控制要求,重复性高达99.9%以上,可以清晰地分辨出Pg、Sg震相,并且PmP和SmS震相可辨.     

基于时变窄带滤波技术提取可控震源扫频信号方法研究

精密主动地震监测为我们主动探测地下介质结构, 并监视其动态变化提供了一条可能的技术途径。 由于精密控制震源释放的能量强度小, 随着传播距离的增加, 信号的快速衰减, 在离震源较远处有用的震源信号被掩盖在很强的噪声中, 这对于震相的识别与走时的拾取精度有着很大的影响。 本文设计了一种时变窄带滤波器进行更为精细的滤波, 期望进一步提高观测数据信噪比, 再结合匹配滤波方法实现主动震源信号的检测与波形变换。 另外, 由于不同震相有一定的到时差, 使用时变窄带滤波器提取一个震相波形信息时, 将压制其他震相的波形信号, 从而实现了震相分离的技术。 仿真计算与实际资料处理显示了该方法在提高观测资料质量、 震相识别分辨率及震相分离方面具有一定的有效性与优越性。 通过对广东省新丰江库区精密可控震源试验数据进行时变窄带滤波方法的处理, 在震中距为200 km左右处的台站记录中检测到主动源信号, 体现了对远台弱信号的提取能力。     

Joint inversion of travel times and waveforms by means of annealing algorithms. Application to a seismic refraction experiment

An algorithm of annealing is applied to a joint inversion of travel times and waveforms belonging to a synthetically generated seismic refraction experiment. The medium (crust and upper mantle) is modelled by a set of plane stratified layers and a halfspace. The obtained structure (elastic parameters and depth of layers) shows that, in spite of contaminating seismic noise and poor knowledge of the seismic source, annealing methods are a good tool in these kinds of inversion problems. We think that many characteristics of the annealing process described here could be used with real data and more sophisticated media for the crust and upper mantle of the earth than the present example.     

Constraining the anisotropy structure of the crust by joint inversion of seismic reflection travel times and wave polarizations

In the context of wide-angle seismic profiling, the determination of the physical properties of the Earth crust, such as the elastic layer depth and seismic velocity, is often performed by inversion of P- and/or S-phases propagation data supplying the geometry of the medium (reflector depths) or any other structural parameter (P- or S-wave velocity, density...). Moreover, the inversion for velocity structure and interfaces is commonly performed using only seismic reflection travel times and/or crustal phase amplitudes in isotropic media. But it is very important to utilize more available information to constrain the non-uniqueness of the solution. In this paper, we present a simultaneous inversion method of seismic reflection travel times and polarizations data of transient elastic waves in stratified media to reconstruct not only layer depth and vertical P-wave velocity but also the anisotropy feature of the crust based on the estimation of the Thomsen’s parameters. We carry out a checking with synthetic data, comparing the inversion results obtained by anisotropic travel-time inversion to the results derived by joint inversion of seismic reflection travel times and polarizations data. The comparison proves that the first procedure leads to biased anisotropic models, while the second one fits nearly the real model. This makes the joint inversion method feasible. Finally, we investigate the geometry, P-wave velocity structure and anisotropy of the crust beneath Southeastern China by applying the proposed inversion method to previously acquired wide-angle seismic data. In this case, the anisotropy signature provides clear evidence that the Jiangshan-Shaoxing fault is the natural boundary between the Yangtze and Cathaysia blocks.     

几种常用地方震相对天然地震震源深度测定的误差解析分析及数值对比

本文从误差解析公式及数值模拟计算两种途径讨论了利用常用地方震相Pg,Sg,PmP,Pn,sPL测定震源深度的误差问题,结果表明,两种途径获取的误差值相当。对于上地壳的地震而言,当直达波走时误差处于0.1 s的量级时,若要将误差控制在3 km左右,则应选用震中距为30 km以内的台站;当走时误差处于0.2 s的量级时,若要控制同等误差,则应选用震中距为20 km以内的台站;如果地震位于下地壳,震中距可适当放宽,然而当震中距更大或走时误差更大时,震源深度的误差则近乎成倍增长。PmP,Pn,sPL对上地壳的震源深度测定误差要小于下地壳,同时对误差的控制较好,不会随震中距的增大而快速增大,震中距处于90 km范围以内且走时误差小于0.1 s时的深度误差基本均能控制在3.5 km以内。此外,本文还通过“棋盘格”的方式定量地分析了速度扰动对走时的影响,并以首都圈地区台网布局为基础,分析了加入首波对震源深度测定的改善效果。这两项数值对比结果均表明,在2%的速度扰动下,只要下地壳和莫霍面的速度参数不同时出现过大或过小现象,加入首波后对震源深度的测定误差则基本能控制在3 km以内,且一致性明显地高于单独使用直达波。      

Seismic detectability of anomalous structure at the top of the Earth’s outer core with broadband array analysis of SmKS phases

I have examined precisely the differential travel times and waveforms of SmKS seismic phases propagating under the southern Indian Ocean obtained from African broadband seismic arrays. The SmKS phases analyzed in this study travel in the mantle with weak heterogeneity confirmed by a global tomographic study for the distance range of 115-135°. The SmKS differential times were obtained from a vespagram (a stack intensity on a time-slowness diagram), and comparison with the vespagram created from synthetic waveforms with PREM gives the travel-time residual for each event-array pair. Although the residuals of S3KS-S2KS times exhibit apparently a systematic dependence on epicentral distance, this is likely due to small-scale heterogeneity beneath the Oceania where is covered by the SmKS ray entering points at the CMB. Waveform modeling was applied to a record section with a small travel-time residual that suggests a small effect from the mantle heterogeneity on the data set, I found that a low-velocity zone in the outermost 50 km in the core rather than PREM can explain an additional arrival detected just after the S3KS phase. This result is still inconclusive because of the small number of data and non-uniqueness of the model and ambiguity due to mantle structure. However, accumulation of the precise measurement described in this study may help the reduction of uncertainty and trade-offs.     

长乐—南澳断裂带两侧地壳结构差异的地震-重力联合反演

长乐一南澳断裂带出露于福建沿海地区,由于海陆过渡带的特殊性,地球物理探测受到许多限制,难以获得由陆到海的清晰而准确的深部构造形态.2014年福建省地震局采集了横跨长乐一南澳断裂带的广角反射/折射剖面(HX-6),由于观测系统的缺陷和原始资料信噪比等问题,单纯使用地震数据反演长乐一南澳断裂带的深部地壳结构有很强的不确定性,无法解答断裂带两侧地壳结构存在何种差异,影响了对断裂带构造属性和区域构造演化的正确认识.基于岩石波速和密度有良好的对应关系,地震-重力联合反演可以有效降低多解性.本文采用地震走时拟合和重力异常拟合同步进行的方法,利用最新采集的高质量P波地震走时数据与高精度实测重力数据,反演得到了连城一厦门一金门外海剖面的二维地壳波速-密度结构模型.联合反演结果显示:长乐一南澳断裂带两侧地壳厚度差约3 km,壳内分层结构和上地幔顶部波速-密度无显著变化,推断长乐一南澳断裂带是华南地块正常陆壳和台湾海峡减薄陆壳的分界,本研究结果为进一步研究该区深部构造环境和长乐一南澳断裂带的地球动力学意义提供了新的地球物理学证据.     

精密控制震源理论波形的计算方法

精密控制震源是精密主动监测技术的关键设备, 它为整个系统发射精确的、 高度可重复的震动信号, 是精密主动监测的基础。 其理论波形对于观测数据处理、 提取监视地下介质的信息是十分重要和关键的。 精密控制震源发出的震动信号由于受到重力的影响, 并不是标准的正弦或余弦波形。 本文研究了重力对震源发射波形的影响, 推导出在重力影响下精密控制震源理论波形的计算公式, 给出了该计算公式的数字递推求解方法和误差控制方法, 可以准确地计算出精密控制震源在重力影响下发射信号的理论波形。     

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.     

宁化—大田—惠安地壳构造与速度结构特征

福建地处欧亚大陆东南缘,新构造活动强烈,区内北东向断裂带异常发育,是华南震区中、强地震活动的频发区.为深入认识我国东南沿海地壳上地幔速度结构特征及其深部构造背景,福建省地震局联合中国地震局地球物理勘探中心于2010年至2012年在福建陆域实施由18次人工爆破、四条北西向原生纵测线和四条北东向集成纵测线构成的三维人工地震测深实验.本文对该实验中以北西-南东走向近似垂直穿过政和—海丰断裂的宁化—大田—惠安深地震测深测线数据进行处理解释,采用地震射线走时正演构建了该剖面二维地壳速度模型.结果显示,沿剖面地壳厚度由西向东逐渐减薄,其西北侧地壳厚约31.8km,东南侧地壳厚达28.4km.剖面上地壳P波速度从5.90km·s~(-1)逐渐增加至6.20km·s~(-1),上地壳厚度横向变化不大,厚度在16~17km左右,但是下地壳厚度由西向东减薄较为明显.地壳以政和—海丰断裂为界,东西两段具有明显不同的速度结构,呈西段速度偏低、东段速度较高的特性,且西段在上下地壳分界面下方存在一个低速层.研究表明,剖面不同区段呈现出的速度结构差异与该区大地构造单元的划分基本吻合,剖面解释结果和以往远震接收函数研究结果均印证了作为闽西南坳陷带和闽东火山断陷带分界线的政和—海丰断裂是一条切割至下地壳底部的深大断裂.     

SsPmp震相地壳探测方法

SsPmp波是远震S波经地表反射转换的P波在莫霍面发生反射后被地表台站接收得到的震相.震中距在30°~50°之间的远震S波震相经地表反射转换的P波射线参数较大,在莫霍面发生全反射,使得台站接收的SsPmp波具有较强的能量,能够从地震记录中清楚地识别出来,为探测台站附近的莫霍面形态提供新的途径.本文通过合成理论地震图分析了SsPmp震相与地壳厚度、射线参数和Pn波速度之间的关系.结果表明:对于水平界面,地壳厚度只影响SsPmp与Ss波之间的相对到时差;Pn波速度只影响SsPmp的相位;射线参数既对SsPmp波的相对到时有影响,也会引起SsPmp波的相位变化.对于复杂的界面,SsPmp反映的深度与速度梯度最大的深度接近,而反映的Pn波速度与实际的Pn波速度一致.     

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.     

Three-dimensional velocity model of the crust of the Bohemian Massif and its effects on seismic tomography of the upper mantle

We have compiled a representative three-dimensional P-velocity model of the crust of the Bohemian Massif (BM) to provide a basis for removing effects of the crustal structure in teleseismic tomography of the upper mantle. The model is primarily based on recently published 2D velocity models from findings of wide-angle refraction and near-vertical reflection seismic profiles of CELEBRATION 2000, ALP 2002, and SUDETES 2003 experiments. The best fitting 3D model of the BM crust (NearNeighbour model) is complemented by velocities according to the reference Earth model at sites where data are sparse, which precludes creating artificial heterogeneities that are products of interpolation method. To test the model, we have performed tomographic inversions of the P-wave travel times measured during the BOHEMA II experiment and compared the results obtained with and without crustal corrections. The tests showed that the presented crustal model decreases magnitudes of velocity perturbations leaking from the crust to the mantle in the western part of the BM. The tomographic images also indicated a highvelocity anomaly in the lower crust or just beneath the crust in the Brunovistullian unit. Such anomaly is not described by our model of the crust since no seismic profile intersects this part of the unit. The tests also indicated that crustal corrections are of the great importance especially for interpretations of the uppermost mantle down to depths of about 100 km.     

Vertically Inhomogeneous Models of the Upper Crustal Structure in the West-Bohemian Seismoactive Region Inferred from the Celebration 2000 Refraction Data

As part of the international refraction measurements in Central Europe in the year 2000, three profiles traversed the region of earthquake swarms in West-Bohemia/Vogtland. The shots were also recorded at the permanent stations of the local seismic networks. The travel times of P-waves, observed in the West-Bohemian region, are discussed and interpreted in the present paper. In general, significantly lower P-wave velocities were found in the Saxothuringian (northern) part of the studied area than in the adjacent southern parts. The observed travel times are interpreted separately for the individual geological units, in particular for the plutons, crystallinicum, and the Mariánské Lázn (Marienbad) Complex. After smoothing the selected data using rational approximations, the Wiechert-Herglotz method was used to compute vertically inhomogeneous velocity models. The characteristic features of the derived models are relatively low P-wave velocities at the surface and prominent velocity increases within the uppermost crust down to a depth of about one kilometre.     

Operational Processing of Hydroacoustics at the Prototype International Data Center

— The Prototype International Data Center (PIDC) has designed and implemented a system to process data from the International Monitoring System's hydroacoustic network. The automatic system detects and measures various signal characteristics that are then used to classify the signal into one of three categories. The detected signals are combined with the seismic and infrasonic detections to automatically form event hypotheses. The automatic results are reviewed by human analysts to form the Reviewed Event Bulletin (REB). Continuous processing of hydroacoustic data has been in place since May 1997 and during that time a large database of hydroacoustic signals has been accumulated. For a two-year period, the REB contains 13,582 T phases that are associated to 8,437 events. This is roughly 25% of REB events after taking station downtime into account. Predicted travel times used in locations are based on the arrival time of the peak a mplitude mode calculated from a normal mode propagation model. Global sound velocity and bathymetry databases are used to obtain reliable 2-D, seasonally dependent, travel-time tables for each hydroacoustic station in the PIDC. A limited number of ground-truth observations indicate that the predicted travel times are good to within 5 seconds for paths extending to over 7,000?km – corresponding to a relative error of less than 0.1%. The ground truth indicates that the random errors in measuring arrival times for impulsive signals are between 1 and 6 seconds. This paper describes and evaluates the automatic hydroacoustic processing compared to the analyst reviewed results. In addition, special studies help characterize the overall performance of the hydroacoustic network.     

Scanning anisotropy parameters in horizontal transversely isotropic media

The horizontal transversely isotropic model, with arbitrary symmetry axis orientation, is the simplest effective representative that explains the azimuthal behaviour of seismic data. Estimating the anisotropy parameters of this model is important in reservoir characterisation, specifically in terms of fracture delineation. We propose a travel‐time‐based approach to estimate the anellipticity parameter η and the symmetry axis azimuth ? of a horizontal transversely isotropic medium, given an inhomogeneous elliptic background model (which might be obtained from velocity analysis and well velocities). This is accomplished through a Taylor's series expansion of the travel‐time solution (of the eikonal equation) as a function of parameter η and azimuth angle ?. The accuracy of the travel time expansion is enhanced by the use of Shanks transform. This results in an accurate approximation of the solution of the non‐linear eikonal equation and provides a mechanism to scan simultaneously for the best fitting effective parameters η and ?, without the need for repetitive modelling of travel times. The analysis of the travel time sensitivity to parameters η and ? reveals that travel times are more sensitive to η than to the symmetry axis azimuth ?. Thus, η is better constrained from travel times than the azimuth. Moreover, the two‐parameter scan in the homogeneous case shows that errors in the background model affect the estimation of η and ? differently. While a gradual increase in errors in the background model leads to increasing errors in η, inaccuracies in ?, on the other hand, depend on the background model errors. We also propose a layer‐stripping method valid for a stack of arbitrary oriented symmetry axis horizontal transversely isotropic layers to convert the effective parameters to the interval layer values.     

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.     

以和静地震为例分析sPn震相在震源深度测定中的应用

利用新疆台网波形数据,对2012年6月30日新源、和静交界区MS6.6地震以及2012年6月30日和静MS4.2地震波形中sPn震相进行分析。通过sPn震相与Pn震相的到时差,并结合2015年新疆模型和3400走时表模型对两个地震进行震源深度测定。结果表明,采用两个模型得出的震源深度与中国地震台网中心测定深度的误差分别为1.7km与1.5km,2.4km与2.1km,其结果真实可靠。     

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.     

上月壳中的散射引起月震尾波的数值模拟研究

在阿波罗月震记录中普遍存在着强烈持久的尾波信号,这样的波形特征无法用均匀分层月球模型解释.一个普遍被接受的解释是月震尾波由月球浅层结构对月震波的散射引起.我们采用基于交错网格的伪谱和有限差分混合方法模拟研究非均匀上月壳对月震波的散射效应,在此基础上解释月震尾波的形成机制,并估计出上月壳速度扰动的强度.我们发现,在均匀分层模型基础上,进一步考虑上月壳中的非均匀结构对月震波的散射效应,能有效地解释月震信号中强烈持久的尾波.我们认为月震尾波可能是由上月壳中的低波速、低衰减和散射这三个因素的共同作用所引起.采用不同的扰动标准差模拟上月壳的非均匀性,并比较模拟波形与真实月震图的相似程度,我们发现上月壳中速度扰动的标准差应该在3%到5%之间,很可能接近于3%.