由地方台网记录确定远震信号的方位角和慢度

介绍一种新的适合于地方遥测台网数据处理的方法，估算远震信号到达台站的方位角和慢度．该方法是根据信号在各台站上的到时与台站位置矢量在信号传播方向上的投影之间的相关性的原理．实际分析结果表明，应用该方法对地方台网的记录进行处理时可以准确地计算出信号的方位角和慢度，并能准确快捷地对地方台网记录的远震信号进行解释．      

上海地震台阵的地震定位方法

上海地震台阵数据处理软件系统的地震定位方法,采用台阵的聚束方法得到地震方位角和视慢度,根据统计得到的视慢度——震中距表推算震中距。并结合了地震台网的定位方法,由单台记录的各类主要震相从JB走时表得到震中距,然后进行地震定位。该定位方法可对近震、远震进行定位处理,并由深震相得到震源深度。     

Slowness and azimuth determination for Bucovina array (BURAR) applying multiple signal techniques

The BURAR seismic array, located in Northern Romania (Bucovina region), is designed to monitor events located in an area poorly covered by other existing seismic stations. In order to use the BURAR array for single-station locations, it is crucial to calibrate the azimuth and slowness parameters, which are currently used in array techniques to locate earthquakes, blasts or nuclear explosions. The goal of this study is to apply “f–k” and plane wave fit techniques in order to constrain the slowness and azimuth parameters at BURAR for teleseismic, regional and local events. The analysis was carried out using P and S waves recorded for events occurred between 2004 and 2008 within a radius of 50° around BURAR. The azimuth values obtained applying both methods strongly deviated from the theoretical values for regions like Central Turkey, Bulgaria, Dodecanese Islands and other parts of Greece, while the ray parameter deviations with respect to a 1-D IASP91 reference model are less significant. For the local events, the anomalies are smaller, except the particular case of Vrancea intermediate-depth earthquakes for which strong azimuth deviations (33.5°), both positive and negative, are observed. We investigate how these systematic deviations in azimuth are explained by the structure lateral heterogeneities which characterize the study region.     

格尔木地震台阵勘址数据分析与台阵布局设计

为了增强西部地区的地震监测能力, 中国拟在格尔木地区建立一个小孔径地震台阵。 本文对台阵勘址数据进行了噪声与信号的相关性分析, 得到地震监测的最佳台间距。 结果表明, 对于近震和区域震的监测, 该台阵子台间距最好小于500 m, 对于远震的监测, 台站组合间距应为1500～2000 m。 最后将勘址布设的台阵作为初选台阵进行了台阵响应计算, 计算显示, 台阵响应的主瓣在NW-ES方向较窄, 表明对来自该方向事件的慢度分辨率较高; 由于呈“L”型分布, 该初选台阵确定某些方向地震的方位角较好, 但检测其他方向事件的方位角精度不高, 这可以通过台阵校正进行改善; 台阵响应中出现的多处侧瓣是由于子台间距较大造成的。     

地震台阵上信号方位角和慢度的时、频域估计方法比较

地震台阵常用频率-波数分析法来估计地震信号的后方位角和慢度。尽管有关的算法是众所周知的,但不同的实现在某些情况下可能导致不同的结果。如应用在宽频带的澳大利亚WRA台阵的记录时,标准的f-k分析方法往往给出不正确的结果。本文发现错误是由FFT的频谱泄漏效应引起的,如在进行FFT之前对原始数据进行高通滤波则可以有效地避免这样的错误。进一步对在时域中直接估算信号方位角和慢度的一种算法进行了分析,比较发现时域方法的计算速度和精度都不低于频域方法,且在某些低信噪比的情况下,前者可以给出更可靠的结果。     

Seismic detections of the 15 February 2013 Chelyabinsk meteor from the dense ChinArray

