模板匹配滤波技术在地震数据处理中的应用

模板匹配滤波技术(Matched Filter Technique,MFT)是一种近几年发展迅速的可用于检测遗漏信号的算法,在图像识别、行星探测等众多领域均有广泛应用。在地震学中,该技术通过扫描连续波形找出模板事件附近发生的、与模板事件相关性较好的遗漏地震事件,可用于完善地震目录,实现对地震活动性以及区域应力状态的更好刻画。本文综述了模板匹配滤波技术在地震数据处理中的应用,首先概述基本原理,其次介绍该方法用于探测各类地震事件的实例,包括前震和余震、远程动态触发地震、诱发地震和非火山震颤等,最后展望了用于提高算法效率及探测精度的一些改进方法。     

摩擦状态-速率依从的区域地震触发模型研究

设计一个具有摩擦状态-速率依从机制模拟大地震后区域触发地震的模型, 该模型由一定数量表示地震成核断层的函数组成。 大震对区域静应力场的影响由应力阶跃载荷表现, 利用摩擦状态-速率依从的一维自由度弹簧滑块模型结果计算地震发生时间提前。 结果表明, 施加幅度为地震应力降几十分之一至百分之一的应力阶跃, 可以引起地震发生速率的明显变化, 该变化通过统计检验, 与大地震后的区域触发地震活动相似。     

On the correlation between seismicity characteristics and S-wave attenuation in the ring structures that appear before large earthquakes

We discuss seismicity characteristics in the source zones of two great earthquakes: the December 26, 2004 Sumatra (M_w = 9.0) and the November 14, 2001 Kunlun (M_w = 7.8) events. Ring structures of low magnitude seismicity have been forming prior to these earthquakes for several decades. We studied the short period shear-wave attenuation field in the area of these ring structures. The method we used is based on the analysis of the rate of attenuation for the early Sn and Lg codas to detect attenuation inhomogeneities in the uppermost mantle. We show that the ring structures have comparatively high attenuation of shear waves compared with the crustal volumes inside the rings. The fact that there is no recent volcanism in the area of the seismicity rings shows that this effect is due to a high content of free fluids in the uppermost mantle. Proceeding by analogy with our results, we identified a zone in northern Tien Shan that is anomalous for these parameters; the zone may be related to the precursory process of a large earthquake. We discuss the geodynamic mechanisms that may be responsible for fluid concentration in the seismicity rings.     

Numerical simulation of the earthquake generation process

A numerical algorithm is proposed for the simulation of the earthquake process during a seismic cycle. The algorithm is based on a heterogeneous discrete model of the fault plane and assumes there are two kinds of seismicity: background crack-like earthquakes and asperity-like events. An active zone of the fault contains an asperity distribution with a characteristic elementary area. The background seismicity randomly develops shear stress-free surfaces which tend to surround the asperities as in a 2D percolation process. The model parameters are taken from observations on the Vrancea (Romania) intermediate depth seismic region. The results emphasize the significant role of the geometry in the mechanism of the seismic failure. The algorithm predicts the nonlinear behavior in the frequency-magnitude distribution, the decrease of theb-slope associated with the asperity-like events, the magnitude range of major earthquakes, and their recurrence times.     

Effects of seismic declustering on seismic hazard assessment: a sensitivity study using the Turkish earthquake catalogue

