Analysis of seismic noise in the Romanian-Bulgarian cross-border region

We present the first detailed investigation of the background seismic noise recorded in the Romanian-Bulgarian cross-border region over 3 years (2012–2015). We used the power spectral densities probability density functions (PSD PDFs) to study the noise variations in the period domain (0.025–1 s) as well as in the secondary microseism band (2–10 s). Strong diurnal variations and an increase of the noise levels during working days were observed at high frequencies at all stations, thus confirming the anthropic origin of the noise at low periods. The noise variations observed at longer periods (>?1 s) are relatively small among the stations and are related to season changes. The dominant feature in the noise spectra between 2 and 10 s is the double-frequency peak (DFP) whose amplitude increases and changes during winter. For a specific interval, from 25th to 27th of January 2014, when a storm was reported in the Black Sea area, the maximum of the DFP shifted from larger periods (~?5.5 s) at stations far from the Black Sea towards smaller periods (~?1.8 s) at stations located on the coastline. The polarization analysis showed that the short period double-frequency microseisms originating from the Black Sea dominate during the winter month. Finally, we showed that site conditions vary due to noise variations related to weather conditions in the Black Sea or to changes in anthropogenic noise sources.     

2013年芦山地震震源区地壳介质地震波速变化的特征分析

为探究芦山M7.0级地震后5年多来,震源区龙门山断裂带西南段介质波速的变化规律,本文基于2012年4月至2018年4月共6年的连续波形数据,运用移动窗互谱与频域偏振等分析方法,结合背景噪声源的特性,对不同深度范围内的相对波速变化以及震后的恢复过程与机制进行了研究.获得的主要认识包括：（1）年尺度而言,震源区周期为1~20 s的背景噪声场相对稳定,但成分复杂、2~10 s频带内至少存在2个能量相对稳定的噪声源;不同周期噪声的能量,在月变与季节性上的变化特征差异明显.（2）获得了长时间尺度、不同频带内介质相对波速的背景变化水平,1~2 s、2~4 s的波动幅度（约为±0.04%）与季节性变化规律强于4~10 s、10~20 s的,结合与降雨量相关的地下水位模型能很好地解释其变化规律.（3）震源区的同震波速降低现象清晰,降幅约为0.08%~0.1%;空间上,波速下降最为显著的区域主要集中在龙门山断裂带两侧约70 km范围内,其中四川盆地一侧平均约为0.1%,略高于青藏高原（0.08%）一侧;在断裂带内的降速不显著.对不同子频带进行测量的结果显示,震后除10~20 s外,其余3个子频带的相对波速在震后较短时间内（约20天左右）均出现较大幅度的波速降低现象,其中4~10 s的平均降速最大（约为0.08%）,分析认为主震及大量余震的松弛效应是引起介质波速下降的主要原因.（4）震后大约1年左右,波速变化基本恢复到震前水平,且至2018年4月前未观察到大幅的波速变化现象,总体上各频带内的结果均沿零线小幅波动.     

Analysis of the primary and secondary microseisms in the wavefield of the ambient noise recorded in northern Poland

Seasonal changes of the primary and secondary microseisms were analysed in the wavefield of the ambient noise recorded during the entire 2014 at the “13 BB star” array located in northern Poland, composed of thirteen, symmetrically arranged, broadband seismic stations. To that, spectral analysis, seismic interferometry, surface scalar wind speed distribution, and beamforming were used. Spectral analysis allowed to observe that a splitting of the secondary microseism peak was present in winter and autumn, and that the primary microseism peak was visible in spring, summer and autumn. Using seismic interferometry, the long-term characteristics of the noise wavefield were recognized. The seasonal variations of the secondary microseism source were described by means of the analysis of the surface scalar wind speed for each month. The splitting of the secondary peak was attributed to the interaction of a strong wind blowing from the North Sea with a weak wind blowing from the Baltic Sea. The seasonal variations of the primary microseism peak were characterized through the frequency-domain beamforming. The peak was identified during spring, summer and autumn, when the generated wavefield was coming from the Baltic Sea. The velocity of the wavefield was evaluated within the 2.0–5.0 km/s range. The described mechanism of generation of the microseisms, based on the interaction of the nearby winds, was found to be consistent with the models reported in the literature.     

Microseismic noise in northern European Russia: Seismic recordings in the town of Vorkuta

Microseismic noise was studied in the frequency range 0.5–30 Hz in the conditions prevailing in the town of Vorkuta. A seismic noise model was developed consisting of power spectral densities of ground motion velocity separately for daytime and nighttime in different frequency ranges. The absolute noise level for frequencies of 1–5 Hz in Vorkuta varies between −140 and −150 dB in daytime and from −152 to −158 dB in nighttime, with the ranges for the 8-15 Hz noise being −140 to −155 dB in daytime and −155 to −165 dB in nighttime. Well-pronounced daily variations in noise amplitude were observed in the frequency ranges 1.5–3 Hz and 14–17 Hz. The noise amplitude varies by 7 dB over 24 hours, with the amplitude of horizontal component variation being 5 dB above that of the vertical component. The power spectral densities of ground motion velocity in the microseismic noise involve several spectral peaks, whose central frequencies did not vary by more than 0.15 Hz during the entire period of instrumental observation. The seismic background in the town of Vorkuta contains seismic events due to distant earthquakes and local industrial explosions.     

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.     

汶川地震区地壳速度相对变化的环境噪声自相关研究

2006年中国地震局地质研究所地震动力学国家重点实验室在川西地区(100°E～105°E, 26°N～32°N)布设了由297个宽频带地震台组成的密集流动地震观测台阵.本文利用川西流动地震台阵29°N以北地区的137个台站2007年1月至2008年10月的连续三分量地震环境噪声记录,研究了汶川地震震前震后地壳速度变化特征.借助川西台阵的密度优势,我们针对单台三分量噪声自相关函数和单台不同分量间的互相关函数,利用互谱移动窗技术,在0.33～1 Hz频带范围内测量了经过50天滑动平均的相关函数与长时间平均参考相关函数的走时变化率,进而求得地壳浅部速度随时间的相对变化,并得到其空间分布特征.本文结果表明,利用环境噪声自相关方法可以得到与互相关方法基本一致的同震速度变化分布图像,同震速度变化分布与同震体应变具有明显相关性.具有更高空间分辨率的噪声自相关研究进一步发现,在鲜水河断裂和龙门山断裂交汇区存在同震速度增加区,这个区域与同震库伦应力变化和地表形变观测预测的周边断层库伦应力增加区一致.我们的研究发现,同震库伦应力增加效应的持续时间大约为2个月左右,在此之后,区域应力场逐渐恢复为普遍的应力下降.     

Coal Mining Induced Seismicity in the Ruhr Area, Germany

Over the last 25 years mining-induced seismicity in the Ruhr area has continuously been monitored by the Ruhr-University Bochum. About 1,000 seismic events with local magnitudes between 0.7 ≤ M _L ≤ 3.3 are located every year. For example, 1,336 events were located in 2006. General characteristics of induced seismicity in the entire Ruhr area are spatial and temporal correlation with mining activity and a nearly constant energy release per unit time. This suggests that induced stresses are released rapidly by many small events. The magnitude–frequency distribution follows a Gutenberg–Richter relation which is a result from combining distributions of single longwalls that themselves show large variability. A high b-value of about 2 was found indicating a lack of large magnitude events. Local analyses of single longwalls indicate that various factors such as local geology and mine layout lead to significant differences in seismicity. Stress redistribution acts very locally since differences on a small scale of some hundreds of meters are observed. A regional relation between seismic moment M ₀ and local magnitude M _L was derived. The magnitude–frequency distribution of a single longwall in Hamm was studied in detail and shows a maximum at M _L = 1.4 corresponding to an estimated characteristic source area of about 2,200 m². Sandstone layers in the hanging or foot wall of the active longwall might fail in these characteristic events. Source mechanisms can mostly be explained by shear failure of two different types above and below the longwall. Fault plane solutions of typical events are consistent with steeply dipping fracture planes parallel to the longwall face and nearly vertical dislocation in direction towards the goaf. We also derive an empirical relation for the decay of ground velocity with epicenter distance and compare maximum observed ground velocity to local magnitude. This is of considerable public interest because about 30 events larger than M _L ≥ 1.2 are felt each month by people living in the mining regions. Our relations, for example, indicate that an event in Hamm with a peak ground velocity of 6 mm/s which corresponds to a local magnitude M _L between 1.7 and 2.3 is likely to be felt within about 2.3 km radius from the event.     