Chin Array is a dense portable broadband seismic network to cover the entire continental China, and the Phase I is deployed along the north-south seismic belt in southwest China. In this study, we analyze seismic data recorded on the Chin Array following the February 15,2013 Chelyabinsk(Russia) meteor. This was the largest known object entering the Earth's atmosphere since the1908 Tunguska meteor. The seismic energy radiated from this event was recorded by seismic stations worldwide including the dense Chin Array that are more than 4000 km away. The weak signal from the meteor event was contaminated by a magnitude 5.8 Tonga earthquake occurred *20 min earlier. To test the feasibility of detecting the weak seismic signals from the meteor event, we compute vespagram and perform F-K analysis to the surface-wave data. We identify a seismic phase with back azimuth(BAZ) of 329.7° and slowness of 34.73 s/deg, corresponding to the surface wave from the Russian meteor event(BAZ *325.97°). The surface magnitude(MS) of the meteor event is 3.94 ± 0.18. We also perform similar analysis on the data from the broadband array F-net in Japan, and find the BAZ of the surface waves to be316.61°. With the different BAZs of Chin Array and F-net,we locate the Russian meteor event at 58.80°N, 58.72°E.The relatively large mislocation(*438 km as compared with 55.15°N, 61.41°E by others) may be a result of thebending propagation path of surface waves, which deviates from the great circle path. Our results suggest that the dense Chin Array and its subarrays could be used to detect weak signals at teleseismic distances.     

那曲、和田台阵对汶川地震序列的定位误差校正

We use the slowness-azimuth station correction （SASC） method to improve the location accuracy of the Wenchuan aftershocks recorded by the Nagqu and Hotan seismic arrays. The results show that the standard deviations of back-azimuth and slowness errors of Wenchuan aftershocks recorded by the Nagqu array decreased by 32% and 58% respectively after correction. The decrease is 38 % and 71% for the Hotan array. After the correction, the location accuracy of all Wenchuan aftershocks recorded by the Nagqu array is improved. For the Hotan array, the accuracy is improved in the slowness estimation for 78 % of aftershocks and in back-azimuth estimation for all aftershocks.     

Slowness Corrections — One Way to Improve IDC Products

—?The first step to identify and locate a seismic event is the association of observed onsets with common seismic sources. This is especially important in the context of monitoring the Comprehensive Nuclear-Test-Ban Treaty (CTBT) at the International Data Center (IDC) being developed in Vienna, Austria. Well-defined slowness measurements are very useful for associating seismic phases to presumed seismic events.¶Shortly after installation of the first seismic arrays, systematic discrepancies between measured and theoretically predicted slowness values were observed, and therefore slowness measurements of seismic stations should be calibrated. The observed slownesses measured with small aperture arrays, some of which will be included in the International Monitoring System (IMS) now being implemented for verifying compliance with the CTBT, show large scatter and deviations from theoretically expected values. However, in this study a method is presented, by which mean slowness corrections can be derived, which show relatively stable patterns specific to each array.¶The correction of measured slowness values of these arrays clearly improved the single array location capabilities. Applying slowness corrections with seismic phases observed by ARCES, FINES, GERES, and NORES, and associated to seismic events in the bulletins of the prototype International Data Center (pIDC) in Arlington, VA, also clearly demonstrates the advantages of these corrections. For arrays with large slowness deviations that are due to the influence of a dipping layer, the corrections were modeled with a sine function depending on the measured azimuth. In addition, the measured values can be weighted with the corresponding uncertainties known from the process of deriving the mean corrections.     

上海地震台阵的标定方法

地震台阵标定可以校正视慢度和后方位角，从而进一步校正震源位置，有助于提高地震定位精度及地震类型鉴别研究等。此外，台阵标定也可校正由于台阵接收系统而导致视慢度和后方位角的偏差。     

A shallow attenuating anomaly inside the ring fracture of the Valles Caldera, New Mexico

Spectral ratios of teleseismic direct and scattered P waves observed in the Valles Caldera, New Mexico, show a systematic pattern of low amplitudes at sites inside the caldera relative to sites on or outside the ring fracture. Waveforms recorded at caldera stations are considerably more complex than those recorded outside the caldera. The data used in this study were collected during a passive seismic monitoring experiment conducted in 1987. Twenty-four teleseismic events were recorded on two linear arrays spanning the caldera. To first order, the pattern of low amplitudes did not vary with source incidence angle or azimuth of approach, and could not be explained by anomalous amplification at the ring fracture. This observation suggests the presence of a shallow, attenuating zone associated with the caldera fill material inside of the ring fracture. We estimated the general features of the caldera's near-surface structure for the two-dimensional vertical cross section beneath the array, using a modification of the Aki-Larner discrete-wavenumber method to forward model the observed amplitude variations. Our results indicate that the caldera fill material must be subdivided into at least two distinct zones: a strongly attenuating lower zone, extending to depths in excess of 4 km, and a mildly attenuating surface layer. To fit the data we had to assign an unrealistically low value to seismic Q in the deeper attenuating anomaly. We attribute this to the inability of the Aki-Larner method to account for strong re-direction of energy away from the caldera due to local heterogeneity that we could not include within the low-Q anomaly. This interpretation is consistent with the pervasive, fractured hydrothermal system that is known to exist in the caldera fill material.     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part I: Critical slowness match

