渭河盆地中土层场地对地震动的放大作用

土层场地对地震动的影响较大,建(构)筑物的选址及其抗震设防必须考虑土层场地的放大作用,以避免或减轻其震害.汶川地震中,布设在渭河盆地中的数字强震动台网共有27个台站(包括2个基岩台站和25个土层台站)获得良好的主震加速度时程.利用这些加速度时程,选择汤峪台做为参考场地,基于考虑几何衰减的传统谱比法分析研究了25个土层场...     

基于递归高通滤波的经验模态分解及其在地震信号分析中的应用

将地震信号分解成包含频谱互不重叠的单主周期的分量有利于地震信号的分析.分析了经验模态分解(EMD)中模态混叠的内在原因和已有的解决方法,梳理了解决模态混叠的思路框架,进而提出了一种新的基于输入递归高通滤波的EMD算法.首先用递归高通滤波器将信号预分解成频率由高到低的多个分量,实现信号的等价带通滤波,再用EMD对各带通分量按频率高低逐级递归筛分,获得完备的经验模态分量.通过合成信号和地震信号的仿真实验表明,该算法较好地克服了模态混叠,获得了频谱互不重叠的单主周期分量,并成功用于震相分离和分析,为地震信号分析提供了一种新思路.     

曹妃甸地震台网天然地震与人工爆破时频分析

搜集2014—2015年曹妃甸测震台网记录的各类爆破与天然地震事件,基于Matlab软件,编制程序实现短时傅里叶变换,对比分析核爆炸、化学爆炸、矿爆炸与天然地震的特征,分析频率随时间的变化规律。研究结果表明,天然地震与爆破在频率范围、频率分布和时频聚集性等方面具有各自特点,主要表现在:天然地震的时频图"多峰"且复杂,频率分布范围更广,频率成分集中在S波段;人工爆破的时频图"少峰"且简单,频率范围较窄,频率成分集中在P波段。本文的研究成果能够为今后区分天然地震和爆破提供依据。     

基于小波包分解和模糊聚类的网格结构损伤诊断方法

本文针对网格结构的特点,利用其在快速正弦扫频激励下的动力响应,基于小波包分解频带能量分析和模糊聚类法,提出了一种适合于网格结构的损伤诊断方法。利用两个传感器损伤信息的联合诊断,有效地解决了使用一个传感器无法检测出对称损伤的难题。利用有限元模拟,应用上述方法对一单层球面网壳结构进行了损伤诊断。结果表明,对网壳杆件较小程度的损伤,子频带能量成分变化敏感,用其构造诊断向量和建立损伤样本库,并进行模糊聚类和目标识别,基本可以正确地识别出网壳杆件的损伤位置和程度。同时该方法测试时间短,使用设备少,有可靠的理论基础,具有一定的工程应用参考价值。     

Analysis of Sacks Body Strain Interference at Beida No.200 Station in Changping and Earthquake Case Study

The relationship between Sacks body strain deformation at Beida No. 200 station in Changping and tidal solids,atmospheric pressure and water level is analyzed in this paper. Sacks body strain deformation data before the M_S8. 0 Wenchuan earthquake is studied based on the analysis of the interference. The short-impending anomaly of the body strain deformation is considered to be reliable. The anomaly characteristics conclude:( 1) The trend anomaly as extensional change of the body strain deformations on a quasi 1 year time scale before the Wenchuan earthquake was recorded, and the accumulative amount was about 4000 × 10~(-9). Correspondingly,the short-term precursor of earthquake was manifested as an extensional abrupt change.( 2) The extensional intermittent anomalous abrupt change was recorded by body strainmeters between March1 and May 7 in 2008.( 3) Four compressional abrupt changes were recorded in the intermittent distortions recorded between April 13 and May 11.( 4) High frequency components were increased in the distortion process in May 1 to 3,5,7,and 9 to 12,caused by slow earthquakes before the Wenchuan earthquake according to wavelet analysis. The abnormal phenomena are summarized and the mechanics discussed in this paper. Strain solid tide distortions in body strain observations,the continuous repeated extensional and compressional abrupt changes accompanying these distortions,and the increase of high frequency components can be regarded as the index of short term and impending earthquake prediction,based on analysis of interference factors such as air pressure and water level.     

Evaluation of simulated ground motions for predicting elastic response of long period structures and inelastic response of structures

