考虑超随机特性的人工地震动合成

文中根据自相关函数随机噪声部分的2个重要性质提出了自相关函数随机噪声部分的迭代生成算法,并给出了以确定性自相关函数为目标函数,考虑随机噪声部分影响的人工合成地震动方法.本文的研究结果表明地震动加速度的平方和（SS）比峰值加速度（PGA）更好地量度了地震地面运动强度.文中用蒙特卡罗方法分析了自相关函数噪声项对人工合成地震动幅值和反应谱的影响：由于自相关函数的随机噪声项对地震动总功率或平方和的贡献为0,因此随机噪声项对合成地震动幅值的影响不大,但由于自相关函数随机噪声项明显改变了地震动功率谱在频域上的分布特征,使原本光滑的平均反应谱出现了明显的波动,从而在一定程度上提高了平均反应谱的峰值,使自相关函数噪声部分对平均反应谱也有一些影响.     

Stochastic seismic response of structures with added viscoelastic dampers modeled by fractional derivative

Viscoelastic dampers, as supplementary energy dissipation devices, have been used in building structures under seismic excitation or wind loads. Different analytical models have been proposed to describe their dynamic force deformation characteristics. Among these analytical models, the fractional derivative models have attracted more attention as they can capture the frequency dependence of the material stiffness and damping properties observed from tests very well. In this paper, a Fourier-transform-based technique is presented to obtain the fractional unit impulse function and the response of structures with added viscoelastic dampers whose force-deformation relationship is described by a fractional derivative model. Then, a Duhamel integral-type expression is suggested for the response analysis of a fractional damped dynamic system subjected to deterministic or random excitation. Through numerical verification, it is shown that viscoelastic dampers are effective in reducing structural responses over a wide frequency range, and the proposed schemes can be used to accurately predict the stochastic seismic response of structures with added viscoelastic dampers described by a Kelvin model with fractional derivative.     

SEISMIC PHASE UNWRAPPING: METHODS,RESULTS, PROBLEMS1

Six known methods of seismic phase unwrapping (or phase restoration) are compared. All the methods tested unwrap the phase satisfactorily if the initial function is a simple theoretical wavelet. None of the methods restore the phase of a synthetic trace exactly. An initial validity test of the phase-unwrapping method is that the sum of the restored wavelet phase spectrum and the restored pulse-trace phase spectrum (assuming the convolutional model for the seismic trace) must be equal to the restored phase spectrum of the synthetic trace. Results show that none of the tested methods satisfy this test. Quantitative estimation of the phase-unwrapping accuracy by correlation analysis of the phase deconvolution results separated these methods, according to their efficiency, into three groups. The first group consists of methods using a priori wavelet information. These methods make the wavelet phase estimation more effective than the minimum-phase approach, if the wavelet is non-minimum-phase. The second group consists of methods using the phase increment Δø(Δω) between two adjacent frequencies. These methods help to decrease the time shift of the initial synthetic trace relative to the model of the medium. At the same time they degrade the trace correlation with the medium model. The third group consists of methods using an integration of the phase derivative. These methods do not lead to any improvement of the initial seismic trace. The main problem in the phase unwrapping of a seismic trace is the random character of the pulse trace. For this reason methods based on an analysis of the value of Δø(Δω) only, or using an adaptive approach (i.e. as Δω decreases) are not effective. In addition, methods based on integration of the phase derivative are unreliable, due to errors in numerical integration and differentiation.     

强地震动幅值谱随机涨落特征的研究

本文从工程应用的观点，将强地面运动幅值谱?Ａ（ｆ）│表示为强度受到调制的随机过程即│Ａ（ｆ）│＝│Ａ（ｆ）│〔１＋ｎ（ｆ）〕其中，│Ａ（ｆ）│为随频率变化的确定性函数亦是│Ａ（ｆ）│的均值，ｎ（ｆ）为均值为零的随机过程。在随机过程ｎ（ｆ）为各态历经的假定下，研究了ｎ（ｆ）的功率谱密度函数的特征及其衰减形式。根据地震学的有关理论，提出了用四个参数控制幅值谱的均值随频率变化的函数形式并对拟合这些参数的     

THE PHASE PROPERTY OF A FINITE REALIZATION OF RANDOM NOISE AND ITS INFLUENCE ON DECONVOLUTION *

A finite realization of a discrete random noise process may be considered as a one-sided energy signal. Its phase property can then be described by means of the center position. The samples of such a realization are the components of a random signal vector and the center position is therefore a random variable. A statistical analysis shows that the expected value of the center position equals half the time duration of the realization. This implies that the Z-transform of the realization may be expected to have an equal number of poles and zeros inside and outside the unit circle. The standard deviation from the expected value of the center position is shown to depend on the time duration of the realization and on the autocorrelation of the process. It follows that, for processes that can be described by the convolution of a white series and a disturbance wavelet, the center position is independent of the phase property of the wavelet. A conclusion based on these results is that the homomorphic technique of wavelet estimation through cepstrum stacking must give questionable outcomes. Another conclusion is that the super-position of a realization of random noise on a minimum phase wavelet will in general give a mixed phase resulting signal. It is pointed out that schemes for the derivation of deconvolution filters do not take account of this phenomenon.     

