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

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

基于脉冲频响函数的MATLAB动力方程求解方法

提出一种基于脉冲频响函数并利用MATLAB求解的方法，由于脉冲频响函数利用卷积公式进行动力微分方法计算，因此，从理论上讲，它适合于任意外力下的动力方程求解，同时，MATAL又提供了卷积函数conv，所以，在MATLAB下求解动力方程将变得十分容易。     

最不利地震动在网壳结构抗震设计中的应用

利用文献提供的最不利地震动，对网壳结构等大跨空间结构进行了分析，分别考察了仅水平输入和三向地震波输入的两种情况，对最不利地震动在网壳结构抗震设计中的适用性进行了探讨，并在一实际工程设计中进行了应用。大跨度空间结构在水平地震动输入时，宜按照的最不利地震动进行设计；在三向地震动同时输入时，应同时分别按所推荐最不利设计地震动和按常规选用的地震动，对大跨度空间结构进行抗震分析。     

Exploration and Study on Deep Crust Structural Characteristics in the Jiashi-Artux Seismic Region

The data from two deep seismic sounding profiles was processed and studied comprehensively. The results show that crnst-mantle structures in the investigated region obviously display layered characteristics and velocity structures and tectonic features have larger distinction in different geological structure blocks. The boundary interface C between the upper and lower crust and Moho fluctuate greatly. The shallowest depths of C (30.0km) and Moho (45.5km) under Jiashi deepen sharply from Jiashi to the western Kunlun mountain areas, where the depths of C and Moho are 44.0km and 70.0km, respectively. The higher velocity structures in the Tarim massif determine its relatively “stable“ characteristics in crust tectonics. The phenomenon in the Jiashi region, where the distribution of earthquake foci mostly range from 20kin to 40kin in depth, may infer that the local uplift of C and Moho interface, anomalonsly lower velocity bodies and deep large faults control earthquake occurrence and seismogenic processes in the Jiashi strong earthquake swarm.     

The inversion of anelastic coefficient,source parameters and site respond using genetic algorithm

It gradually becomes a common work using large seismic wave data to obtain source parameters, such as seismic moment, break radius, stress drop, with completingof digital seismic network in China (Hough, et al, 1999; Bindi, et al, 2001). These parameters are useful on earthquake prediction and seismic hazard analysis.Although the computation methods of source parameters are simple in principle and the many research works have been done, it is not easy to obtain the parameters accurately. There are two factors affecting the stability of computation results. The first one is the effect of spread path and site respond on signal. According to the research results, there are different geometrical spreading coefficients on different epicenter distance. The better method is to introduce trilinear geometrical spreading model (Atkinson, Mereu, 1992; Atkinson, Boore, 1995; WONG, et al, 2002). In addition, traditional site respond is estimated by comparing with rock station, such as linear inversion method (Andrews, 1982), but the comparative estimation will introduce some errors when selecting different stations. Some recent research results show that site respond is not flat for rock station (Moya, et al, 2000; ZHANG,. et al, 2001; JIN, et al, 2000; Dutta, et al, 2001). The second factor is to obtain low-frequency level and corner frequency fromdisplacement spectrum. Because the source spectrum model is nonlinear function,these values are obtained by eye. The subjectivity is strong. The small change of corner frequency will affect significantly the result of stress drop.     

Theoretical Study on Stochastic Modeling of Combined Gravity-Magnetic-Electric-Seismic Inversion and Its Application

As gravity field,magnetic field,electric field and seismic wave field are all physical fields,their object function,reverse function and compound function are certainly infinite contiuously differentiable func-tions which can be expanded into Taylor (Fourier) series within domain of definition and be further reduced in-to solving stochastic distribution function of series and statistic inference of optimal approximation,This is the basis of combined gravity-magnetic-electric-seismic inversion built on the basis of separation of field and source gravity-magnetic difference-value(D-value)trend surface,taking distribution-independent fault sys-tem as its unit,depths of seismic and electric interfaces of interests as its corresponding bivariate compound re-verse function of gravity-magnetic anomalies and using high order polynomial(high order trigonometric func-tion)approximating to its series distribution,The difference from current dominant inversion techniques is that,first,it does not respectively create gravity-seismic,magnetic-seismic deterministic inversion model from theoretical model,but combines gravity-magnetic-electric-seismic stochastic inversion model from stochastic model;second,after the concept of equivalent geological body being introduced,using feature of independent variable of gravity-magnetic field functions,taking density and susceptibility related to gravity-magnetic func-tion as default parameters of model,the deterministic model is established owing to better solution to the con-tradictioc of difficulty in identifying strata and less test analytical data for density and susceptibility in newly explored area;third,under assumption of independent parent distribution,a real modeling by strata,the prob-lem of difficult plane closure arising in profile modeling is avoided,This technology has richer and more detailed fault and strata information than sparse pattern seismic data in newly explored area,successfully inverses and plots structural map of Indosinian discontinuty in Hefei basin with combined gravity-magnetic-electric-seismic inversion,With development of high precision gravity-magnetic and overall geophysical technology,it is certain for introducing new methods of stochastic modeling and computational intelligence and promoting the develop-ment of combined gravity-magnetic-electric-seismic inversion to open a new substantial and promoting the develop-ment of combined gravity-magnetic-electric-seismic inversion to open a new substantial path.     

团树断层场地地震反应的数值模拟

团树断层场地在多次地震中出现了烈度异常。为了弄清烈度异常的原因，应用显式有限元方法，数值模拟了在脉冲地震波斜入射情形下断层场地的地震反应。结果表明：在基岩山脚附近土层场地地震动放大比基岩山区平均高4倍、比远离基岩山脚土层场地平均高2倍，与历史震害分布有较好的对应。     

前驱地震波的试验研究及波谱分析

研究了含预置裂纹的有机玻璃试样在加载实验过程中，大破裂之前出现的长周期、小振幅低频事件，并与大地震前可能出现的同类现象（即前驱地震波）进行了对比研究，根据波形记录和频谱分析得出：前驱地震波是地壳介质的长周期波动，是一种客观存在的前兆特征，具有短临前兆意义。     

沂沭断裂带中部地区地震波速及地壳厚度的测定

选取1990～2001年发生在沂沭断裂带中部地区的23次ML2．3～4．0级地震为基础资料，对沂沭断裂带中段及附近地区地震波速进行了研究。结果表明，研究区-↑Vp=5．83km／s，-↑Vs=3．54km／s。同时，根据反射波走时方程，得出该地区反射波速度Vp11=6．21km／s，Vs11=3．65km／s和平均地壳厚度为37km。     

检层法钻孔S波速测量中的震相问题

在岩土工程测量中，Ｓ波速度是一个重要参数，正确地识别其震相，对提高Ｓ波速度的测试精度非常重要。郭铁拴先生等发表关于“检层法鲇孔波速测量中的震相分析”一文中提出在有套管的检层Ｓ波测量中，深度大于４０ｍ时，Ｓ波震相淹没在折射波中。