P/S震相幅值比判据对低震级地震事件的适用性检验c

基于中国新疆及中亚地区大量天然地震和地下爆炸，检验了较大震级情况下得到的区域P/S震相幅值比识别判据对低震级地震事件的适用性. 分析结果表明，对较大震级地震事件具有较好识别效果的P/S震相幅值比判据对低震级事件同样适用，而且对低震级事件也具有与较大震级事件大致相当的识别效果. 据WMQ、BLK、MUL和MAK台的25个识别效果较好的P/S震相幅值比判据统计， 大于ML3.5和ML3.5以下地震事件的误识率之差约为2个百分点.      

Observed Characteristics of Regional Seismic Phases and Implications for P/S Discrimination in the European Arctic

—?In this paper, we use data from seismic stations operated by NORSAR, the Kola Regional Seismological Centre (KRSC) and IRIS to study the characteristics of regional phases in the European Arctic, with emphasis on the P/S ratio discriminant. While the detection and location capability of the regional station network is outstanding, source classification of small seismic events has proved very difficult. For example, the m _b ?=?3.5 seismic event near Novaya Zemlya on 16 August, 1997 has been the subject of extensive analysis in order to locate it reliably and to classify the source type. We consider the application of the P/S discriminant in the context of this event and other events observed at regional distances in the European Arctic. We show that the P/S ratios of Novaya Zemlya nuclear explosions measured in the 1–3?Hz filter band scale with magnitude, indicating a need for caution and further research when applying P/S discriminants. Using mainly data from the large NORSAR array, we note that observed P/S amplitude ratios in the European Arctic show large variability for the same source type and similar propagation paths, even when considering closely spaced observation points. This effect is most pronounced at far regional distances and relatively low frequencies (typically 1–3?Hz), but it is also significant on closer recordings (around 10 degrees) and at higher frequencies (up to about 8?Hz). Our conclusion from this study is that the P/S ratio at high frequencies (e.g., 6–8?Hz) shows promise as a discriminant between low-magnitude earthquakes and explosions in the European Arctic, but its application will require further research, including extensive regional calibration and detailed station-source corrections. Such research should also focus on combining the P/S ratio with other short-period discriminants, such as complexity and spectral ratios.     

振幅衰减特性在地震与爆破识别中的应用

基于小震级地震与爆破事件,从快速识别要求出发,分析了P波初动振幅AI和P波最大振幅AP分别与S波最大振幅AS的幅值比判据及其对事件的识别能力．考虑到传播路径对地震与爆破各振幅的影响,选用合理的振幅随距离衰减公式,分别研究了P波、 S波各振幅随震中距的衰减特征,结果显示爆破振幅衰减比地震快; 在100 km处进行衰减校正后重新计算P波与S波幅值比,得到经过衰减校正后的幅值比AI/AS正确识别率从84%提高到98%,AP/AS正确识别率从92%提高到100%,表明经衰减校正后的幅值比判据可以更好地应用于小震级地震与爆破的识别中．      

Discrimination of earthquakes and quarry blasts in the eastern Black Sea region of Turkey

In recent years, a large number of quarry blasts have been detonated in the eastern Black Sea region. When these blasts are recorded by seismic stations, they contaminate the regional earthquake catalog. It is necessary to discriminate quarry blast records from the earthquake catalogs in order to determine the real seismicity of the region. Earthquakes and quarry blasts can be separated through different methods. These methods should be applied concurrently in order to safely distinguish these events. In this study, we discriminated quarry blasts from earthquakes in the eastern Black Sea region of Turkey. We used 186 seismic events recorded by the Karadeniz Technical University and Bogaziçi University Kandilli Observatory Earthquake Research Institute stations which are Trabzon, Espiye, Pazar, Borçka, Ayd?ntepe, and Gümü?hane between years of 2002 and 2010. For the discrimination of quarry blasts from earthquakes, we used both, statistical methods (calculation of the maximum ratio of S to P waves (S/P), complexity (C)) and spectral methods (spectrogram calculation). These methods included measuring the maximum amplitude S/P, C, spectral ratio, and time-frequency analysis. We especially relied on two-dimensional time-frequency analysis methods to discriminate quarry blasts from earthquakes in Turkey. As a result of this study, 68 % of the examined seismic events were determined to be quarry blasts and 32 % to be earthquakes. The earthquakes occurring on land are related to small faults and the blasts are concentrated in large quarries. Nearly 40 % of the earthquakes occurred in the Black Sea, most of them are related to the Black Sea thrust belt, where the largest earthquake was observed in the time period studied. The areas with the largest earthquake potential in the eastern Black Sea region are in the sea.     

不同地区人工爆炸与天然地震记录特征及识别研究

本文分析了河北怀来多次爆炸、河北三河采石场多次爆炸和低震级天然地震事件的记录特征和时频差异。结果显示:河北怀来爆炸的P波能量强、衰减快、S波发育弱;河北三河采石场爆炸的P波、S波主频均低于怀来爆炸,S波与面波混淆,不同震中距的台站记录低频发育明显;而天然地震的有效频带更宽,频率成分更为复杂。将Pg/Sg谱比判据应用于小震级地震与爆炸的识别中,探索交叉频带谱比对不同地区爆炸的识别。结果表明:高频（＞5 Hz）Pg/Sg谱比判据可将研究数据中的爆炸与小震级地震完全区分;与Sg低频（0—2 Hz）有关的交叉频带谱比可对两个不同地区的爆炸进行识别,交叉频带的谱比判据较传统的单一频带谱比判据能够更好地反映出不同类型事件的特征差异。      

乌兰浩特地震台典型地震事件震相特征分析

选取乌兰浩特地震台记录到的典型天然地震事件及爆破事件,根据断裂带分布及波形记录特征,分析地震及爆破事件的震相特征。结果表明,研究区域近震波形Pn、Pg、Sn、Sg震相及Sm面波较清晰,Pg、Sg波走时差一般不大于23 s,与中国地震台网中心发布的M_S震级间的偏差一般小于0.6级;远震P、S、pP、sP震相可较清晰识别,Pm、L面波记录较明显,P、S波走时差一般不小于25 s,震级偏差一般小于0.4级,且通过震相特征及实地考察发现,爆破多为霍林郭勒区域爆破事件。     

乌兰浩特地震台典型地震事件震相特征分析

选取乌兰浩特地震台记录到的典型天然地震事件及爆破事件,根据断裂带分布及波形记录特征,分析地震及爆破事件的震相特征.结果表明,研究区域近震波形Pn、Pg、Sn、Sg震相及Sm面波较清晰,Pg、Sg波走时差一般不大于23 s,与中国地震台网中心发布的MS震级间的偏差一般小于0.6级;远震P、S、pP、sP震相可较清晰识别,Pm、L面波记录较明显,P、S波走时差一般不小于25 s,震级偏差一般小于0.4级,且通过震相特征及实地考察发现,爆破多为霍林郭勒区域爆破事件.     

Calibration of P/S amplitude ratios for seismic events in Xinjiang and its adjacent areas based on a Bayesian Kriging method

Correction maps of P/S amplitude ratios for seismic events distributed in Xinjiang, China and its adjacent areas were established using a Bayesian Kriging method for the two seismic stations WMQ and MAK. The relationship between correction maps and variations of along-path features was analyzed and the validity of applying the cor- rection maps to improve performances of P/S discriminants for seismic discrimination was investigated. Results show that obtained correction maps can generally reflect event-station path effects upon corresponding P/S dis- criminants; and the correction of these effects could further reduce scatters of distance-corrected P/S measurements within earthquake and explosion populations as well as improve their discriminating performances if path effects are a significant factor of such scatters. For example, as corresponding Kriging correction map was applied, the misidentification rate of earthquakes by Pn(2~4 Hz)/Lg(2~4 Hz) at MAK was reduced from 16.3% to 5.2%.     

利用克里金技术标定新疆及附近地区P/S震相幅值比及其在地震事件识别中的应用

针对中国的乌鲁木齐台（WMQ）和哈萨克斯坦的马坎奇台（MAK），采用贝叶斯克里金技术建立了新疆及附近地区地震事件的震相幅值比的校正曲面，并分析了校正曲面与地震波传播路径的关系，以及校正曲面对于改善震相幅值比判据识别效果的有效性.分析结果表明，建立的校正曲面基本合理，结合传播路径差异和地震波传播规律，能够较好地解释校正曲面的起伏特征. 而且，在传播路径差异较大，使得震相幅值比经震中距校正后离散度仍然较高的情况下，利用校正曲面修正震相幅值比，可以进一步降低其离散度，从而改善判据的识别效果.经校正曲面修正后，对天然地震的误识率由16.3%下降为5.2％.      

朝鲜2009年和2013年两次核爆的地震学特征对比研究