陆地地震勘探随机噪声统计特性

在地震勘探随机噪声压制领域,噪声通常被假设为平稳、高斯随机过程的信息.然而,在某些情况下,这样的假设并不准确.本文应用现代统计检验方法对地震勘探随机噪声的平稳性、高斯性和线性进行了研究.结果表明地震勘探随机噪声并不是传统意义上认为的平稳随机过程,其平稳性受到噪声时长和采集环境复杂程度的影响.发现噪声时间越长,采集环境越复杂,随机噪声的平稳性越差,但是对于短时长随机噪声而言,其可以近似认为是平稳的.同时,采集环境的复杂程度也影响着随机噪声的高斯性和线性特性,环境条件越复杂,随机噪声高斯性越好,线性特性越差,但总的来说随机噪声可以归为线性非高斯随机过程.

     

A response spectrum method for random vibration analysis of mdf systems

A response spectrum method for stationary random vibration analysis of linear, multi-degree-of-freedom systems is developed. The method is based on the assumption that the input excitation is a wide-band, stationary Gaussian process and the response is stationary. However, it can also be used as a good approximation for the response to a transient stationary Gaussian input with a duration several times longer than the fundamental period of the system. Various response quantities, including the mean-squares of the response and its time derivative, the response mean frequency, and the cumulative distribution and the mean and variance of the peak response are obtained in terms of the ordinates of the mean response spectrum of the input excitation and the modal properties of the system. The formulation includes the cross-correlation between modal responses, which is shown to be significant for modes with closely spaced natural frequencies. The proposed procedure is demonstrated for an example structure that is subjected to an ensemble of earthquake-induced base excitations. Computed results based on the response spectrum method are in close agreement with simulation results obtained from time-history dynamic analysis. The significance of closely spaced modes and the error associated with a conventional method that neglects the modal correlations are also demonstrated.     

二维MT快速曲线对比反演方法的可行性研究

最近开发了一种针对二维大地电磁野外数据进行处理解释的新反演方法。该方法以加入阻抗相位信息的一维大地电磁连续介质曲线对比法为基础 ,把一维反演得到的电阻率和相位的数据集作为二维反演的初始模型 ,使用二维有限单元法做正演模拟。在程序的后继迭代中 ,深度方向上用一维反演修改模型的电阻率和深度值 ,沿测线方向由二维有限元作修改 ,反演结束可得到一个接近真实电性分布的电阻率数据集 ,并绘制成电阻率断面图。对模型的反演实验结果显示 ,该反演方法能够较真实地反映地下电性分布 ,而且避免了偏导数矩阵的计算 ,其原理简单 ,计算速度快 ,表明该反演方案是可行的     

EFFECTS OF DATA ERRORS AND ACCURACY IN SECOND DERIVATIVE COMPUTATION BY RATIONAL APPROXIMATION*

The effects of systematic (constant) and random errors in the observed data have been investigated analytically for rational approximation method of computing second derivative involving a summation of the products of the averages of the gravity field with the corresponding weight coefficients, both in numerator as well as in denominator. A theoretical gravity anomaly over three spheres has been analyzed to demonstrate the high accuracy in the approximation. Since the sums of the weight coefficients in numerator and denominator are zero and one respectively, the regional gravity anomaly, even though approximated by a constant value over the entire area under computation, can produce substantially large error in the calculated derivative value. This is happening because of the contribution of the regional field in the denominator. Thus, inspite of the high accuracy in rational approximation, the method has limited application to field cases where a combined gravity field consisting of regional and residual anomalies is usually used. Master curves are presented for the constant and random errors by which a rough estimate of the percentage of error in second derivative computation can be made provided one has some idea of the magnitudes of the regional field and random error.     

Simulation of artificial random field considering the temporal-spatial variation

Introduction For the seismic design of special structures such as nuclear power station, marine platform, long-span bridge and dam, generally the time-history response analysis of the structure under seismic excitation is imperative, which was coded in most seismic design codes. The earthquake records suitable for the seismic situation and site condition are necessary to be used as the seismic input in the dynamic analysis of structures. As a result of the limited observational condition of st…     

一种神经网络改进小波的地震数据随机噪声去除方法

地震资料的有效信号反射弱,且易受多次波的影响,不可避免地存在随机噪声干扰。提出一种基于神经网络改进小波的地震数据随机噪声去除方法,采用神经网络模型,识别出随机噪声信号,对该信号进行小波包分解,获取多类别随机噪声信号,采用级联BP神经网络模型提取出多类别随机噪声信号,实现地震数据的随机信号压制。实验结果显示,这种改进小波方法对地震数据随机噪声信号的去噪效果较好,在复杂沉积地质结构被探测介质的地震数据随机噪声压制方面具有较强的适用性。     

地震动随机函数模型研究(Ⅰ)——模型建立

