强震动记录自适应基线校正方法:汶川地震应用

基线校正是强震动记录数据处理非常关键的环节之一。现有方法大都依赖主观经验,得到的永久位移值存在较大的不确定性,不能实现自动化处理。为此,提出一种自适应的强震动记录基线校正方法,以汶川地震中汶川卧龙台站(051WCW)获取的强震动记录为例叙述新方法的每个环节,并与已有方法进行比较,验证了新方法的计算准确性。最后,将新方法运用于汶川地震主震中,得到了37个强震动台站的永久位移场。通过与Wang等的研究结果比较,发现除3个台站由于数据不全无法比较外,新方法给出的24个台站的最大永久位移无论是数值还是对应方位角都与其一致。表明新方法具有一定的可靠性,同时通过对比也发现新方法存在一定的局限性,还有待进一步改进。     

2017年九寨沟MS7.0地震场地特征分析

场地条件对地震动具有较大影响, 研究不同场地条件下的地震动特征对地震动的校正具有重要意义。 本文以九寨沟地震为例, 收集了66个台站的198条三分量强震观测记录和SRTM(Shuttle Radar Topography Mission)公里格网的DEM(Digital Elevation Model)数据, 从多个角度对场地特征进行分析。 首先根据坡度法使用DEM数据对九寨沟特征进行了场地分类, 然后讨论了不同场地类型下的加速度时程、 加速度反应谱、 地震动衰减等地震动特征。 研究表明, 缺少实测资料的情况下, 地形坡度可以作为Vs30(地表以下30 m范围的平均剪切波速度)的一种替代指标, 利用坡度法可以较快速地对场地进行分类; 利用强震观测记录能够从多角度对台站的场地特征进行分析, 不同场地类型对地震动影响不同, 其中, 土层对地震动具有明显的放大效应。 该结论可以为地震动结果的校正提供依据。     

Strong motion observations and recordings from the great Wenchuan Earthquake

The National Strong Motion Observation Network System （NSMONS） of China is briefly introduced in this paper. The NSMONS consists of permanent free-field stations, special observation arrays, mobile observatories and a network management system. During the Wenchuan Earthquake, over 1,400 components of acceleration records were obtained from 460 permanent free-field stations and three arrays for topographical effect and structural response observation in the network system from the main shock, and over 20,000 components of acceleration records from strong aftershocks occurred before August 1, 2008 were also obtained by permanent free-field stations of the NSMONS and 59 mobile instruments quickly deployed after the main shock. The strong motion recordings from the main shock and strong aftershocks are summarized in this paper. In the ground motion recordings, there are over 560 components with peak ground acceleration （PGA） over 10 Gal, the largest being 957.7 Gal. The largest PGA recorded during the aftershock exceeds 300 Gal.     

利用固定台站分析长江激发气枪信号特征

“地学长江计划”安徽实验是以气枪震源为核心的大型主动源探测实验。通过在长江安徽段20个固定点定点激发气枪震源,结合109个固定台站、11条流动测线组成的观测网络,首次利用主动源实现了对长江流域安徽段约6万km2面积的三维地下结构探测。本文利用固定台站对长江激发气枪信号进行了分析,结果表明,长江中气枪信号激发效果良好,固定台记录中气枪信号可识别的最远距离达300km。对气枪信号绝对振幅的研究结果表明：① 50km处的气枪信号约为10nm量级,200km处的气枪信号小于1 nm;② 气枪信号强度的空间分布存在一定的方位各向异性,可能与长江的几何形状有关;③ 台站背景噪声对于提取气枪信号至关重要,高质量的固定台网为识别nm量级气枪信号提供了可能。     

基于强震动记录确定的场地卓越周期

