随机地震动模拟的时间序列法及其工程应用

系统研究了AR和ARMA时间序列方法生成给定加速度功率谱人工地震波的新方法，根据相应理论编制了程序，并生成了人工地震加速时程曲线，然后用生成的人工地震波进行结构动力时程分析，与实际地震波分析的结果进行了对比，计算结果表明，本文用于合成人工地震波的方法是可行的，生成的人工波可以用于实际结构的地震反应分析。     

大型双槽渡槽地震反应分析

为了得到比较精确的大型双槽渡槽结构的地震反应，根据渡槽结构复杂、影响动力因素较多等特点，采用8结点空间块体单元对渡槽槽身和槽墩进行有限元离散，分别用耦合自由度和多个弹簧单元分别模拟渡槽槽身横向拉杆、加劲肋和盆式橡胶支座，建立了大型双槽渡槽结构地层反应分析的有限元模型，计算了南水北调水利工程中的双洎河大型双槽渡槽的振动特性，并分别用反应谱法和时程分析法对该渡槽进行了地震反应分析，计算结果可为该渡槽的抗震设计提供参考。     

高维PP时间序列分析在地震预报中的应用

将投影寻踪(PP)与高维时间序列分析结合起来，建立了地震PP综合预测模型。并选取祁连山地区作研究区，做了未来三个月内最大震级的短期预测，经展望式检验，合格率(预测误差≤0．5或≤10％)≥80％。中根据实际需要提出建立一定震级门限(M≥4．0)的预测具有更好的效果。本预测模型还可增加自变量、维数和改变时间尺度，并运用到其他领域中。     

Weibull分布的信息熵在地震预报中的应用研究

根据统计分布与信息熵理论，定义了震级信息熵Hm和地震间隔时间信息熵Ht，并推导了它们的计算公式。通过时空扫描和计算，发现强震前1～3年Hm和Ht出现低值异常，与地震有较好的对应关系，可以作为一组中期或中短期预测指标。     

面波群速度频散计算中滤波的数值实验研究

面波群速度计算中容易出现长周期部分时间域分辨率低的问题，因此，如何提高时间域的分辨率很重要。本文将引入等误差滤波的方法就改进群速度计算中滤波效果进行数值试验。结果表明，采用等误差滤波能明显提高群速度长周期部分的分辨率，对准确地测定长周期部分的群速度有重要意义。     

Discussion on the feature of strong earthquake: Orderly distribution in time, space and intensity before the Western Kunlun Mountain Pass M=8.1 earthquake

Introduction The Western Kunlun Mountain Pass M=8.1 earthquake occurred on November 14, 2001 is the other M=8 earthquake occurred 50 years after Dangxiong, Tibet M=8.0 earthquake in Chinese mainland. The earthquake has caused the attention of the seismologists in the following aspects: 1) The fracture length is more than 400 km, which is far away from the estimated length by the statistic empirical function between the magnitude and the fracture length (WANG, et al, 2002); 2) The aftersh…     

Tensor time domain electromagnetic resistivity measurements at Ngatamariki geothermal field, New Zealand

Experimental measurements in the Ngatamariki geothermal field, North Island, New Zealand were made to test the applicability of the time domain electromagnetic method for detailed investigation of the resistivity structure within a geothermal field. Low-frequency square wave signals were transmitted through three grounded bipole current sources sited about 8 km from the measurement lines. Despite high levels of electrical noise, transient electric field vectors could be determined reliably for times between 0.02 and 3.3 s after each step in the source current. Instantaneous apparent resistivity tensors were then calculated. Apparent resistivity pseudosections along the two measurement lines show smooth variations of resistivity from site to site. Over most of the field the images consistently show a three-layer resistivity structure with a conductive middle layer (3–10 Ωm) representing the conductive upper part of the thermal reservoir. A deep-seated region of low resistivity in the northwest of the field may indicate a conductive structure at about 1 km associated with a deeper diorite intrusion. Measurements sited closer than about 100 m to drillholes appear to have been disturbed by metallic casing in the holes. A change in resistivity structure in the east of the field may indicate a major geological or hydrothermal boundary.     

Scattering of elastic waves by a general anisotropic basin. Part 2: a 3D model