Kinematical characteristics of reflected waves in anisotropic elastic media play an important role in the seismic imaging workflow. Considering compressional and converted waves, we derive new, azimuthally dependent, slowness-domain approximations for the kinematical characteristics of reflected waves (radial and transverse offsets, intercept time and traveltime) for layered orthorhombic media with varying azimuth of the vertical symmetry planes. The proposed method can be considered an extension of the well-known ‘generalized moveout approximation’ in the slowness domain, from azimuthally isotropic to azimuthally anisotropic models. For each slowness azimuth, the approximations hold for a wide angle range, combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. We consider two cases for the wide-angle ray: a ‘critical slowness match’ and a ‘pre-critical slowness match’ studied in Parts I and II of this work, respectively. For the critical slowness match, the approximations are valid within the entire slowness range, up to the critical slowness. For the ‘pre-critical slowness match’, the approximations are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. The critical slowness match is particularly effective when the subsurface model includes a dominant high-velocity layer where, for nearly critical slowness values, the propagation in this layer is almost horizontal. Comparing the approximated kinematical characteristics with those computed by numerical ray tracing, we demonstrate high accuracy.     

On the rotation of teleseismic seismograms based on the receiver function technique

The receiver function (RF) technique is a well-established method to investigate the crustal and upper mantle structures based on three-component seismograms of teleseismic events. In the present study, we propose a modified automatic procedure to determine the back azimuth and polarization angles of a teleseismic event based on the RF technique. The method is tested for the recording of 3 permanent and 3 temporary broadband seismic stations located in the vicinity of Poland. Additionally, the analysis of Rayleigh wave polarization is conducted to show that the new procedure is not sensitive to incorrect seismometer orientation. The synthetic modelling of RF by a modified ray-tracing method for 2.5D models beneath each seismic station down to a depth of 60 km is performed to show the effectiveness of the proposed method in the calculation of RF for a complex structure with dipping layers.     

远程事件次声监测中的信号参数二维子空间计算方法

越来越多的观测发现,在地震、火山爆发、泥石流等重大自然灾害发生前,常产生异常的次声信号,这为地震及其他自然灾害的预报工作增加了一种可能的信息;同时,次声还是监测大气层、浅地表爆炸的有效手段.在自然灾害和爆炸事件次声监测中,慢度和方位角等参数对于源信号传播、定位以及源性质识别等工作具有重要意义.然而,目前的慢度和方位角等参数的算法——频率波数(FK)分析法,尚存在精度和分辨率不高等问题,特别是对多源次声信号的识别能力较差.为提高次声信号的监测精度,基于次声信号和噪声的子空间不相关性,构建了次声信号慢度和方位角二维子空间计算模型,并在此基础上提出了一种高分辨率次声信号二维子空间算法,仿真实验和实际数据的对比分析结果表明:本文提出的方法在精度和分辨率方面明显优于FK法,且能够更好地分离多源次声信号.     

上海佘山地震台阵远震记录的压缩感知分析

将压缩感知方法应用于上海佘山地震台阵远震定位,对于该台阵记录的M_S 5.5以上全球地震事件,根据震级大小、地震波走时、事件分布,筛选得到45个远震事件记录,采用分析压缩感知及传统方法,计算最优反方位角和慢度值,发现压缩感知方法在地震台阵的远震定位中表现良好;对于震相较为复杂的地震,在求取能量最大及超过最大能量95%以上的点,得到最终源信号,也就是震源位置,压缩感知方法具有更高的分辨率。     

接收函数的克希霍夫２D偏移方法