Earthquakes trigger other earthquakes and build up clusters in space and time that in turn create a bias in seismic catalogues. Therefore, declustering is considered as a prerequisite in seismic studies, particularly for probabilistic seismic hazard analysis, not only to eliminate the bias but also to decouple mainshocks and triggered events. However, a declustering process is not a straightforward task due to the complex nature of earthquake phenomena. There exist several declustering methods that mostly employ subjective rules to distinguish between background seismicity and offsprings. Eventually, the final declustered catalogues usually deviate significantly according to the employed method. This issue is raising some concerns, such as how to select the most suitable declustering algorithm, or to assess how this selection affects seismic hazard assessment. In consequence, the main goal of this paper is to quantify the sensitivity of seismic hazard assessments to different declustering techniques. Accordingly, the recently compiled Turkish earthquake catalogue was declustered by making use of three declustering algorithms. A total of six declustered catalogues, two catalogues per method, one by implementing the default input parameters, and one by altering the free input parameters of the employed methods, were produced. The clusters of selected earthquakes were studied in terms of the spatial–temporal distribution of earthquake sequences. A sensitivity analysis was conducted through the major steps of seismic hazard assessment for Istanbul metropolitan city. The seismicity of Istanbul and surroundings was modeled on the basis of four areal source zones. Comparative studies showed that, while the selected declustering algorithm did not significantly affect the completeness periods of moderate to large size earthquakes, it considerably altered those of small magnitude events (e.g. M_w 4.3–5.2) and consequently the recurrence parameters of the source zones. Depending on the declustering algorithm and input parameters, the activity rate was observed to vary up to a factor of two. The differences in the declustered catalogues obtained from different declustering approaches resulted in considerable variations in seismic hazard estimations. The hazard maps at return periods of 475 and 2475 years indicated that peak ground acceleration values may vary up to 20% at some locations. Moreover, the differences in 5% damped elastic spectral accelerations at T = 0.2 for the return periods of 475 and 2475 years are about 18 and 12%, respectively, on the southern shores of Istanbul where the highest hazard levels are observed.     

Triggered and tectonic driven earthquakes in the Koyna–Warna region,western India

Two strong M?>?5.0 earthquakes within a span of six months occurred in a triggered seismicity environment in the Koyna–Warna region in western India in 2000. The region is experiencing continued seismicity since the last five decades indicating that this region is close to critical stresses and minor perturbations in the stresses due to reservoir loading and unloading can trigger earthquakes. In the present study we applied the technique developed for identification of prognostic anomalies for tectonic earthquakes to the Koyna–Warna catalogue prior to these two earthquakes with an aim to study the process of source preparation for triggered earthquakes. In case of tectonic earthquakes, unstable conditions in a source zone develop gradually leading to a metastable zone which shows variations in certain seismicity parameters known as prognostic anomalies. Our results indicate that the variations in seismicity parameters before the two strong earthquakes in the Koyna region have a pattern of prognostic anomalies typical of tectonic earthquakes. We conclude that initiation of failure in a metastable zone can be caused both, by external impacts, reservoir loading and unloading in our case, and internal processes of avalanche-like failure development.     

The precursory phenomena before the Karymskii seismovolcanic crisis in parameters of seismicity in a wide range of energy

We report extensive anomalies identified in seismicity parameters at different energy levels which were observed during the precursory process of the Karymskii seismovolcanic crisis of January 1–2, 1996. The seismicity of different energies includes earthquakes contained in the Kamchatka regional catalog and seismic noise (amplitudes of 10^?9–10^?12 m, frequencies of a few tens of hertz), which is a manifestation of the seismic process in the lowest energy range. The parameters of background seismicity are considered in retrospect using techniques for analyzing the dynamics of the seismic process: RTL and the Z function. Microseismicity is examined using these authors’ own method based on monitoring the response of high frequency seismic noise to tidal excitation     

Prediction of magnitude of the largest potentially induced seismic event

We propose a method for determining the possible magnitude of a potentially largest induced seismic event derived from the Gutenberg–Richter law and an estimate of total released seismic moment. We emphasize that the presented relationship is valid for induced (not triggered) seismicity, as the total seismic moment of triggered seismicity is not bound by the injection. The ratio of the moment released by the largest event and weaker events is determined by the constants a and b of the Gutenberg–Richter law. We show that for a total released seismic moment, it is possible to estimate number of events greater than a given magnitude. We determine the formula for the moment magnitude of a probable largest seismic event with one occurrence within the recurrence interval (given by one volumetric change caused by mining or injecting). Finally, we compare theoretical and measured values of the moment magnitudes of the largest induced seismic events for selected geothermal and hydraulic fracturing projects.     

基于经验模型和物理模型研究2013 M_S7.0芦山地震余震序列