基于随机过程的随机函数解释,发展了一种随机地震动模型.模型包括两个部分,分别是基于理论物理关系的Fourier幅值谱模型和基于经验关系的Fourier相位谱模型.基于拟加速度谱概念建立的幅值谱模型能够反映强震记录的幅值特性,而Fourier相位谱模型则利用累积相位谱的单值连续特性,以线性多项式与三角级数组合的形式给出.模拟记录与实测记录对比结果显示两者不仅在波形结构上相似性良好,而且反应谱层次亦吻合良好.     

基于随机Fourier谱的地震动合成研究

本文基于与物理随机模型相对应的随机Fourier谱,通过对地震动非平稳性及其所受影响机制的分析,建立了非平稳地震动合成的新方法。通过引入基本相位差谱,并将初始相角视为随机变量,以相位差谱的分布特性、随机地震动的统计特征以及对结构随机反应的影响为原则,利用计算机程序生成相位差谱,提出了基于随机Fourier谱的合成地震动方法。同时,利用快速FFT技术提高合成精度。根据本文提出的合成方法获得的地震动具有非平稳特性,将为后续研究工作提供合理的地震动加速度时程。     

随机弹性介质中地震波散射衰减分析(英文)

地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。     

考虑时间-空间变化的人工随机场模拟

采用随频率变化的视波速代替随意给定的视波速,并在随机场的模拟中引入相位差谱来考虑地震动的频率含量非平稳性．用谱表示法按不同的设计烈度生成了时间-空间变化的非平稳人工随机场,可用于大跨度空间结构多点激励的地震动输入．     

Generation of multivariate autoregressive sequences with emphasis on initial values

Certain aspects of data generation are studied through multivariate autoregressive (AR) models. The main emphasis is on the preservation of certain desired moments and the effect of initial values on these moments. The problem of preservation of moments is approached in a nontraditional way by starting with the initial values. For this purpose, general AR processes with a random start and with time-varying parameters are introduced to lay a foundation for the analysis of all types of AR processes, including the periodic cases. It is shown that an AR process with a random start and with parameters obtained from the moment equations is capable of generating jointly multivariate normal vectors with any specified means and covariance matrices, and with any specified autocovariance matrices up to a given lag. With a random start, there is no transition period involved for achieving these moments. A simple solution is proposed for matrix equations of the form BB^T = A which appear in the moment equations. The aggregation properties of general AR process are also studied.

A more detailed analysis is given for the two-period first-order periodic autoregressive model, PAR₂(1). For the PAR₂(1) process with a random start and with parameters obtained from the moment equations, it is shown that the autocovariance function depends only on the period and the lag, and therefore the process is periodic (covariance) stationary. The PAR₂(1) process with a fixed start is also studied. It is shown that the moments of this process depend on the absolute time, in addition to the period and the lag, and therefore the process is not periodic stationary. This dependence diminishes with time, and periodic stationarity is realized if the AR parameters satisfy certain conditions. In that case, the PAR₂(1) process with a fixed start converges to that with a random start, but only after a certain transition period. This proves the superiority of a random start over a fixed start.     

Determination of the phase projection of a dynamic system from strain field observations

The paper analyzes the possibility of extracting a geophysical field component from a time series considered as a phase projection of a dynamic system. Based on model calculations, the method of averaging of phase trajectories is shown to be applicable for the recovery of time variations of a dynamic system and more reliable determination of its quantitative characteristics if the variations are superimposed by a high-dimension process (a purely random process) and a low-dimension process (linear interpolation). The analysis of the extraction of dynamic system variations from strain records obtained at various points of seismically active regions of China provided the following estimates of dynamic parameters: the phase space dimension is m = 5–7, and the correlation dimension varies, depending on the observation point, within the range d ₂ = 2.6–3.5.     

Spatial dispersion of surface waves at the free surface of a general anisotropic elastic medium

A new technique relates the wave velocity of the surface waves in anisotropic elastic medium to its elastic constants. Anisotropic propagation of surface waves is studied in a half-space occupied by a general anisotropic elastic solid. The phase velocity expressions of quasi-waves, in three-dimensional space, are used to derive the secular equation of surface waves. The complex secular equation is resolved, analytically, into real and imaginary parts and is then solved, numerically, for phase velocity along a given phase direction on the surface. The complete procedure is thus analogous to the one used for conventional Rayleigh waves in isotropic medium. A non-linear equation relates the ray direction of the surface waves to its phase direction on the (plane) surface of the medium. The analytical differentiation of secular equation yields the directional derivative of phase velocity. This derivative is used to calculate the wave velocity of surface waves. Spatial variations of phase velocity, wave velocity and ray direction over the free plane surface are plotted for the numerical models of crustal rocks with orthorhombic, monoclinic and triclinic anisotropies.     

Theoretical research on intensity envelope function of strong ground motion

（金星，瘳振鹏）Ｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｎｉｎｔｅｎｓｉｔｙｅｎｖｅｌｏｐｅｆｕｎｃｔｉｏｎｏｆｓｔｒｏｎｇｇｒｏｕｎｄｍｏｔｉｏｎ￥ＸｉｎｇＪＩＮａｎｄＺｈｅｎ－ＰｅｎｇＬＩＡＯ（ＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈ...