安徽数字测震台网新参评台站背景噪声分析

计算并分析安徽数字测震台网9个新参评台站数字化记录背景噪声,得到各台址背景噪声均方根RMS值、有效测量动态范围、噪声功率概率密度谱,按照地噪声水平规定,对各台进行台基噪声分类,数据表明,9个新参评台站有4个Ⅰ类台址、5个Ⅱ类台址。对于9个新参评台站噪声功率概率密度谱及各频点噪声干扰源,分析认为,大多数台站在频率10 Hz附近存在2组概率较高的背景噪声,与白天和夜晚的不同噪声水平相对应。     

首都圈数字地震台网对微弱爆破信号的检测能力

利用首都圈数字地震台网接收人工地震信号,进行地下结构研究具有重要意义.但人工震源释放的能量小,激发的地震波以短周期为主,因此本文较全面地研究了地震台网对短周期微弱信号（1～20 Hz）的检测能力：(1) 分析了台网的背景噪声,结果表明基岩台址的地震台噪声比沉积盖层台址的地震台噪声低约13 dB,这相当于近1个震级的检测阈值;夜间的噪声比白天低约5 dB;噪声有逐年增高的趋势,2006年比2001年噪声提高约4 dB.(2 )分析了在台网内进行的药量为25 kg的陆地井下爆破实验,一次爆破相当于0.69级（M_L）的天然地震,有18个地震台可辨认爆炸产生的Pg、Pm或Pc波;离爆破点218 km的基岩台,仍可以接收到振幅只有1.6 nm 的Pm波,这个结果可为地震勘探实际工作提供参考.(3) 研究了台网外核爆试验的信号特征,2006年发生在朝鲜的地下核试验是一次检验台网检测微弱信号能力的好机会.波形记录经1～5Hz滤波后,台网中噪声小的18个基岩台可以清晰辨认核爆破产生的P波或Lg波,P波平均振幅为16 nm,计算的平均震级为m_b4.3,和NEIC给出的震级相同;分析还表明背景噪声是影响台站信号检测能力的主要因素之一.     

Experimental study of the ground motion on a large scale topographic hill at Kitherion (Greece)

We present the results of an experiment about the effects on the seismic ground motion of a hill of large dimensions (6km long, 3km large and 700m high) near Corinth (Greece). We installed 7 seismometers across the hill, and analyzed the ground motion with different methods: the classical spectral ratios (CSR) and the horizontal to vertical spectral ratios calculated both on noise (HVNR) and earthquake data (RF). The amplification at the resonant frequency (around 0.7Hz) is low, as expected, but one of the two stations installed at the top exhibits large amplification (up to a factor of 10) at 3Hz. The comparison between the H/V results and the CSR ones suggests that the former method is able to point out the fundamental frequencies of a hill. Amplifications and spatial localization predicted by numerical 3D modeling are consistent with the experimental data only at frequencies lower than 1Hz.     

The East Thuringian Seismic Network

The East Thuringian Seismic Network (OTSN) was installed in 1997. It started its operation with five and now consists of six seismic stations, the GRSN (German Regional Seismic Network) station MOX and a control and analysis centre. All stations are equipped with 3-component GÜRALP and short-period seismometers, RefTek 24-bit data acquisition systems (dynamic range 23.5 bit), hard disks, GPS-receivers, modems and communication computers for dial-up purposes. The seismic signals are sampled at 100 Hz and stored on the hard disk. Simultaneously, the signals are processed by a STA/LTA detector which generates an extended event list. The central station calls these event lists once per day, analyses them, produces a list of real seismic events and calls the waveform data for these events only from the single stations. All stations operate completely autonomously and the whole system works automatically, but all operations can also be carried out interactively. The event analysis is performed manually using common seismic analysis programs. The main purpose of installing the seismic network is to investigate the local seismicity, its relation to recent tectonics, the stress field and structure of the upper crust in order to render more precisely the seismic hazard of East Thuringia. A further aim of the network is to improve the seismic monitoring situation for the neighbouring regions, especially the Vogtland/Northern Bohemia and the Western Saxony area.     

北京市测震台网数字地震台站台基背景噪声分析

基于北京市测震台网连续三分量地震计波形数据,计算28个测震台站台基噪声,利用Welch方法计算噪声功率谱密度（PSD）,进而计算地震台台基1-20 Hz地动噪声均方根值（RMS）和观测动态范围。结果表明,依照《地震台站观测环境技术要求》,北京市测震台网中有11个Ⅰ类台、9个Ⅱ类台、6个Ⅲ类台、2个Ⅳ类台。通过分析北京市测震台网数字地震台背景噪声水平,为测震台网的规划建设提供数据支持。     

P- and S-Wave Site Response of the Seismic Network RESNOM Determined from Earthquakes of Northern Baja California, Mexico

—We determined the response to P- and S-wave incidence of the permanent stations of the seismic network of Baja California (RESNOM) using two independent methods. We selected 65 events with magnitudes between 2.2 and 4.8 and hypocentral distances ranging between 5 and 330 km. The site response of the ten stations analyzed was first estimated using average spectral ratios between the horizontal and the vertical components of motion (H/V ratios). As a second approach we performed a simultaneous inversion for source and site. In order to invert the spectral records to determine the site response, we made an independent estimate of the attenuation for two different source-station path regions. Then we corrected the spectral records for the attenuation effect before we made the inversion. Although the average H/V ratio of many sites is inside the error bars of the site response estimated with the spectral inversion, the spectral inversion tends to give higher values. For the S wave some sites show similar frequency of predominant peak when comparing the responses obtained with both methods. In contrast, for the P waves the H/V ratios disagree with the results of the inversion. In general, the site response of the stations is strongly frequency dependent for both P and S waves. We also found that the natural frequency of resonance of the sites is near 0.5 Hz for P and near 0.8 Hz for the S waves.     

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.     

Broadband frequency-dependent amplification of seismic waves across Bucharest, Romania

The determination of seismic amplitude amplification is a fundamental contribution to seismic hazard assessment. While often only high-frequency amplitude variations (>1 Hz) are taken into account, we analyse broadband waveforms from 0.14 to 8.6 Hz using a temporary network of 32 stations in and around the earthquake-prone city of Bucharest. Spectral amplitudes are calculated with an adaptive multiple-taper approach. Across our network (aperture 25 km × 25 km), we find a systematic northwest/southeast-oriented structural influence on teleseismic P-wave amplitudes from 0.14 to 0.86 Hz that can be explained by constructive interference in the dipping Cenozoic sedimentary layers. For higher frequencies (1.4–8.75 Hz), more local site effects prevail and can be correlated partly among neighbouring stations. The transition between systematic and localised amplitude variations occurs at about 1 Hz.     

苏鲁地区剪切波速度结构研究及与地震关系探讨

利用S波纯波形拟合法以及T函数法反演了苏鲁地区壳幔剪切波速度结构,并利用长周期P波T函数反演得到了连云港和莱阳台下方800km深度的速度结构。结果显示:(1)苏鲁地区大部分台站地壳表层及上地壳浅部速度偏高,分别对应高压、超高压物质和古老基底出露地区;(2)沿郯庐断裂带分布的台站均显示明显低速层,并具有北浅南深的特点;(3)连云港和莱阳台超深度反演结果显示两台均在150km深度下出现高速层,反映扬子板块的俯冲深度为100km以下,俯冲板片厚度在100km以上;板片拆离下沉深度甚至达到300km或者更深;(4)地震深度分布与低速层关系密切,沿郯庐断裂和烟台—五莲断裂的中小地震震源深度都比较深,有的甚至达到地壳的底部,反映这两条断裂目前切割深度都比较大,而且地幔物质相对比较活跃。     

利用双台站背景噪声分析2017年墨西哥MW7.1地震震源区的地震波速变化

背景噪声中包含有丰富的地球内部结构信息,可以用来监测地震波速度变化.基于墨西哥中部G.UNM和MX.TLIG两个地震台站2016年4月1日至2018年4月1日的连续波形资料,本文首先对其进行不同频带下的背景噪声互相关处理,然后采用小波变换对获取到的互相关函数进行滤波来解决台站间距过长所引起的互相关函数不稳定问题,最后分别通过Stretching方法和MWCS方法测量得到的走时偏移来提取地震波速度变化.结果表明：利用小波变换能有效提高互相关信噪比和结果精度;采用Stretching方法发现震源区浅部5~10 km的地震波速在地震中上升约0.07%,而在2~5 km的地震波速则有约0.05%的下降趋势;MWCS方法计算结果稳定性稍差,但变化趋势与Stretching结果相同.造成该地震波速变化的主要原因可能是2017年地震所引起的脆性岩层结构破裂.     

松辽盆地地层吸收特性和地震波衰减规律研究

选取有代表性的岩石样本，实验室模拟地下岩层的实际物理状态，对P波的品质因素Q和波速进行了测试，结果表明，松辽盆地地下岩层的P波Q和波速皆随压力的增中呈线性增加，随温度的增加呈线性减小；且Q与波速符合幂律关系，粘弹介质地震波场正演模拟表明，导致地震波衰减的主要是第四系低速带的影响，通过对震测井资料估算，按目前的地震采集方式，即潜水面（低速带）以下激发，检波器地面接收，从激发到接收，地震波在第四系中高频成分将损失170Hz，在最理想情况下，爆炸震源激发的初始子波频宽为402Hz，经过第四系和下伏地层吸收，第四系以下的反射波频带宽度不会超过232Hz，这是松辽盆地受地层吸收限制的高频采集极限。     