利用布设在长白山地区临时地震台站接收到的朝鲜核爆的波形资料,对2009年5月25日和2013年2月12日两次朝鲜核试验的地震学特征进行比较.震中距范围从145km到420km.采用P/S型谱比值方法识别朝鲜核爆,通过与2009年3月20日长春地震和2013年1月23日沈阳地震事件的比较,表明在频率大于3 Hz时P/S型谱比值能够有效识别发生在中朝边境地区的地下核试验.选定参考台站,利用区域震相Pg波的振幅谱比值计算朝鲜核爆至各台站路径上的相对衰减.结合介质速度模型,在一定程度上反映了长白山地区衰减情况,为进一步研究长白山地区衰减层析成像提供初始模型.     

Seismic Discrimination of the May 11, 1998 Indian Nuclear Test with Short-period Regional Data from Station NIL (Nilore, Pakistan)

v—vRegional seismic discriminants for the May 11, 1998 Indian underground nuclear test(s) and earthquakes recorded at station NIL (Nilore, Pakistan) provide new data to test strategies that can be used to monitor the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Three categories of regional discriminants (ratios of P- and/or S-wave energy) were measured on short-period (0.5–6uHz) seismograms: P/S amplitude ratios (phase ratios) measured in the same frequency band, P- and S-wave spectral ratios (i.e., low frequency to high frequency for the same phase) and P/S cross-spectral ratios (i.e., low frequency S-wave to high frequency P-wave). The P/S amplitude ratios show good separation of the Indian nuclear test and regional earthquakes for Pn/Lg and Pn/Sn, however Pg/Lg does not discriminate as well. Pn/Lg and Pn/Sn discriminate well at frequencies as low as 0.5–2uHz, especially after accounting for path effects. This observation differs from previous studies that report poor separation of earthquakes and explosions at lower frequencies. The P/S amplitude ratios do not show any magnitude dependence, suggesting that forming the ratios in a fixed frequency band cancels the effects of source size-corner frequency scaling. Spatial variability of the observed discriminants arises from variations in crustal waveguide and/or attenuation structure (path propagation effects). Grouping amplitude ratios for earthquakes with paths similar to the Indian test greatly improves discrimination. Removing distance trends does not generally improve discrimination. Accounting for path effects with Bayesian kriging significantly improves discrimination. Spectral ratios (e.g., Pn [0.5–1uHz]/Pn [4–6uHz]) and cross-spectral ratios (e.g., Lg [1–2uHz]/Pn [4–6uHz]) show distance and magnitude dependence. We developed a technique for simultaneously removing the effects of distance and source size-corner frequency scaling on amplitude and spectral ratios. The technique uses a gridsearch to find several parameters that characterize the observed distance and magnitude dependence. Discrimination of the Indian test improved dramatically after the distance and magnitude trends were removed from the spectral and cross-spectral ratio data.     

朝鲜地下核试验的地震学观测

自2006年至2017年,朝鲜民主主义人民共和国在中朝边界地区的试验场进行了6次地下核试验.本文综合报道根据东北亚地区的宽频带数字地震资料利用地震学方法对这六次地下核爆炸的研究.结果表明,朝鲜地下核试验在区域台网产生的地震记录具有典型浅源爆炸的特征.针对上述资料发展了处理核爆数据的方法并据此得出各次朝鲜核爆的地震学参数,包括事件识别、当量测定、以及震中相对定位等.对6次核爆和4次天然地震P/S类型谱振幅比的统计分析表明,2 Hz以上台网平均谱振幅比可以正确地将朝鲜核爆从天然地震中识别出来,从而有效监测在朝鲜半岛进行的当量大于0.5 kt的地下核试验.同时也发现,建立在体波-面波震级比之上的识别方法不适用于朝鲜核试验场.通过建立中朝边界地区基于Lg波的体波震级系统,计算了各次朝鲜核试验的体波震级m_b（Lg）,并由此估计了它们的地震学当量,其值介于0.5 kt至60 kt之间.由于缺少爆炸埋藏深度的数据,上述当量有可能被低估,因而有必要对深度影响做进一步研究.以第一次爆炸的位置为参考震中,利用Pn波相对走时数据和高精度相对定位方法获得了各次核爆在试验场中的精确定位.     

Seismic Event Identification of Earthquakes and Explosions in Germany Using Spectral Lg Ratios