本文介绍了3种根据场地强震动记录获取场地卓越周期的方法: ① 地表记录的傅里叶谱分析法; ② 地表水平/垂直傅里叶谱比法; ③ 地表/地下傅里叶谱比法. 基于日本强震动观测台网KiK-net中两个基岩台站和两个Ⅲ类场地台站获取的数百条强震动记录, 分别使用上述3种方法确定场地的卓越周期, 并对比分析各种方法的优缺点及其适用情况. 结果表明: 对于基岩场地, 由于记录信息复杂, 局部场地条件对地震动影响较小, 地震动自身特性突出, 导致3种方法所得到的结果均比较分散; 对于土层场地, 场地条件影响比较显著, 3种方法基本都可以得到一个较为准确的数值; 但对于某些场地, 地表水平/垂直傅里叶谱比法所得结果存在不确定性, 相比之下, 地表/地下傅里叶谱比法则能给出一个更为准确的场地卓越周期值.      

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.     

The implications of regional variations in rainfall for reconstructing rainfall patterns using tree rings

In Australia, multidecadal periods of floods and droughts have major economic consequences. Due to the short duration of Australian instrumental precipitation records, it is difficult to determine the patterns of these multidecadal periods. Proxy records can be used to create long‐term rainfall reconstructions for regions that are lacking instrumental data. However, the spatial extent over which single‐site proxy records can be applied is poorly understood. Southeast Queensland (SEQ) is an area where tree rings can be used to reconstruct long‐term rainfall patterns, but their regional representation is unknown. In this study, the spatial variability in rainfall across SEQ is investigated from 1908 to 2007 using 140 instrumental rainfall stations. Pearson correlation analysis between stations is used to create groups at the r = 0.80, 0.85, and 0.90 correlation levels, and then annual deviations from the mean are determined. These patterns indicate that rainfall is not uniform across SEQ but can be broken into 2 main spatially consistent groups. Each of these groups is broken down into several subgroups with higher correlation levels. Long‐term streamflow records are found to be correlated to rainfall patterns local to the streamflow stations, indicating that analysis of extreme events should consider spatial precipitation variability. Finally, the only currently available proxy rainfall reconstruction for the region, a 140‐year Toona ciliata tree ring width record from Lamington National Park, is compared to rainfall groups at different correlation levels across all of SEQ. The correlation between the reconstruction and the rainfall station groupings is best for the groups within which the tree‐ring record is spatially located, and this correlation improves as rainfall group correlation increases. Correlation is nearly nonexistent for groupings located at a distance from the tree‐ring site. These results demonstrate the importance of assessing the spatial variability of precipitation so that the spatial applicability of proxy records can be assessed.     

基于实测记录分析场地因素对海底地震动特性影响

目前基于海底实测记录的分析发现海底与陆地地震动特性存在明显差异,但难以进一步确定海底地震动特性的影响因素。在以往研究的基础上对比同次地震中相邻海底台站间地震动特性的差异,并分析造成差异的原因。以日本K-NET地震台网中6个海底强震台站及其相邻不同场地条件陆地台站监测的8次强震记录为研究对象,通过分析强震记录的峰值加速度、水平放大系数谱、竖向与水平反应谱的比谱等,对比分析不同海底台站地震动的特性,以及海底台站与相邻不同场地条件陆地台站地震动特性的差别。研究发现：（1）海底与陆地竖向地震动存在明显差异;（2）不同海底台站间地震动特性亦存在较大差异和明显的规律性,海底场地条件、地形等场地因素对海底地震动特性的影响较大;（3）海底水平向地震动反应谱的特征周期较大,谱特性介于陆地中硬土与软土场地间。     

基于谱比法的江苏强震动台站场地特征分析

选取江苏省强震动台网近年来获得的强震动记录,基于水平竖向谱比法(HVSR),计算了22个强震动台站的速度反应谱平均谱比曲线,采用谱比形状一致法进行场地分类,并结合典型台站场地的工程地质资料进行对比分析。研究结果表明,平均谱比曲线的卓越周期、场地放大等特征与台站场地地质资料相吻合,与传统场地分类方法相比,基于谱比法的场地特征分类方法能够突出不同场地的地震动反应特征,较好地体现了台站的场地反应特征。     

一种实时场地校正方法及其在地震预警中的应用