山东数字测震台站台基背景噪声分析

介绍了山东数字地震台网基本情况,计算了40个测震台站台基背景噪声,利用Welch方法计算噪声功率谱密度（PSD）,进而计算地震台台基1—20 Hz地动噪声均方根值（RMS）和有效动态观测范围。根据计算结果,依照《地震台站观测环境技术要求》,对山东测震数字台网40个参评测震台站进行背景噪声级别分类,并分析不同台站背景噪声水平较低的原因,以期为测震台网的优化建设提供数据支持。     

Signal Processing for Indian and Pakistan Nuclear Tests Recorded at IMS Stations Located in Israel

v—vIn compliance with the Comprehensive Nuclear-Test-Ban-Treaty (CTBT) the International Monitoring System (IMS) was designed for detection and location of the clandestine Nuclear Tests (NT). Two auxiliary IMS seismic stations MRNI and EIL, deployed recently, were subjected to detectability, travel-time calibration and discrimination analysis. The study is based on the three recent 1998 underground nuclear explosions: one of India and two of Pakistan, which provided a ground-truth test of the existing IMS. These events, attaining magnitudes of 5.2, 4.8 and 4.6 correspondingly, were registered by many IMS and other seismic stations.¶The MRNI and EIL broadband (BB) stations are located in Israel at teleseismic distances (from the explosions) of 3600, 2800 and 2700ukm, respectively, where the signals from the tests are already weak. The Indian and the second Pakistan NT were not detected by the short-period Israel Seismic Network (ISN), using standard STA/LTA triggering. Therefore, for the chosen IMS stations we compare the STA/LTA response to the results of the more sensitive Murdock-Hutt (MH) and the Adaptive Statistically Optimal Detector (OD) that showed triggering for these three events. The second Pakistan NT signal arrived at the ISN and the IMS stations in the coda of a strong Afghanistan earthquake and was further disturbed by a preceding signal from a local earthquake. However, the NT signal was successfully extracted at EIL and MRNI stations using MH and OD procedures. For comparison we provide the signal analysis of the cooperating BB station JER, with considerably worse noise conditions than EIL and MRNI, and show that OD can detect events when the other algorithms fail. Using the most quiet EIL station, the most sensitive OD and different bandpass filters we tried in addition to detect the small Kazakh chemical 100-ton calibration explosion of 1998, with magnitude 3.7 at a distance approaching 4000ukm. The detector response curve showed uprising in the expected signal time interval, but yet was low for a reliable decision.¶After an NT is detected it should be recognized. Spectra were calculated in a 15-sec window including P and P-coda waves. The spectra for the first Pakistan NT showed a pronounced spectral null at 1.7uHz for all three components of the EIL station. The effect was confirmed by observation of the same spectral null at the vertical component of the ISN stations. For this ground-truth explosion with a reported shallow source depth, the phenomenon can be explained in terms of the interference of P and pP phases. However, the spectral null feature, considered separately, cannot serve as a reliable identification characteristic of nuclear explosions, because not all the tests provide the nulls, whereas some earthquakes show this feature. Therefore, the multi-channel spectral discrimination analysis, based on a spectral ratio of low-to-high frequency energy (in the 0.6–1uHz and 1–3uHz bands), and a semblance of spectral curves (in the 0.6–2uHz band), was conducted. Both statistics were calculated for the vertical component of the ISN stations as well for the three components of the EIL station. The statistics provided a reliable discrimination between the recent NT and several nearby earthquakes, and showed compliance with the former analysis of Soviet and Chinese NT, where nuclear tests demonstrated lower values of energy ratio and spectral semblance than earthquakes. ¶Accurate location of NT requires calibration of travel time for IMS stations. Using known source locations, IASPEI91 travel-time tables and NEIC origin times we calculated expected arrival time for the P waves to the EIL and MRNI stations and showed that the measured arrival time has a delay of about 4 sec. Similar results were obtained for the nearby Pakistan earthquakes. The analysis was complimented by the P travel-time measurements for the set of Semipalatinsk NT, which showed delays of about 3.7usec to the short-period MBH station which is a surrogate station for EIL. Similar delays at different stations evidence a path- rather than site-effect. The results can be used for calibration of the IMS stations EIL and MRNI regarding Asian seismic events.     

Regionalized difference of the 660 km discontinuity beneath Izu-Bonin