v—vAt the German NDC initial work on seismic event identification has focused on the application of spectral amplitude ratios for Lg in order to discriminate naturally occurring seismic events from other events associated with mining and quarry activities. Only about 10% of all seismic events occurring in Germany and adjacent areas are due to natural seismicity and are mostly constrained to the Alpine regions and areas along the Rhinegraben, Rhenish massif, Swabian Jura, and the Bohemian massif (Vogtland region). Using data from the broadband GRSN network, spectral amplitude ratios are calculated from maximum trace amplitudes in the 1–2uHz and 6–8uHz frequency bands, which are within the passbands of the deployed STS-2 instruments and the recorded 20uHz data streams. These amplitude ratios then must be corrected with an appropriate attenuation model in order to remove propagation paths effects. For event identification, a scoring scheme is applied across the GRSN network, based on station-dependent scoring thresholds. In a case study aimed at testing the identification scheme, events are investigated from a quarry in southern Germany that provided ground-truth information for six events in 1997 to demonstrate the suitability of this identification approach. Except for one event with a rather strong earthquake signature, i.e., a low spectral Lg ratio, these events could be screened out from the earthquake population by their large Lg ratios. In a second step, aimed at applying the identification scheme, all events in Germany and neighboring areas that occurred in 1995 were processed, with approximately 800 out of more than 1200 events showing explosion-type Lg ratios, while only 10% remain in the earthquake population. However, specific mining areas appear to consistently produce earthquake-type spectral ratios indicative of particular blasting practices.     

Comparison of the attenuation properties for two different areas in eastern Anatolia,Turkey

In this study, the attenuation properties of the crust and the quality factor of S wave in eastern Anatolia (Turkey) were determined by local earthquakes for two different areas, Oltu and Erzurum. Seismic wave attenuation can be changed with high pressure or structural effects. Therefore, we argued that the estimation of attenuation coefficient in seismic active zones in Eastern Anatolia is a very useful tool to determine seismic activities. It uses regional waveform data set from two stations, OLT and ERZ, for 95 events that occurred in these regions between 2001 and 2005. The attenuation has been determined using the Chobra–Alexeev model based on the epicenter distance–amplitude relations. This model allows for investigation of the effects of variations in attenuation properties for different areas. We introduced a new magnitude formula for these areas using the amplitude normalization methods for reference values M_L=4, so as to correct effects of the magnitudes. We also determined velocity of seismic waves. The average attenuation coefficient (α), average quality factor (Q_s) and P and S waves velocities were obtained with normalized amplitude values for Erzurum (ERZ) and Oltu (OLT) as 0.0135 km⁻¹, 37, 6.20 km/s and 3.38 km/s and 0.0151, 34, 6.13 and 3.48.     

A Seismic Model of Casing Failure in Oil Fields

—We develop a seismic model that characterises the sudden tensional failure of oil-well casings. The energy released by the rupture of a well casing is transformed into heat and seismic energy. The upper bound of the seismic efficiency of this process is estimated at about 3%. The static situation at the completion of a casing failure episode is modelled by calculating the static displacement field generated by two opposing forces separated by an arm. The azimuthal patterns of these displacements and the change in the strain and stress fields caused by the force couple are described. The dynamics of the failure episode are modelled as a dipole with a seismic moment equivalent to the product of the average drop in shear stress, the failure surface, and an arm. The radiated P and S waves have mean-square radiation pattern coefficients of 1/5 for P waves and 2/15 for S waves. The displacement field as a function of time during rupture and the spectral properties in the far field are derived. The most promising seismic parameters that can be used for distinguishing between casing failure events and other possible events are polarisation properties of S waves and S/P amplitude ratios. S-wave polarisation distinguishes between shear events and casing failure events. S/P amplitude ratios distin guish between tensile events and casing failure events.     

Amplitudes of long-period PcP and the core-mantle boundary

The amplitudes of the core reflection PcP are sensitive to the wave velocities and densities in the neighborhood of the core-mantle boundary (CMB). We study the amplitude ratio of the long-period phases PcP and P from two South American deep-focus earthquakes with favorable fault-plane solution, depth and magnitude, as recorded by WWNSS and CSN stations in North America.Comparison is made with long-period PcP/P amplitude ratios, derived from theoretical seismograms for a variety of CMB models. Models from previous studies, which were mainly derived from short-period PcP observations and which are characterized by discrete layers above the CMB, are almost all inconsistent with the long-period data. The data also discriminate against low nonzero S velocities below the CMB. Simple first-order-discontinuity models of the CMB, for instance according to the Jeffreys-Bullen earth model or according to recent models based mainly on free oscillations, explain the data reasonably well.Model improvements are attempted by varying the P-velocity gradient above the CMB. The best amplitude fit is obtained for a rather strong decrease in P velocity with depth in this zone which, however, gives no acceptable traveltime fit for PcP. The scatter in body-wave amplitudes is considerable even for long-period waves and may prevent the correct assessment of that part of the amplitude variation of a phase with distance that is due to the variation of velocities and densities with depth alone.     

基于贝叶斯原理的爆炸识别判据综合技术研究

基于贝叶斯原理，建立了地震事件识别判据的综合方法（简称贝叶斯方法）， 并针对中亚地区大量地震事件，采用该方法综合了不同台站和不同频带的P/S震相幅值比判据. 该方法的主要优点是能够给出事件属于爆炸的概率，而且充分考虑了判据之间的相关性. 理论计算和模拟数据分析结果表明， 贝叶斯方法是最优的综合方法， 其识别效果比任何单个判据好，也好于普通线性组合方法的识别效果. 针对中亚地区大量地震事件综合不同台站和不同频带的P/S震相幅值比判据的结果也表明，贝叶斯方法的误识率低于单个判据和线性组合方法的误识率.      

Multivariate discrimination technique based on the Bayesian theory

A multivariate discrimination technique was established based on the Bayesian theory. Using this technique, P/S ratios of different types (e.g., Pn/Sn, Pn/Lg, Pg/Sn or Pg/Lg) measured within different frequency bands and from different stations were combined together to discriminate seismic events in Central Asia. Major advantages of the Bayesian approach are that the probability to be an explosion for any unknown event can be directly calculated given the measurements of a group of discriminants, and at the same time correlations among these discriminants can be fully taken into account. It was proved theoretically that the Bayesian technique would be optimal and its discriminating performance would be better than that of any individual discriminant as well as better than that yielded by the linear combination approach ignoring correlations among discriminants. This conclusion was also validated in this paper by applying the Bayesian approach to the above-mentioned observed data.     

Features of different types of seismic events in China’s Capital Region

Seismic records produced by different seismic sources vary. In this study, we compared the waveform records and time-frequency characteristics of tectonic earthquakes, artificial explosions, and mine collapses in China’s Capital Region. The results show that tectonic earthquakes are characterized by stronger S-wave energy than P-wave energy, obvious high-frequency components, and wide frequency bands of P and S waves. Artificial explosions are characterized by greater P-wave amplitude than S-wave amplitude and near-station surface wave development. Mine collapses are characterized by lower overall frequency, more obvious surface waves, and longer duration. We extracted quantitative discriminants based on the analysis of different event records, with 31 feature values in 7 categories (P/S maximum amplitude ratio, high/low frequency energy ratio, P/S spectral ratio, corner frequency, duration, the second-order moment of spectrum, and energy strongest point). A comparison of the ability of these feature values to recognize distinct events showed that the 6–17 Hz P/S spectral ratio was able to completely distinguish artificial explosions from the other two types of events. The S-wave corner frequency performed relatively well in identifying all three types of events, with an accuracy of over 90%. Additionally, a support vector machine was used to comprehensively distinguish multiple features, with an accuracy for all three types of events reaching up to 100%.     

Detection of earthquake clusters on the basis of waveform similarity: An application in the monferrato region (Piedmont, Italy)

In August 2000 and July 2001 two seismic sequences, characterized by mainshocks with Ml (local magnitude) respectively 5.1 and 4.8, occurred in the Monferrato region (Italy). The regional seismic network of North-western Italy (RSNI) recorded more than 250 foreshocks and aftershocks. The routine locations, obtained from the Hypoellipse code, show a seismic activity concentrated in a circular area, of a radius of about 15 km, located near Acqui Terme, and randomly distributed in depth. Location errors, due to an uneven azimuthal station distribution of the regional seismic network, prevent recognition of the geometry of the active zone. Waveform analysis revealed the presence of several multiplets. In order to discriminate and successively relocate them, an automatic cross-correlation procedure was applied. Normalized cross-correlation matrixes, for the RSNI stations which recorded almost 90% of considered events, on the basis of a signal to noise ratio analysis, were computed using only S wave time windows. The use of a relocation procedure, based on the double-difference method which incorporates ordinary absolute travel-time measurements and/or cross-correlation differential travel-times, allowed us to overcome the constraints of the uneven distribution of stations giving a quite different frame of seismicity. The improved locations showed that the seismic activity is mainly arranged along a NE-SW oriented volume, at a depth range of 8–20 km, involving the basement crystalline units. This orientation is confirmed by the analysis of the focal mechanisms: most focal solutions show a strike slip component with one of the nodal planes consistent with the main orientation of the seismic events.