地震预警系统对地震数据处理的实时性要求极高。其系统数据来源除布设在基岩的测震台站外,还有大量非基岩场地的强震动台站和地震烈度仪台站,其场地影响不容忽视。为了考虑震级估算和地震动场预测中的场地影响,需实时对各种场地条件下的地震波形进行校正。目前处理一般使用某个标量来表征场地放大效应。本文采用一种实时的场地校正方法,首先计算目标场地与参考场地的谱比,然后通过最小二乘法、双线性变换将谱比转化为因果递归的无限冲激响应(IIR)时域滤波器,之后可以应用该滤波器进行实时场地校正。该方法考虑了场地放大系数的频率依赖性,相比于标量校正提高了准确度。应用我国四川和日本部分强震动台站记录,验证并讨论了这种实时场地校正方法在地震预警中的应用效果。     

一种利用单台强震记录计算任意场地三分量地面运动的方法

本文利用近场理论格林函数反褶积的方法计算土层的地震反应.解答中包含了产生地面运动的各种反射、透射和转换波震相,也计入了土层与基底介质的相互作用.将计算的土层响应与已获得的近场加速度记录褶积,可算出不同土层场地上垂直向和水平向三分量的地面运动.文中利用1982年卢龙地震的近场加速度记录,检验了方法的可行性.     

一种利用单台强震记录计算任意场地三分量地面运动的方法

本文利用近场理论格林函数反褶积的方法计算土层的地震反应.解答中包含了产生地面运动的各种反射、透射和转换波震相,也计入了土层与基底介质的相互作用.将计算的土层响应与已获得的近场加速度记录褶积,可算出不同土层场地上垂直向和水平向三分量的地面运动.文中利用1982年卢龙地震的近场加速度记录,检验了方法的可行性.     

Identification of accelerometric stations in ITACA with distinctive features in their seismic response

In this work, we analyze the records of the Italian strong motion database (ITACA, ) with the aim of identifying stations affected by site effects that are not captured by standard seismic classification schemes. In particular, we consider four different site classifications, two of them based on geological/geophysical characteristics and two driven by data. For each classification we develop a ground motion prediction equation using a random effect approach to isolate the between-station and within-station distribution of errors. The site coefficients obtained for the different classes confirm that site amplification effects are significant for both the horizontal and vertical components. The between-station error normalized to the standard deviation of the between-station error distribution is then used to identify stations characterized by large errors, attributable to site effects not accounted for by the classification schemes. The results show that large errors can affect the predictions when the site effects are not uniquely related to the reduction of the seismic impedance in the uppermost layers. For example, amplifications of ground motion over the long period range are observed for stations installed within alluvial closed-shape basins, as consequence of locally generated surface waves. For these stations, classifications based on the horizontal to vertical response spectra ratio are not reliable, since amplifications are also affecting the vertical component. Another interesting feature which emerges from the analysis is the significant de-amplification of short period spectral ordinates that seems to be related to stations typically set in at the foundation level of massive structures. To increase the usefulness of the data set, the most important distinctive features of the strong motion stations are documented in the ITACA database reports containing the instrument information and the available geological-geotechnical data.     

Characterization of site conditions (soil class,V<Subscript>S30</Subscript>, velocity profiles) for 33 stations from the French permanent accelerometric network (RAP) using surface-wave methods

Data provided by accelerometric networks are important for seismic hazard assessment. The correct use of accelerometric signals is conditioned by the station site metadata quality (i.e., soil class, V_S30, velocity profiles, and other relevant information that can help to quantify site effects). In France, the permanent accelerometric network consists of about 150 stations. Thirty-three of these stations in the southern half of France have been characterized, using surface-wave-based methods that allow derivation of velocity profiles from dispersion curves of surface waves. The computation of dispersion curves and their subsequent inversion in terms of shear-wave velocity profiles has allowed estimation of V_S30 values and designation of soil classes, which include the corresponding uncertainties. From a methodological point of view, this survey leads to the following recommendations: (1) perform both active (multi-analysis surface waves) and passive (ambient vibration arrays) measurements to derive dispersion curves in a broadband frequency range; (2) perform active acquisitions for both vertical (Rayleigh wave) and horizontal (Love wave) polarities. Even when the logistic contexts are sometimes difficult, the use of surface-wave-based methods is suitable for station-site characterization, even on rock sites. In comparison with previous studies that have mainly estimated V_S30 indirectly, the new values here are globally lower, but the EC8-A class sites remain numerous. However, even on rock sites, high frequency amplifications may affect accelerometric records, due to the shallow relatively softer layers.     

Average crustal travel times in Sweden re-examined

A long-range Scandinavian seismic profile of August 1979 is subjected to a preliminary evaluation by means of the records of the permanent Swedish network of stations. The resulting travel times confirm those earlier adopted in the calculation of Swedish earthquakes for the period 1951–1976.     

Towards an Automatic Monitoring System of Infrasonic Events at Mt. Etna: Strategies for Source Location and Modeling

Active volcanoes characterized by open conduit conditions generate sonic and infrasonic signals, whose investigation provides useful information for both monitoring purposes and studying the dynamics of explosive processes. In this work, we discuss the automatic procedures implemented for a real-time application to the data acquired by a permanent network of five infrasound stations running at Mt. Etna volcano. The infrasound signals at Mt. Etna consist in amplitude transients, called infrasound events. The adopted procedure uses a multi-algorithm approach for event detection, counting, characterization and location. It is designed for an efficient and accurate processing of infrasound records provided by single-site and array stations. Moreover, the source mechanism of these events can be investigated off-line or in near real-time by using three different models: (1) Strombolian bubble; (2) resonating conduit and (3) Helmholtz resonator. The infrasound waveforms allow us to choose the most suitable model, to get quantitative information about the source and to follow the time evolution of the source parameters.     

由多台定点长水管倾斜资料构建形变场时空演化图像的初步尝试及应用

尝试利用多台定点地倾斜形变资料构建形变场时空动态演化图形的方法,对中国大陆47个定点形变台水管倾斜仪资料进行分析处理,构建了地倾斜形变场时空演化分布图像.以2001年11月14日昆仑山口西8.1级地震为例,研究了地震前后水管倾斜资料的时空演化特征.结果显示在地震前后较多台站水管倾斜仪的日均值都出现了趋势上的较大变化,东西向日均值、南北向日均值和南北向潮汐因子的空间分布也都出现了大范围的变化,而东西向M2波的潮汐因子没有出现较明显的变化.随着定点形变台站的建设和资料的积累,该方法可以构建每月或者每周形变场的变化,有望用于地震预报工作.     

H/V spectral ratio analysis and Rayleigh modelization in Eastern Thuringia,Germany

We use records from the East Thuringian Seismic Network (OTSN, Ostthüringer Seismisches Netzwerk) to characterize the site response for each station, and to analyze the scope and limits of the Rayleigh modeling for H/V spectral ratio. The stations considered in this work can be classified by their seismic response as hard rock sites or as sites with some site effect. From these results we propose velocity models based on Rayleigh modeling (theoretical Rayleigh wave ellipticity).Our results show that for locations affected by site effects the H/V spectral ratio can be modeled by the theoretical ellipticity of layered velocity models. For hard rock sites the spectral ratio is rather flat and the modeling with the theoretical ellipticity was not very clear. This may be explained by the fact that for hard rock sites the conditions for a clear fundamental frequency associated with S-wave resonance, and therefore with Rayleigh wave ellipticity, are not fulfilled.     

2017年8月8日四川九寨沟7.0级地震强震记录及特征分析

2017年8月8日四川九寨沟县发生7.0级地震,中国数字强震动台网布设在四川、甘肃、陕西、宁夏的66个强震台获得主震加速度记录。本文首先对198条三分向强震记录进行常规处理,计算出近场强震记录的加速度峰值随震中距的分布情况;再根据2个典型台站的加速度时程记录,通过计算其加速度反应谱并与设计反应谱比较,分析本次地震的基本特征;然后将实际观测数据与意大利新一代地震动衰减公式对比,分析峰值加速度（PGA）及谱加速度的衰减关系;最后结合已有的工程场地钻探资料,采用H/V谱比法对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给出的震级相同;分析还表明背景噪声是影响台站信号检测能力的主要因素之一.