本文将地震反射资料成像的Kirchhoff偏移方法引入到接收函数成像研究，提出了一种接收函数的Kirchhoff偏移方法，以适应介质速度的横向变化，提高接收函数成像的精度和分辨率. 模型检验表明，与传统的水平分层的共转换点偏移叠加方法相比，Kirchhoff偏移方法能够实现转换波的有效归位，消除水平分层共转换点偏移叠加造成的假象. 一个很重要的现象是，倾斜或弯曲界面的成像明显受控于远震的入射方向，其成像精度的提高有赖于不同方向接收函数的联合成像.     

Automatic onset phase picking for portable seismic array observation

Automatic phase picking is a critical procedure for seismic data processing, especially for a huge amount of seismic data recorded by a large-scale portable seismic array. In this study is presented a new method used for automatic accurate onset phase picking based on the proporty of dense seismic array observations. In our method, the Akaike＇s information criterion （AIC） for the single channel observation and the least-squares cross-correlation for the multi-channel observation are combined together. The tests by the seismic array observation data after triggering with the short-term average/long-term average （STA/LTA） technique show that the phase picking error is less than 0.3 s for local events by using the single channel AIC algorithm. In terms of multi-channel least-squares cross-correlation technique, the clear teleseismic P onset can be detected reliably. Even for the teleseismic records with high noise level, our algorithm is also able to effectually avoid manual misdetections.     

用接收函数方法研究青华地震台下方地壳厚度结构

利用小湾电站水库诱发地震台网之青华地震台记录到的93个远震地震事件,挑选其中效果较好的55个远震地震事件进行接收函数的计算,研究同一个台站下方地壳厚度随反方位角的变化情况。结果表明:①在青华地震台下方地壳平均厚度约为40.5 km;②不同方位的远震事件反映出台站下方不同方位的地壳厚度存在差异;③青华地震台下方地壳厚度存在由南向北逐渐增厚的特点,而东西向的横向变化不明显;④青华地震台反方位角在112°附近区域地壳厚度变化异常明显。     

甘肃省测震台网地震台站地震计方位角检验与校正

利用2007年8月—2013年9月甘肃省测震台网记录的549个MS≥6远震事件的P波资料,对所属44个地震台站分量方位等问题进行检核计算。考虑到地震计维修、更换等因素,提供了甘肃省测震台网地震台站按时间序列计算的方位角偏差及方位校正变化情况,以保障测震台网数据的连续性与可靠性。研究结果表明,反演的方位偏差与文献[1]的结果基本一致,表明中国"十五"数字地震网络确实存在部分台站方位偏差较大等方面问题。因此,在进行现代地震学研究中应充分考虑台站地震计方位误差较大和研究时间段内方位角变动等因素的影响。     

上海地震台阵的F-K分析

频率－波数（F－K）谱分析是地震台阵数据处理基本方法之一。采用该分析方法，可以从地脉动背景中提取有用的地震信号，提高地震的定位精度。本文列举了F－K功率谱的不同计算方法，并将该方法应用于上海地震台阵所记录的地震事件。对地震和地脉动的F－K分析结果得出，地震的方位角和慢度值分别对记录长度的分布呈线性特征，而地脉动的方位角和慢度值则呈不规则分布。     

中国数字地震台网(CDSN)单台站P波偏振分析

偏振分析可以定量描述地震波的质点运动。 P波质点运动在水平方向上发生偏振, 使得传播方向会偏离台站与地震震源之间大圆弧的方向。 P波偏振分析可以用来约束地下结构以及反映地震各向异性, 提供与剪切波分裂等手段不同的证据。 本文介绍并使用主成分分析(Principle Component Analysis, PCA)的方法, 计算了单一地震台的P波偏振, 同时, 利用谐和分析, 对台站下方的各向异性进行了分析。 将P波偏振分析应用到中国数字地震台网(China Digital Seismograph Network, CDSN)的四个台站10年左右的数据, 筛选并分析了震相清晰、 信噪比高的远震三分量初至P波的质点运动, 拟合出四个台站接收到的地震事件后方位角与P波水平偏离角度之间的三角函数曲线, 计算出拉萨台、 昆明台下方上地幔各向异性的偏振方向分别为66°和57°, 大致呈现北东东方向。