Introduction The detail structure of the Earths interior is a major field of geophysics and the existences and properties of mantle discontinuities are its important content. Based on the seismic tomographic results (Fukao, et al, 1992, Van der Hilst, et al, 1991) and the distribution of earthquake hypocenters in the Circum-Pacific region and other regions, the effects of subducting slab at different regions on the 660 km discontinuity are different: some subducting slabs penetrate through; …     

Evidence for fault-related directionality and localized site effects from strong motion recordings of the 2003 Boumerdes (Algeria) earthquake: Consequences on damage distribution and the Algerian seismic code

The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (M_w) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (M_w) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.     

Detailed gravity survey to help seismic microzonation: Mapping the thickness of unconsolidated deposits in Ottawa, Canada

In this study, measurements of gravity were made to map and model the thickness of Quaternary deposits (sand and clay) overlying Ordovician limestones in a suburb of Ottawa (Orléans, Ontario). Because ground motion amplification is partly related to the thickness of unconsolidated deposits, this work helps refine the assessment of the earthquake damage potential of the area. It also helps the mapping of clay basins, which can locally exceed 100 m in thickness, where ground motion amplification can occur. Previous work, including well log data and seismic methods, have yielded a wealth of information on near-surface geology in Orléans, thereby providing the necessary constraints to test the applicability of gravity modeling in other locations where other methods cannot always be used. Some 104 gravity stations were occupied in an 8 × 12 km test area in the Orléans. Stations were accurately located with differential GPS that provided centimetric accuracy in elevation. Densities of the unconsolidated Quaternary deposits (Champlain Sea clay) determined on core samples and densities determined on limestone samples from outcrops were used to constrain models of the clay layer overlying the higher density bedrock formations (limestone). The gravity anomaly map delineates areas where clay basins attain > 100 m depth. Assuming a realistic density for the Champlain Sea clays (1.9-2.1 g/cm³), the thickness over the higher density bedrock formations (Ordovician carbonate rocks) was modeled and compared with well logs and two seismic reflection profiles. The models match quite well with the information determined from well logs and seismic methods. It was found that gravity and the thickness of unconsolidated deposits are correlated but the uncertainties in both data sets preclude the definition of a direct correlation between the two. We propose that gravity measurements at a local scale be used as an inexpensive means of mapping the thickness of unconsolidated deposits in low-density urban areas. To obtain meaningful results, three conditions must exist. Firstly, elevations of gravity stations must be measured accurately using differential GPS; secondly, that the regional gravity field must be well defined, and thirdly, that the local geology be simple enough to be realistically represented with a two-layer model.     

利用数字化地震波形资料研究山西中南部地区的非弹性衰减系数、震源参数和场地响应

利用遗传算法反演山西中南部地区14个地震震源谱的低频水平、拐角频率和7个台站的场地响应。结果表明,Q(f)与频率的关系是:Q(f)=299.4f0.563;7个基岩台中岢岚、东山、阳城台在3 Hz~6 Hz附近有6~8倍的放大、离石台在2 Hz~10 Hz有2~3倍的放大。标量地震矩M。与ML震级线性相关较好,震源半径在230 m~508 m之间;拐角频率和地震矩之间的关系表现为拐角频率随地震矩的减小而增大;应力降和地震矩之间存在着应力降随地震矩的增大而增大的关系。     

内蒙古区域背景噪声特征分析

应用PDF方法,计算了内蒙古现运行48个测震台站0.01~20 Hz频带范围内的功率谱密度(PSD)和1~20 Hz频带范围内噪声均方根(RMS)值,定量分析了内蒙古区域背景噪声水平。结果显示:平均噪声水平属于Ⅰ类的台站有45个,Ⅱ类有3个;台站背景噪声在1 Hz以上频段内,主要受公路和人为影响;在0.6~1 Hz频段内背景噪声水平差异较小;在低频段,水平向受温度和湿度影响大于垂直向,山洞台受影响小于地面台。     

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.     

Frequency dependence of site amplification factor in the Messina Strait area

The spectral ratio method is used here for evaluating the frequency dependence of a site amplification factor in the Messina Strait area. All stations exhibit the same features of spectral ratios evidencing two peaks at about 10 Hz and 16 Hz. We relate this observation to the same geological structures. In fact, all the stations (except the referenced one) are situated over pleistocenic sediments with a similar grain size. This causes the disappearing of any site effect when the average spectrum is used as a reference. MES station, situated in the city of Messina, presents a clear site effect at a different frequency (6 Hz) which cannot be related to any geological structure because no substantial difference is observed between the sediments on which this station is situated and the sediments on which the other stations are situated. MT1 station exhibits a spectral peak at about 2 Hz which can be explained with oscillations of a little sedimentary basin.