利用噪声研究地下介质速度结构及其变化在中国大陆的应用

背景噪声层析成像技术已被广泛应用于地壳和上地幔速度结构的研究,这种方法不依靠地震的发生和人工源爆破,只需记录连续的噪声信号而无需产生信号,因为噪声穿过地下介质时会携带信息,然后通过利用台站记录到的连续背景噪声数据进行互相关计算和叠加,即可得到台站间的经验格林函数,从而获取对地下结构的认识。这种方法已经很好地应用于中国的东北地区、华北克拉通、青藏高原以及华南地区,并成功地揭示了这些地区地壳与上地幔顶部的速度结构。此外近年来,一些学者还利用噪声互相关技术研究地下介质地震波速度随时间的变化,通过对汶川大地震前后连续噪声记录的研究发现,大震发生后呈现同震波速降低和震后波速逐渐恢复的特点,这表明可以通过观测地震波特性的变化来监测地下应力的变化,从而为大震的预测预防工作提供科学依据。本文主要综述了近些年来背景噪声技术及其在中国大陆地区的应用。     

利用背景噪声互相关研究汶川地震震源区地震波速度变化

利用2007年3月至2009年3月四川数字地震台网的宽频带连续波形资料,通过计算地震背景噪声互相关提取台站对间的经验格林函数,在0.1～0.5 Hz频带下测量每天经验格林函数与参考经验格林函数的走时偏移,进而得到各台站对在该时段内的相对地震波速度变化.结果表明,2008年5月12日汶川M_s8.0级地震造成了震源区地震波速度的急剧降低,最大降幅达0.4%;大致以安县为界,余震带西南部地区在汶川主震后波速降即达到最大值,而东北部地区的最大波速降一般出现在主震后的1~4个月,相对地震波速度变化的这种分段特性与地震序列的时空分布特征有较好的对应关系;在震源区外围的四川盆地也观测到了震后波速降低,而川西高原内部则没有出现显著的波速变化.进一步的分析和计算结果表明主震的静态应力变化和强地面运动引起的地表破坏都不能很好地解释震后波速的急剧降低,地震导致的断层区内部结构破坏和周边介质应力状态改变可能是波速变化的主要原因.     

Semiannual velocity variations around the 2008 M w 7.9 Wenchuan Earthquake fault zone revealed by ambient noise and ACROSS active source data

We continuously monitor the long-term seismic velocity variation of one of the major ruptured faults of the devastating 2008 M _w 7.9 Wenchuan earthquake in China from July 2009 to January 2012, jointly using accurately controlled routinely operated signal system active source and seismic noise-based monitoring technique. Our measurements show that the temporal velocity change is not homogeneous and highly localized in the damaged fault zone and the adjacent areas. Velocity variations from the active and passive methods are quite consistent, which both are characterized by ±0.2 % seasonal variation, with peak and trough at winter and summer, respectively. The periodic velocity variation within fault zone exhibits remarkably positive correlation with barometric pressure with stress sensitivity in the order of 10^?6 Pa^?1, suggesting that the plausible mechanism might be the crack density variation of the shallow subsurface medium of the damaged fault zone in response to the cyclic barometric pressure loading.     

利用时移层析成像方法揭示与2014年云南鲁甸MS6.5地震有关的P波速度变化

基于双差层析成像的时移层析成像方法能够获得不同时间段之间的地下介质的速度变化,且不同的数据分布和质量对结果影响较小.在本研究中我们使用云南地震区域台网所记录到的P波绝对到时及相对到时数据,利用时移层析成像方法得到了2014年鲁甸M S6.5地震发生前后震源区高精度的P波速度变化的时空分布.结果表明:2014年鲁甸地震发生后,震源区同震期P波速度下降,但没有下降至最大,而是在震后1—4个月内下降至最大,接着P波速度开始上升,震源区开始愈合,愈合过程从浅层逐渐发展至深层,逐渐恢复至震前水平.同时发现在空间上P波速度变化与余震分布变化相一致,鲁甸地震同震应力变化与速度变化之间也有较好的一致性,所以认为地震引起的应力变化是造成鲁甸地震速度变化的一个重要原因,余震的动态和静态应力造成了震源区介质物理属性的改变从而影响了P波速度的变化.     

Preliminary analysis on the source properties and seismogenic structure of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>7.0 Jiuzhaigou earthquake

At GMT time 13:19, August 8, 2017, an M_s7.0 earthquake struck the Jiuzhaigou region in Sichuan Province, China, causing severe damages and casualties. To investigate the source properties, seismogenic structures, and seismic hazards, we systematically analyzed the tectonic environment, crustal velocity structure in the source region, source parameters and rupture process, Coulomb failure stress changes, and 3-D features of the rupture plane of the Jiuzhaigou earthquake. Our results indicate the following: (1) The Jiuzhaigou earthquake occurred on an unmarked fault belonging to the transition zone of the east Kunlun fault system and is located northwest of the Huya fault. (2) Both the mainshock and aftershock rupture zones are located in a region where crustal seismic velocity changes dramatically. Southeast to the source region, shear wave velocity at the middle to lower crust is significantly low, but it rapidly increases northeastward and lies close to the background velocity across the rupture fault. (3) The aftershock zone is narrow and distributes along the northwest-southeast trend, and most aftershocks occur within a depth range of 5–20 km. (4) The focal mechanism of the Jiuzhaigou earthquake indicates a left-lateral strike-slip fault, with strike, dip, and rake angles of 152°, 74° and 8°, respectively. The hypocenter depth measures 20 km, whereas the centroid depth is about 6 km. The co-seismic rupture mainly concentrates at depths of 3–13 km, with a moment magnitude (M_w) of 6.5. (5) The co-seismic rupture also strengthens the Coulomb failure stress at the two ends of the rupture fault and the east segment of the Tazang fault. Aftershocks relocation results together with geological surveys indicate that the causative fault is a near vertical fault with notable spatial variations: dip angle varies within 66°–89° from northwest to southeast and the average dip angle measures ~84°. The results of this work are of fundamental importance for further studies on the source characteristics, tectonic environment, and seismic hazard evaluation of the Jiuzhaigou earthquake.     

Imaging heterogeneity of the crust adjacent to the Dead Sea fault using ambient seismic noise tomography

For the purpose of studying the Earth’s crust by means of tomography, we investigated cross-correlation functions emerging from long-term observations of propagating ambient seismic noise at pairs of broadband stations in Israel and Jordan. The data was provided by the eight permanent broadband stations of the Israel Seismic Network evenly distributed over Israel and the 30 stations of the DESERT2000 experiment distributed across the Arava Fault (South of the Dead Sea basin). To eliminate the influence of earthquakes and explosions, we have applied the methodology of Bensen et al. (Geophys J Int 169:1239–1260, 2007), including bandpass filtering and amplitude normalization in time and frequency domain. The cross-correlation functions estimated from continuous recordings of a few months were used to extract Rayleigh waves group velocity dispersion curves using automatic version of the frequency–time analysis procedure. Subsequently, these curves have been converted into the Rayleigh wave group velocity maps in the period range 5–20 s and S waves velocity maps in the depth range 5–15 km. In these maps, four velocity anomalies are prominent. Two of them are outlined by the previous reflection-refraction profiles and body wave tomography studies, i.e. a low velocity anomaly corresponds to the area of the extremely deep (down to 14 km) sedimentary infill in the Southern Dead Sea Basin and a high velocity anomaly in the Southern Jordan corresponds to the area of the Precambrian crystalline rocks of the Nubian Shield on the flanks of the Red Sea. The two other anomalies have not been reported before - the high velocity zone close to the Beersheba city and the low velocity anomaly in the region of Samaria-Carmel mountains - Southern Galilee. They have relatively low resolution and hence need further investigations for approving and contouring. The highest contrast between the average Rayleigh wave group velocity (2.7 km/s) and the anomalies is 10–13 %, comparable, however, to the level of noise in the data. The results have been verified by modeling the revealed anomalies which showed that all the four zones mentioned above could be detected by the tomography study.     

Comparison of source parameters estimated in the frequency and time domains for seismic events at the Rudna copper mine,Poland

Source parameters estimated in the frequency domain for 100 selected seismic events from the Rudna copper mine, with moment magnitude ranging from 1.4 to 3.6, were collected to study their scaling relations and to compare them with the parameters estimated in the time domain. The apparent stress and static stress drop, corrected for the limited bandwidth recording, increase slightly in a similar manner with increasing seismic moment. The ratio of apparent stress to static stress drop, a measure of radiation efficiency, is practically constant and its mean value is close to 0.1. For 37 seismic events, with moment magnitude between 1.9 and 3.4, source parameters were estimated in the time domain from relative source time functions, that displayed unilateral rupture propagation, and their rupture velocity could be estimated. It ranges from 0.23 to 0.80 of shear wave velocity and is almost independent of seismic moment. The fault length, estimated from the average source pulse width and rupture velocity, is clearly dependent on seismic moment and is smaller than the source radius estimated from the corner frequency on the average by about 25 percent. There is no correlation between the values of static stress drop estimated in the frequency and time domains, but the time domain stress drop is in general similar to that estimated in the frequency domain. The apparent stress increases with increasing rupture velocity, and the ratio of apparent stress to static stress drop seems also to depend on rupture velocity.     

Ambient noise levels and detection threshold in Norway

Ambient seismic noise is caused by a number of sources in specific frequency bands. The quantification of ambient noise makes it possible to evaluate station and network performance. We evaluate noise levels in Norway from the 2013 data set of the Norwegian National Seismic Network as well as two temporary deployments. Apart from the station performance, we studied the geographical and temporal variations, and developed a local noise model for Norway. The microseism peaks related to the ocean are significant in Norway. We, therefore, investigated the relationship between oceanic weather conditions and noise levels. We find a correlation of low-frequency noise (0.125–0.25 Hz) with wave heights up to 900 km offshore. High (2–10 Hz) and intermediate (0.5–5 Hz) frequency noise correlates only up to 450 km offshore with wave heights. From a geographic perspective, stations in southern Norway show lower noise levels for low frequencies due to a larger distance to the dominant noise sources in the North Atlantic. Finally, we studied the influence of high-frequency noise levels on earthquake detectability and found that a noise level increase of 10 dB decreases the detectability by 0.5 magnitude units. This method provides a practical way to consider noise variations in detection maps.     

Retrieval of super‐virtual refraction by cross‐correlation

Topography and severe variations of near‐surface layers lead to travel‐time perturbations for the events in seismic exploration. Usually, these perturbations could be estimated and eliminated by refraction technology. The virtual refraction method is a relatively new technique for retrieval of refraction information from seismic records contaminated by noise. Based on the virtual refraction, this paper proposes super‐virtual refraction interferometry by cross‐correlation to retrieve refraction wavefields by summing the cross‐correlation of raw refraction wavefields and virtual refraction wavefields over all receivers located outside the retrieved source and receiver pair. This method can enhance refraction signal gradually as the source–receiver offset decreases. For further enhancement of refracted waves, a scheme of hybrid virtual refraction wavefields is applied by stacking of correlation‐type and convolution‐type super‐virtual refractions. Our new method does not need any information about the near‐surface velocity model, which can solve the problem of directly unmeasured virtual refraction energy from the virtual source at the surface, and extend the acquisition aperture to its maximum extent in raw seismic records. It can also reduce random noise influence in raw seismic records effectively and improve refracted waves’ signal‐to‐noise ratio by a factor proportional to the square root of the number of receivers positioned at stationary‐phase points, based on the improvement of virtual refraction's signal‐to‐noise ratio. Using results from synthetic and field data, we show that our new method is effective to retrieve refraction information from raw seismic records and improve the accuracy of first‐arrival picks.     

2014年鲁甸MS6.5地震前后介质波速变化特征分析

本文利用研究区内13个地震台站2014年1月至12月的连续波形数据, 采用背景噪声互相关方法, 提取了2014年8月3日鲁甸M_S6.5地震的同震波速变化。 结果表明, 鲁甸地震对介质波速的影响具有明显的空间差异性, 较大幅度的同震波速变化主要集中于则木河断裂—小江断裂及莲峰断裂区域。 其中, 莲峰断裂SW段和则木河断裂NW段区域的波速出现了较为明显的同震降低, 则木河断裂和小江断裂交界处则出现了明显的波速升高。 通过对比研究区内介质波速的同震变化与鲁甸地震对周边主要活动断裂应力积累的影响后发现, 波速变化与应力变化在空间分布上具有较高的一致性, 且两者的变化幅度也成正相关关系。 由此可以看出, 鲁甸地震造成的应力变化可能是地下介质波速变化的主要原因。 采用相同的方法对2014年4月5日永善M_S5.1地震研究后发现, 永善地震后鲁甸地震震源区及其附近区域波速明显升高。     

The effect of baric disturbances in the atmosphere on microseismic processes in the crust

The effect of baric variations of different origins on characteristics of seismic noise is analyzed in the frequency range 0.03–20 Hz. Long period variations in atmospheric pressure caused by cyclones, whose period T ranges from half a day to a few days, are shown to increase the microseismic background amplitude by two to four times in the frequency range 0.03–1 Hz (the coefficient of linear correlation between time variations in the amplitude and atmospheric pressure is K = 0.65 at a significance level of r = 0.95). Short-period baric variations with T ～ 5–30 min associated with the passage of cold fronts lead to a tenfold increase in the microseismic background amplitude in the frequency range 4–8 Hz (K = 0.67 at r = 0.95). In this case, disturbances of seismic background are recorded for 20–60 min after the passage of an atmospheric front and display an exponential drop in the amplitude. In distinction to cyclones, an atmospheric front increases the number of impulsive microseismic events of the resonance type.     

Linking ground motion measurements and macroseismic observations in France: a case study based on accelerometric and macroseismic databases

Comparison between accelerometric and macroseismic observations is made for three M _w?=?4.5 earthquakes, which occurred in north-eastern France and south-western Germany in 2003 and 2004. Scalar and spectral instrumental parameters are processed from the accelerometric data recorded by nine accelerometric stations located between 29 and 180 km from the epicentres. Macroseismic data are based on French Internet reports. In addition to the single questionnaire intensity, analysis of the internal correlation between the encoded answers highlights four predominant fields of questions bearing different physical meanings: (1) “vibratory motions of small objects”, (2) “displacement and fall of objects”, (3) “acoustic noise” and (4) “personal feelings”. Best correlations between macroseismic and instrumental observations are obtained when the macroseismic parameters are averaged over 10-km-radius circles around each station. Macroseismic intensities predicted by published peak ground velocity (PGV)–intensity relationships agree with our observed intensities, contrary to those based on peak ground acceleration (PGA). Correlation between the macroseismic and instrumental data for intensities between II and V (EMS-98) is better for PGV than for PGA. Correlation with the response spectra exhibits clear frequency dependence for all macroseismic parameters. Horizontal and vertical components are significantly correlated with the macroseismic parameters between 1 and 10 Hz, a range corresponding to both natural frequencies of most buildings and high energy content in the seismic ground motion. Between 10 and 25 Hz, a clear lack of correlation between macroseismic and instrumental observations exists. It could be due to a combination of the decrease in the energy signal above 10 Hz, a high level of anthropogenic noise and an increase in variability in soil conditions. Above 25 Hz, the correlation coefficients between the acceleration response spectra and the macroseismic parameters are close to the PGA correlation level.     

浙江珊溪水库地震震源参数特征研究

利用2014年1月-2015年12月浙江省区域数字地震台网记录的1 360次珊溪水库M_L0.5-4.4地震资料,采用Brune模式,将台站获取的速度记录进行几何扩散校正、介质衰减校正及仪器校正后得到速度谱,对速度谱在频率域进行傅里叶积分得到震源位移谱,之后使用遗传算法计算得到拐角频率和震源位移谱零频极限值,进而计算出该区的地震矩为1010-1014 N·m,震源破裂半径为33-550 m,拐角频率为2.4-39.7 Hz,地震应力降为0.04-6.74 MPa,视应力为0-2.75 MPa.在此基础上,对各参数特征及其关系进行系统性分析,结果表明:各震源参数之间表现出一定的对数线性或半对数线性关系;空间上,应力降和拐角频率在库区西北段较高,东南段较低;应力降与地震矩在空间上呈反向关系,拐角频率与破裂半径在空间上呈负相关.