基于中国地震台网中心2013 MS7.0芦山地震余震数据我们首先确定了余震空间分布范围并根据G-R关系计算了主震后半小时内的完备震级Mc=3.5,并且得到了ML≥3.5和ML≥3.0的地震在2001年至芦山地震前的背景场地震发生率.通过Omori-Ustu经验定律和两种Dieterich模型对芦山地震余震发生率的拟合,我们发现阶梯型Dieterich模型只能模拟p=1的情况,从而造成了模拟曲线与观测数据的差别;前人研究表明震后滑移同样是产生余震的原因,如果假设余震序列由主震静态剪应力Δτ和震后滑移共同作用所产生,我们数值模拟得到的对数型Dieterich模型能够较好地推断余震发生率R随时间t增加而衰减的趋势,能够从物理机制上解释MS7.0芦山地震余震序列衰减指数大于1这一现象.通过对数型Dieterich模型的拟合并结合Andrews的方法,我们还得到MS7.0芦山地震Aσ约为0.155 MPa,ta约为8.4年,这一值与前人研究结果十分接近.     

2010年智利8.8级地震在北京房山岩体附近的动态触发活动

本文采用波形分析和β统计方法,基于首都圈地区数字测震台网和部分流动地震台的观测资料,对2010年2月27日M_w8.8智利大地震在北京房山岩体附近地区的动态触发活动进行了分析讨论.结果表明,在背景地震活动相对较弱的北京房山岩体附近,我们检测出至少5次小震活动事件在智利地震面波到达时发生.这些触发活动显著地改变了房山岩体地区这一地震活动相对平静地区的地震活动性,但是对北京地区地震活动性的影响并不明显.NKY地震台记录的智利地震触发活动的最大动态应力与之前的研究结果相比要小,约为7 kPa.这可能与房山岩体附近地区的背景地震在智利地震前一直较为平静有关.此外,相对有利的面波入射方向,以及在对蹠点上的前两组面波叠加的效应等因素综合作用,使得本研究能在约2万公里极远处观测到远震触发小震活动.然而,在随后的多次面波叠加期间并没有观测到明显的触发地震活动.在智利地震10 h前发生的琉球M_w7.0地震也未在房山地区触发相关的微震活动,这可能与它们引起的动态应力变化太过微弱有关.     

中强地震活动区地震年平均发生率确定方法的比较

Two kinds of methods for determining seismic parameters are presented,that is,the potential seismic source zoning method and grid-spatially smoothing method. The Gaussian smoothing method and the modified Gaussian smoothing method are described in detail, and a comprehensive analysis of the advantages and disadvantages of these methods is made. Then,we take central China as the study region,and use the Gaussian smoothing method and potential seismic source zoning method to build seismic models to calculate the mean annual seismic rate. Seismic hazard is calculated using the probabilistic seismic hazard analysis method to construct the ground motion acceleration zoning maps. The differences between the maps and these models are discussed and the causes are investigated. The results show that the spatial smoothing method is suitable for estimating the seismic hazard over the moderate and low seismicity regions or the hazard caused by background seismicity; while the potential seismic source zoning method is suitable for estimating the seismic hazard in well-defined seismotectonics. Combining the spatial smoothing method and the potential seismic source zoning method with an integrated account of the seismicity and known seismotectonics is a feasible approach to estimate the seismic hazard in moderate and low seismicity regions.     

RELOCATION OF THE BACKGROUND SEISMICITY AND INVESTIGATION ON THE BURIED ACTIVE FAULTS IN SOUTHEASTERN CHINA

Most of the regions in southeastern China are covered by thick Cenozoic sediments, or are the mountainous areas, so it is difficult to find and locate the active faults using the conventional geologic methods. The precisely relocated background seismicity in the seismically active region can be used to identify the buried active structure. In this paper, we investigated the relationship between regional tectonics and background seismicity, and interpreted the possible buried active faults in southeastern China using the relocated background seismicity. We relocated the background seismicity occurring in the region from 106°E to 122°E and from 22°N to 35°N between 1990 and 2014 using the doubble difference earthquake location algorithm. More than 51000 small earthquakes were relocated. In general, the relocated background seismicity corresponds well to the tectonics, showing the zonation features with typical seismicity pattern in each tectonic regime. It is observed that in the weakly active tectonic regime, the seismicity distributes dispersely or even scarcely, while in the strongly active tectonic region, the seismicity is highly clustered and organized to lineation pattern showing the same direction as the strike of the dominating fault zone. We interpreted the buried active faults using the lineation of seismicity. The inferred active faults are observed in the southeast coast region, the northwest Guangxi Province, the southeast boundary region of the Sichian Basin, and around the Huoshan Fault, many of which were not found by previous studies. The relocated hypocentral depth varies greatly in southeastern China. The shallowest earthquakes between 0 and 15km mainly distribute in the central region, indicating that the brittle deformation process only occurred in the upper crust, while the middle and lower crust are to be half-ductile and ductile deformation. There are earthquakes occurred in lower crust in the southeast coast region. The maximum depths distribute in the southeast boundary region of the Sichuan Basin, some are greater than 40km, indicating that the crust depth is larger than other places and the lower crust still sustains brittle deformation, which corresponds to the lower geothermal value and high crustal strength.     

华北地区中强震前中小地震活动的动态演化特征研究

应用地震活动加速指数AI算法对华北地区1970年以来M_L≥6.0地震进行了回溯性研究。 利用该方法统计中强震前的中小地震活动频度, 分析结果显示在空间上能够反应出中小震活动频度的异常。 另外, 应用有关地震活动的频度面积S值方法, 对M_L6.0级以上强震前后的中小地震活动在时间分布上进行了研究, 此方法对AI算法的计算结果在时间上给予了支持。 结果表明: AI指数能够量化地给出中强震前中小地震活动的相对增强或减弱程度, 中强震的S值在震前表现为高值或低值异常, 震后下降或上升。 S值的异常点与AI扫描值的异常结果在时空分布上相辅相成。 在地震目录完整的情况下, 中强震震中区及其附近地区震前1 a的中小地震活动与之前3 a、 5 a的活动背景相比, 出现明显的“加速”或“减速”现象。 利用这两个参量的变化, 可为该地区中强震的预测分析提供一定的判定指标。     

统计预测模型对三峡库区诱发地震预测的应用

通过实地考察并结合前人的研究成果,在地震地质背景研究的基础上,将长江三峡工程坝址至巫山库段划分为5个预测区共12个预测单元,选用8个诱震因素,运用统计预测模型预测了水库诱发地震的可能性及诱发地震震级。结果表明,库水蓄至最高水位,三峡工程坝址至巫山库段中多数预测单元诱发微震和不发生地震的概率值高,无震可能性较大。但九湾溪断层沿线、仙女山断层库段具备诱发3.0～4.5级或4.5～6.0级地震的可能;考虑到龙会观5.1级地震与高桥断裂的关系,预测结果显示高桥断裂沿线具有诱发强烈水库地震(≥6.0级)的可能性,其预测概率值为0.23     

华北与加利福尼亚中部沿海地区地震活动性C值的对比(英文)

为了解加州中部沿海的地震活动特征,本文利用了C值方法分析了过去36年的加州中部沿海地区地震资料,并与b值和地震频度结果进行了对比。分析了36年地震活动的基本情况。华北和加州中部沿海地区的C值资料分析表明,地震的空间分布不是随机的也不服从泊松分布。C值、b值和地震频度这三个定量的参数能够描述过去36年里加州中部沿海地区的地震活动的基本特征。     

地下水动态趋势和短临异常

在无其他因素引起水量增、减条件下,地下水动态大于一年的趋势变化是地下水动态的异常变化。识别地下水动态异常变化的条件是调查观测井的水文地质环境、水动力条件及水量有无变化。在趋势异常发展背景下地下水动态趋势的转折与地震活动有关。列举了河间等三口井近年来的水位观测结果,表明短临异常在震前均有趋势变化背景,因此短临异常的判别要以趋势变化分析为基础。地下水动态趋势异常和短临异常随时间的分布是有序的,用序列异常方法分析异常与地震活动形势,研究序列异常发生、发展、演变的规律与地震活动的关系,对地震发生的时间和强度有较好的预报意义。     