Scattering of plane harmonic waves by a three‐dimensional basin of arbitrary shape embedded within elastic half‐space is investigated by using an indirect boundary integral equation approach. The materials of the basin and the half‐space are assumed to be the most general anisotropic, homogeneous, linearly elastic solids without any material symmetry (i.e. triclinic). The unknown scattered waves are expressed in terms of three‐dimensional triclinic time harmonic full‐space Green's functions. The results have been tested by comparing the surface response of semi spherical isotropic and transversely isotropic basins for which the numerical solutions are available. Surface displacements are presented for a semicircular basin subjected to a vertical incident plane harmonic pseudo‐P‐, SV‐, or SH‐wave. These results are compared with the motion obtained for the corresponding equivalent isotropic models. The results show that presence of the basin may cause significant amplification of ground motion when compared to the free‐field displacements. The peak amplitude of the predominant component of surface motion is smaller for the anisotropic basin than for the corresponding isotropic one. Anisotropic response may be asymmetric even for symmetric geometry and incidence. Anisotropic surface displacement generally includes all three components of motion which may not be the case for the isotropic results. Furthermore, anisotropic response strongly depends upon the nature of the incident wave, degree of material anisotropy and the azimuthal orientation of the observation station. These results clearly demonstrate the importance of anisotropy in amplification of surface ground motion. Copyright © 2003 John Wiley & Sons, Ltd.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Simplified analysis and supplemental damping system design

Investigated is the accuracy in estimating the response of asymmetric one‐storey systems with non‐linear viscoelastic (VE) dampers by analysing the corresponding linear viscous system wherein all non‐linear VE dampers are replaced by their energy‐equivalent linear viscous dampers. The response of the corresponding linear viscous system is determined by response history analysis (RHA) and by response spectrum analysis (RSA) extended for non‐classically damped systems. The flexible and stiff edge deformations and plan rotation of the corresponding linear viscous system determined by the extended RSA procedure is shown to be sufficiently accurate for design applications with errors generally between 10 and 20%. Although similar accuracy is also shown for the ‘pseudo‐velocity’ of non‐linear VE dampers, the peak force of the non‐linear VE damper cannot be estimated directly from the peak damper force of the corresponding linear viscous system. A simple correction for damper force is proposed and shown to be accurate (with errors not exceeding 15%). For practical applications, an iterative linear analysis procedure is developed for determining the amplitude‐ and frequency‐dependent supplemental damping properties of the corresponding linear viscous system and for estimating the response of asymmetric one‐storey systems with non‐linear VE dampers from the earthquake design (or response) spectrum. Finally, a procedure is developed for designing non‐linear supplemental damping systems that satisfy given design criteria for a given design spectrum. Copyright © 2003 John Wiley & Sons, Ltd.     

Optimum multiple tuned mass dampers for structures under the ground acceleration based on the uniform distribution of system parameters

The five MTMD models, with natural frequencies being uniformly distributed around their mean frequency, have been recently presented by the first author. They are shown to have the near‐zero optimum average damping ratio (more precisely, for a given mass ratio there is an upper limit on the total number, beyond which the near‐zero optimum average damping ratio occurs). In this paper, the eight new MTMD models (i.e. the UM‐MTMD1～UM‐MTMD3, US‐MTMD1～US‐MTMD3, UD‐MTMD1 and UD‐MTMD2), with the system parameters (mass, stiffness and damping coefficient) being, respectively, uniformly distributed around their average values, have been, for the first time here, proposed to seek for the MTMD models without the near‐zero optimum average damping ratio. The structure is represented by the mode‐generalized system corresponding to the specific vibration mode that needs to be controlled. Through minimization of the minimum values of the maximum dynamic magnification factors (DMF) of the structure with the eight MTMD models (i.e. through the implementation of Min.Min.Max.DMF), the optimum parameters and values of Min.Min.Max.DMF for these eight MTMD models are investigated to evaluate and compare their control performance. The optimum parameters include the optimum mass spacing, stiffness spacing, damping coefficient spacing, frequency spacing, average damping ratio and tuning frequency ratio. The six MTMD models without the near‐zero optimum average damping ratio (i.e. the UM‐MTMD1～UM‐MTMD3, US‐MTMD1, US‐MTMD2 and UD‐MTMD2) are found through extensive numerical analyses. Likewise, the optimum UM‐MTMD3 offers the higher effectiveness and robustness and requires the smaller damping with respect to the rest of the MTMD models in reducing the responses of structures subjected to earthquakes. Additionally, it is interesting to note, by comparing the optimum UM‐MTMD3 with the optimum MTMD‐1 recently investigated by the first author, that the effectiveness and robustness for the optimum UM‐MTMD3 is almost identical to that for the optimum MTMD‐1 (without inclusion of the optimum MTMD‐1 with the near‐zero optimum average damping ratio). Recognizing these performance benefits, it is preferable to employ the optimum UM‐MTMD3 or the optimum MTMD‐1 without the near‐zero optimum average damping ratio, when installing the MTMD for the suppression of undesirable oscillations of structures under earthquakes. Copyright © 2003 John Wiley & Sons, Ltd.