The response of linear elastic and non-linear hysteretic systems having a single degree of freedom to recorded and simulated ground motions is studied. The objective is to evaluate whether the commonly used simulated motions are appropriate for predicting inelastic response of structures and elastic response of long period structures. Eight simulated motions were generated to model properties of horizontal ground motions recorded during four earthquakes. The simulated motions are sample functions of a stationary, Gaussian white noise process, multiplied by a temporal intensity function and passed through a linear single-degree-of-freedom filter. Two versions, corresponding to parabolic and ‘standard’ base line corrections (BLC), of each of the simulated and recorded accelerograms were considered. The following general conclusions are deduced. Simulated ground motions should be subjected to the standard BLC, because it results in more reliable ground velocities and displacements, which in turn would lead to more reliable predictions of response of long period structures. Furthermore, the spectral density of the underlying random process, from which the simulated motions are obtained, should be modified to be more representative of the frequency content of recorded motions, especially in the low frequency range. Such an improved model can be expected to lead to better agreement, over a broad range of periods, in the average response spectra of simulated and recorded motions, for elastic as well as inelastic systems.     

Wavelet analysis on transient behaviour of tidal amplitude fluctuations observed by meteor radar in the lower thermosphere above Bulgaria

On the basis of bispectral analysis applied to the hourly data set of neutral wind measured by meteor radar in the MLT region above Bulgaria it was demonstrated that nonlinear processes are frequently and regularly acting in the mesopause region. They contribute significantly to the short-term tidal variability and are apparently responsible for the observed complicated behavior of the tidal characteristics. A Morlet wavelet transform is proposed as a technique for studying nonstationary signals. By simulated data it was revealed that the Morlet wavelet transform is especially convenient for analyzing signals with: (1) a wide range of dominant frequencies which are localized in different time intervals; (2) amplitude and frequency modulated spectral components, and (3) singular, wave-like events, observed in the neutral wind of the MLT region and connected mainly with large-scale disturbances propagated from below. By applying a Morlet wavelet transform to the hourly values of the amplitudes of diurnal and semidiurnal tides the basic oscillations with periods of planetary waves (1.5/20 days), as well as their development in time, are obtained. A cross-wavelet analysis is used to clarify the relation between the tidal and mean neutral wind variability. The results of bispectral analysis indicate which planetary waves participated in the nonlinear coupling with the atmospheric tides, while the results of cross-wavelet analysis outline their time intervals if these interactions are local.     

Spectral characteristics of magnetogravity waves generated by high energy mass source in the equatorial region of the atmosphere,Part II

In the first part of study the possibility of high frequency (HF) and low frequency (LF) modes of magnetogravity waves (MGW) propagation in the equatorial atmosphere in conditions of ionospheric F2 layer is demonstrated taking into account the finite conductivity. Here the spectral features of the excited HF MGW mode is studied. At the second part of study the dispersion relations for LF MGW in the ionospheric equatorial region is received. The spectral features of the disturbed components of medium displacement for LF MGW generated by a mass source are investigated. The results are applied to calculate the values of the spectral components of magnetogravity density, pressure and magnetic field disturbances caused by horizontal meteor sweep. The estimates may be useful in the analysis of ionospheric disturbances excited by high-energy geophysical sources.     

Identification and analysis of wind speed patterns extracted from multi-sensors measurements

The understanding of the vertical as well as the horizontal behaviours of wind speed is of great importance in many applications such as aviation, meteorology and wind energy conversion. In this work, we propose to apply the principal component analysis (PCA) in order to extract probable components of wind speed. The idea behind the use of PCA is to introduce mixed sources signals to PCA algorithm as input in order to obtain a separated patterns as output. Hence, values of wind speed measured at three levels above the ground will be used as three separate sensors in order to extract the horizontal and the vertical components of wind speed. Once the principal components of wind speed separated, a process of recognition and identification is undertaken via the inspection of the statistical as well as the cyclical behaviors of the obtained components. For the examination of the statistical properties of wind speed, we propose to carry a comparison of the probability density of the extracted components with the Weibull distribution (commonly used to fit wind speed distributions). However, the spectral behavior of the obtained patterns is examined using time–frequency analysis rather than the traditional Fourier analysis. The time–frequency analysis has been chosen as it serves the purpose of following the diurnal and seasonal time variation of the wind speed spectrum. As a result, it has been found that the horizontal wind speed component fits the Weibull distribution and it is characterized by synoptic and intra-seasonal oscillations. On the other hand, the wind speed vertical component is better fitted by the extreme value distribution. It has been also found that the diurnal oscillations are the most significant oscillations in the vertical components especially in the summertime period.     

Site response in the Qionghai Basin in the Wenchuan earthquake

Amplification effects of soil site response can significantly impact ground motions, and must be considered in the seismic fortification of buildings/structures to prevent or mitigate this potential seismic hazard. Utilizing acceleration time histories from the main shock of the Wenchuan earthquake recorded at four stations (i.e., one on bedrock and three on soil) in the Qionghai Basin, the site responses from three soil sites are studied by using the traditional spectral ratio method. The bedrock site is selected as a reference site. This study found that peak ground accelerations (PGAs) on the soil sites are much larger than on bedrock, with EW, NS and UD components of 3.96–6.58, 6.27–10.98, and 3.17–6.66 times those of the bedrock site, respectively. The amplification effects of the soil sites on ground motions in the frequency range of 0.1 Hz to 10 Hz are significant, depending on the thickness of the soil layer and the frequency content of the site. A significant amplification occurs with high frequency components of ground motion at shallow soil sites, and low and high frequency components of ground motion at intermediate soil sites.     

Seismic analysis of underground tunnels by the 2.5D finite/infinite element approach

A procedure for the seismic analysis of underground tunnels using recorded free-field earthquakes based on the 2.5D finite/infinite element approach is presented. The near and far fields of the half space are modeled by finite and infinite elements, respectively. Using the 1D wave theory, the nodal force and displacement on the near-field boundary are computed for each spectral frequency of the earthquake. Then, equivalent seismic forces are computed for the near-field boundary for the earthquake spectrum. By assuming the soil-tunnel system to be uniform along the tunnel axis, the 2.5D approach can account for the wave transmission along the tunnel axis, which reduces to the 2D case for infinite transmission velocity. The horizontal and vertical components of the 1999 Chi-Chi Earthquake (TCU068) are adopted as the free-field motions in the numerical analysis. The maximal stresses and distribution patterns of the tunnel section under the P- and SV-waves are thoroughly studied by the 2.5D and 2D approaches, which should prove useful to the design of underground tunnels.     

连续小波变换系数的主成分分析法及其重力固体潮应用

针对连续小波变换中各尺度系数间存在较大相关性的特点,建立协方差矩阵,用主成分分析法去除各尺度系数间的相关性,提取原信号中的主要频率成分。利用给定样本信号对该方法的有效性进行检验,同时,作为应用,利用该方法提取出重力固体潮信号中的半日波分量和日波分量,为信号特征频率提取提供了一种可行和有效的分析方法。     

时程信号的Hilbert-Huang变换与小波分析/

在简要介绍时程信号的小波分析和Hilbert-Huang变换(HHT)理论的基础上，通过地震波和其它时程信号实例，对比分析了小波变换和HHT变换结果. 比较显示:HHT变换和小波变换均能用于对非平稳的信号进行分析，并能捕捉到信号变化的主要特征；与受所选母波影响较大的小波分析不同，HHT变换得到的固有模态函数是直接从原始时程数据中分离出来的，它更能反映原始数据的固有特性；小波分析得到的谱的能量在频率范围内分布较广，而HHT变换的Hilbert谱的大部分能量都集中在一定的时间和频率范围内，能清晰地刻画信号能量随时间、频率的分布. 因此，Hilbert-Huang变换不仅是对非平稳信号进行分析的有效方法，而且也是检测时程信号局部特征的有用工具.      

Shear Wave Splitting Analysis to Estimate Fracture Orientation and Frequency Dependent Anisotropy

Shear wave splitting is a well-known method for indication of orientation, radius, and length of fractures in subsurface layers. In this paper, a three component near offset VSP data acquired from a fractured sandstone reservoir in southern part of Iran was used to analyse shear wave splitting and frequency-dependent anisotropy assessment. Polarization angle obtained by performing rotation on radial and transverse components of VSP data was used to determine the direction of polarization of fast shear wave which corresponds to direction of fractures. It was shown that correct implementation of shear wave splitting analysis can be used for determination of fracture direction. During frequencydependent anisotropy analysis, it was found that the time delays in shearwaves decrease as the frequency increases. It was clearly demonstrated throughout this study that anisotropy may have an inverse relationship with frequency. The analysis presented in this paper complements the studied conducted by other researchers in this field of research.     

应用小波变换提取张北地震的震磁效应

以1998年1月10日的张北地震为例,应用小波变换及傅里叶分析,分析震中附近的6个地磁台站在地震前后共6个月的地磁场三分量分钟值采样数据.发现各个地磁台站记录的地磁场Z分量在张北地震前三个月内出现不同程度的异常现象,这种现象的出现频率随着离发震时间越近而越来越高,出现的地点由南向北迁移.地磁场H和D分量没有此类异常出现.用这6个地磁台出现异常时同一频率所对应的振幅值绘出等值曲线,得到的图形是稳定的四象限分布,震中基本上位于0～3nT的等值线上.结果表明小波分析是提取地震磁效应的有效工具,同时也证实了地震磁效应的客观存在.     

Wavelet filter analysis of splitting and coupling of seismic normal modes below 1.5 mHz with superconducting gravimeter records after the December 26, 2004 Sumatra earthquake

New generation superconducting gravimeters (SGs), which have been demonstrated to be better than the best seismometers STS-1 at frequencies below 1 mHz, can be accepted as the quietest vertical seismometers for observation of long-period earth free oscillations. Wavelet filtering with narrow band-pass frequency response as shown in this paper is very helpful in removing at- mospheric pressure effects from on gravity records in long-period seismic mode frequency bands. The processing of high quality SG records after the great Sumatra earthquake (Dec. 26, 2004) with wavelet filtering leads to clear observations of all coupled toroidal modes below 1.5 mHz except these for 0T5, 0T7 and 1T1; moreover 1T2 and 1T3 are, for the first time, unambiguously revealed in the vertical components of the free oscillations. The three well-resolved splitting singlets of overtones 2S1 are observed from a single SG record for the first time.     

小波包单子带重构提取GPS序列特征信息

GPS观测环境愈来愈复杂,动态观测值包含的影响因素较多,函数关系复杂,影响特征信息的提取和参数模型的解释能力.小波包具有良好的时频分析能力,利用小波包理论对GPS数据序列进行分解与重构过程中有三个基本运算：与小波滤波器卷积、隔点采样、隔点插零,该三项运算产生频率交错和频率折叠等频率混淆现象.为消除频率混淆现象,分解与重构时,每作一次信号与小波卷积后,将其结果作一次快速傅立叶变换,频谱中多余的频率成分的谱值置零,再对置零后的频谱进行傅立叶逆变换,然后继续进行小波包的分解与重构,从而实现单子带重构提取GPS数据序列特征项.通过实例验证了小波包单子带重构提取GPS特征信息的有效性.     

复数子波匹配追踪算法识别薄层砂体

针对东营凹陷滩坝砂岩薄砂体的预测识别问题,研究了复小波匹配追踪方法.从地震信号morlet小波分解及重构分析入手,讨论了重构地震信号时关于频率、振幅、相位及时移参数的选取方法及其具体的选取变量.并就关键的振幅参数,研究了改善的计算方法.常规的匹配追踪算法扫描地震信号的所有时间点和频率.在本文中,为了改善时频分辨率,采用局部扫描主频和子波的时间延迟.在具体实现过程中,采用计算地震信号的平均瞬时频率和地震信号的瞬时包络处的时间,做为采样频率和时间采样点,得到与地震信号最佳相关的子波.对研究的方法从理论模型和实际资料进行了论证分析,证明理论上正确合理,实际应用效果显著.     

Characteristic scales, temporal variability modes and simulation of monthly Elbe River flow time series at ungauged locations

Spatial and temporal patterns of the long-range extreme monthly Elbe River flows across Germany are investigated, using various statistical methods, among others, principal component and wavelet analysis. Characteristic time scales are derived for various time series statistics. The wavelet analysis of the raw river discharge data as well as of the major principal component reveal the main oscillatory components and their temporal behavior, namely low frequency oscillations at interannual (6.9 yr) and interdecadal (13.9 yr) scales. The EOFs at ungauged stations are estimated from the principal components of the observed time series sampled over a limited time span whose length equals the major temporal variability scale (≈7 yr). The EOFs (empirical orthogonal functions) obtained in this way are subsequently used to simulate long-range flows at these locations. A comparison of this method with linear interpolation and ordinary kriging of the EOF shows the superiority of the former in representing the distributional properties of the observed time series. The simulated time series preserve also short and long-memory.     

地震低频伴影的数值模拟与应用

Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-frequency shadow is still unclear. To stud） the mechanism, we performed seismic numerical simulation of geological models with a hydrocarbon-bearing zone using the 2-D diffusive-viscous wave equation which car effectively model the characteristics of velocity dispersion and transform the seismic dat~ centered in a target layer slice within a time window to the time-frequency domain by usinl time-frequency signal analysis and sort the frequency gathers to common frequency cubes. Then, we observe the characteristics of the seismic low-frequency shadow in the common frequency cubes. The numerical simulations reveal that the main mechanism of seismic lowfrequency shadows is attributed to high attenuation of the medium to high seismic frequency components caused by absorption in the hydrocarbon-filled reservoir. Results from a practical example of seismic low-frequency shadows show that it is possible to identify the reservoir by the low-frequency shadow with high S/N seismic data.