Mathematical basis, correlativity and selection of seismicity parameters

Using the analytic expression of seismicity parameters, the dependency and correlativity between the statistical parameters and seismic frequency or seismic intensity are discussed. The statistical parameters are divided into two kinds. One kind is the regional seismicity parameter, 17 parameters are analyzed in this paper. The other kind is the seismicity distribution parameter. They are the distribution parameters of time, space and magnitude. The existent base and rationality of distribution parameters depend on the application of distribution model. We analyze and draw an analogy between the natural probability, Poisson, Weibull distributions and multi-fractal analytic formula in time, space and magnitude. And some examples are given in this paper. The P value and H value of aftershock sequence attenuation, the U value and F value of earthquake swarm sequence and the entropy of information are discussed preliminarily. Another method about analyzing relationships among time-series curves are given. The resemblance relativity degree and the relativity degree of relative change slope can be used as the determining values. At last, some preliminary ideas about sifting and using for seismicity parameters are advanced in this paper.     

A Hidden Markov Model Based Tool for Geophysical Data Exploration

--Unsupervised learning techniques provide a way of investigating scientific data based on automated generation of statistical models. Because these techniques are not dependent on a priori information, they provide an unbiased method for separating data into distinct types. Thus they can be used as an objective method by which to identify data as belonging to previously known classes or to find previously unknown or rare classes and subclasses of data. Hidden Markov model based unsupervised learning methods are particularly applicable to geophysical systems because time relationships between classes, or states of the system, are included in the model. We have applied a modified version of hidden Markov models which employ a deterministic annealing technique to scientific analysis of seismicity and GPS data from the southern California region. Preliminary results indicate that the technique can isolate distinct classes of earthquakes from seismicity data.     

Spatial Separation of Large Earthquakes, Aftershocks, and Background Seismicity: Analysis of Interseismic and Coseismic Seismicity Patterns in Southern California

We associate waveform-relocated background seismicity and aftershocks with the 3-D shapes of late Quaternary fault zones in southern California. Major earthquakes that can slip more than several meters, aftershocks, and near-fault background seismicity mostly rupture different surfaces within these fault zones. Major earthquakes rupture along the mapped traces of the late Quaternary faults, called the principal slip zones (PSZs). Aftershocks occur either on or in the immediate vicinity of the PSZs, typically within zones that are ??2-km wide. In contrast, the near-fault background seismicity is mostly accommodated on a secondary heterogeneous network of small slip surfaces, and forms spatially decaying distributions extending out to distances of ??10?km from the PSZs. We call the regions where the enhanced rate of background seismicity occurs, the seismic damage zones. One possible explanation for the presence of the seismic damage zones and associated seismicity is that the damage develops as faults accommodate bends and geometrical irregularities in the PSZs. The seismic damage zones mature and reach their finite width early in the history of a fault, during the first few kilometers of cumulative offset. Alternatively, the similarity in width of seismic damage zones suggests that most fault zones are of almost equal strength, although the amount of cumulative offset varies widely. It may also depend on the strength of the fault zone, the time since the last major earthquake as well as other parameters. In addition, the seismic productivity appears to be influenced by the crustal structure and heat flow, with more extensive fault networks in regions of thin crust and high heat flow.     

紫坪铺水库蓄水前天然地震活动

紫坪铺水库位于龙门山断裂带。了解和掌握水库蓄水前库坝区及其周围天然地震活动背景水平,可对水库诱发地震的危险性进行前期评价,可为蓄水后的诱发地震活动监测提供可靠的鉴别依据。紫坪铺水库地震台网在水库蓄水前运行逾1年,本文用这批资料提供表征库区及附近和库坝区蓄水前天然地震活动水平确切实用的各项指标。用全国和四川台网地震资料,给出含库区的较大区域地震活动背景和对库区的影响。