巨厚沉积土夹火山岩场地非线性地震反应特性

地震基岩深度和土体动力本构模型的选取对核岛场地地震效应评价结果的合理性具有重要影响。以拟建某沿海核电厂深度470 m沉积土夹火山岩层场地的3个钻孔剖面为研究对象,采用一维等效线性波传分析（ELA）法、基于Matasovic本构模型和Davidenkov-Chen-Zhao（DCZ）本构模型的一维非线性分析（NLA）法,选取不同剪切波速的5个岩土层作为地震基岩,研究了输入地震动特性、地震基岩深度和土体动力本构模型的选取对巨厚沉积土夹火山岩层场地非线性地震反应特性的影响。结果表明：（1）以浅层硬岩夹层或深部土层作为地震基岩,NLA法计算的5%阻尼比的地表谱加速度SA的短周期部分较之ELA法的计算值大,但两者计算的地表SA谱的长周期部分几乎一致;（2）基于Matasovic模型和DCZ模型的NLA法计算的地表SA谱谱形和峰值加速度随深度的变化趋势基本一致;（3）从NLA法计算的地表峰值加速度和累积绝对速度而言,以剪切波速约2 500 m/s的浅层硬岩夹层作为地震基岩是适宜的。     

苏州第四纪深厚沉积层一维等效线性和非线性地震反应对比分析

合理描述土体动力本构关系对场地设计地震动参数取值的准确性有显著影响。以苏州城区200 m深的钻孔剖面为研究对象,对比分析了一维等效线性波传分析法（ELA法）和基于修正Matasovic本构模型的多自由度集中质量非线性分析法（NLA法）给出的深厚场地地震反应,研究了基岩输入地震动特性和地震基岩面的选取对深厚场地地震反应的影响。结果表明：（1）随基岩地震动强度（PGA）的增大,ELA法给出的地表PGA呈单调递增的特征,而NLA法给出的地表PGA呈先快速增大后缓慢减小或几乎不变的趋势;（2）ELA法和NLA法给出的地表加速度反应谱（Sa谱）在短周期范围内存在明显差异,ELA法对基岩高频地震动具有显著的滤波效应,而NLA法对基岩高频地震动的影响表现为随基岩PGA的增大先放大后减小的特征;（3）随地震基岩面深度的增大,地表Sa谱的谱值呈现出稍许增大的趋势,但对周期T<0.1 s部分,NLA法给出的地表Sa谱的谱值则呈现出稍许减小的现象;（4）中、大震作用下,地表地震动持时不仅与基岩地震动特性有关,还与地震反应分析方法和地震基岩面的选取密切相关,NLA法能更合理地反映基岩地震动强度和上覆土层厚度对地表地震动持时的影响。     

Evidences of landslide earthquake triggering due to self-excitation process

The basin-like setting of stiff bedrock combined with pre-existing landslide masses can contribute to seismic amplifications in a wide frequency range (0–10 Hz) and induce a self-excitation process responsible for earthquake-triggered landsliding. Here, the self-excitation process is proposed to justify the far-field seismic trigger of the Cerda landslide (Sicily, Italy) which was reactivated by the 6th September 2002 Palermo earthquake (M _s = 5.4), about 50 km far from the epicentre. The landslide caused damage to farm houses, roads and aqueducts, close to the village of Cerda, and involved about 40 × 10⁶ m³ of clay shales; the first ground cracks due to the landslide movement formed about 30 min after the main shock. A stress–strain dynamic numerical modelling, performed by FDM code FLAC 5.0, supports the notion that the combination of local geological setting and earthquake frequency content played a fundamental role in the landslide reactivation. Since accelerometric records of the triggering event are not available, dynamic equivalent inputs have been used for the numerical modelling. These inputs can be regarded as representative for the local ground shaking, having a PGA value up to 0.2 m/s², which is the maximum expected in 475 years, according to the Italian seismic hazard maps. A 2D numerical modelling of the seismic wave propagation in the Cerda landslide area was also performed; it pointed out amplification effects due to both the structural setting of the stiff bedrock (at about 1 Hz) and the pre-existing landslide mass (in the range 3–6 Hz). The frequency peaks of the resulting amplification functions (A(f)) fit well the H/V spectral ratios from ambient noise and the H/H spectral ratios to a reference station from earthquake records, obtained by in situ velocimetric measurements. Moreover, the Fourier spectra of earthquake accelerometric records, whose source and magnitude are consistent with the triggering event, show a main peak at about 1 Hz. This frequency value well fits the one amplified by the geological setting of the bedrock in correspondence with the landslide area, which is constituted of marly limestones and characterised by a basin-like geometry.     

Evaluating the structural model from ambient vibration recordings and shallow seismic refraction techniques. A case study of Zagazig City,Egypt

This study uses the seismic refraction and noise measurements to investigate the velocity structure of the subsurface and emphasize the advantage of ambient vibration over the conventional seismic refraction technique. Field measurements were carried out at nine sites in and around Zagazig city. Shallow seismic refraction data were interpreted using the delay time method to obtain the two-dimension ground model at each site. Ambient vibration arrays are used to infer the one-dimensional compressional and shear wave velocity profiles through two main steps. The first step is to derive the dispersion curve from the recorded signals using the frequency-wavenumber method. The second is to invert the dispersion curve to obtain the site velocity profiles. The results of the compressional wave velocities obtained from seismic refraction technique showed that the subsurface consists of a number of layers ranging from two to four layers and give a good agreement with the results of the seismic wave velocities obtained from the ambient vibration arrays. The ambient vibration arrays gave a deeper depth of penetration than the other method, providing more information on the subsurface structure without any disturbance to the environment. This work provides reliable estimates of the seismic velocity structures of both shallow and deep sedimentary layers within the area of interest.     

基于 GIS 的智能辅助矿产预测系统研究

本文是在“中大比例尺矿床统计预测智能辅助系统（MILASP系统）”的研制中完成的。其中对GIS 与Internet的关系作了论述‚推导并应用了“扩展证据加权模型”‚对陕西省勉略宁地区进行了实例应用 研究‚初步实现了矿床统计预测工作的计算机化和自动化‚加强了智能辅助支持功能的开发‚集GIS、 RS、ES3种功能为一体‚代表了专用型GIS软件开发的最新方向。     

Equivalent-Linear Seismic Ground Response Analysis of Some Typical Sites in Mumbai

Any earthquake event is associated with a rupture mechanism at the source, propagation of seismic waves through underlying rock and finally these waves travel through the soil layers to the particular site of interest. The bedrock motion is significantly modified at the ground surface due to the presence of local soil layers above the bedrock beneath the site of interest. The estimation of the amplifications in ground response due to the local soil sites is a complex problem to the designers and the problem is more important for mega cities like Mumbai in India, where huge population may get affected due to devastations of earthquake. In the present study, the effect of local soil sites in modifying ground response is studied by performing one dimensional equivalent-linear ground response analysis for some of the typical Mumbai soil sites. Field borelog data of some typical sites in Mumbai city viz. Mangalwadi site, Walkeswar site, BJ Marg near Pandhari Chawl site are considered in this study. The ground responses are observed for range of input motions and the results are presented in terms of surface acceleration time history, ratio of shear stress to vertical effective stress versus time, acceleration response spectrum, Fourier amplitude ratio versus frequency etc. The typical amplifications of ground accelerations considering four strong ground motions with wide variation of low to high MHA, frequency contents and durations are obtained. Results show that MHA, bracketed duration, frequency content have significant effects on the amplification of seismic accelerations for typical 2001 Bhuj motion. The peak ground acceleration amplification factors are found to be about 2.50 for Mangalwadi site, 2.60 for Walkeswar site and 3.45 for BJ Marg site using 2001 Bhuj input motion. The response spectrum along various soil layers are obtained which will be useful for designers for earthquake resistant design of geotechnical structures in Mumbai for similar sites in the absence of site specific data.     

地面地震CT在浅层勘探中的应用

在地表面沿研究区域的四周布置激发点和检波器 ,进行人工地震观测记录。地震能量由激发点到达接收点 ,利用直达波、反射波和折射波等多震相 ,进行 CT反演 ,得到波速地震 CT图像。该图像给出浅层波速横向分布 ,反映了覆盖层分区及基底岩性特点与构造分布。该方法应用于某工程场地的浅层勘探中 ,所得波速 CT图像与地质露头及开挖相吻合。该方法将地震 CT技术从钻孔和巷道中解放出来。     

Basin Structure and Numerical Simulation for the Mechanisms of Seismic Disasters

In the present study seismic wave propagation in heterogeneous media is numerically simulated by using the pseudospectral method with the staggered grid RFFT differentiation in order to clarify the cause for the complicated distribution characteristics of strong ground motion in regions with basin structure. The results show that the maximum amplitudes of simulated ground acceleration waveforms are closely related to the basin structure. Interference of seismic waves in the basin strongly affects the distribution of maximum seismic waveforms, which may result in peak disasters during earthquakes. Peak disasters might be away from basin boundaries or earthquake faults. Seismic energy transmitted into the basin from the bedrock can hardly penetrate the bottom of the basin and then travel back into the bedrock region. The seismic energy is absorbed by basin media, and transferred into the kinematical energy of seismic waves with great amplitude in the basin. Seismic waves between basins may result in seriou     

波的传播及其与介质间断相互作用的数值模拟研究

本文用有限差分法模拟研究二维空间中波的传播、散射及其与介质间断的相互作用。文章第一部分运用所模拟的纵波波源讨论了半无限空间表面台阶地形或槽形裂隙所造成的波的散射。这些间断可模拟具有波长量级的悬崖或狭谷。第二部分用差分法与扰动法结合而成的混合法处理半无限空间表面有另一种物质构成的表面浅障碍时波的传播,以模拟山脊、堑沟等地形的影响。此混合法也被用于讨论半无限空间内含有水平夹层或孔隙的问题。所得研究结果对波的传播及其与介质间断的相互作用的理论研究和地震波破坏效应的分析和隔离等应用研究有一定的意义。     

Three-dimensional seismic modelling of crustal structure in the TESZ region based on POLONAISE''97 data

This paper reports the results of 3-D tomographic modelling of crustal structure in the Trans European Suture Zone region (TESZ) of Poland, eastern Germany and Lithuania. The data are the product of a large-scale seismic experiment POLONAISE'97, which was carried out in 1997. This experiment was designed to provide some 3-D coverage. The TESZ forms the boundary between the Precambrian crustal terranes of the East European Craton (EEC) and the younger Phanerozoic terranes to the southwest. The 3-D results generally confirm the earth models derived by earlier 2-D analyses, but also add some important details as well as a 3-D perspective on the structure. The velocity model obtained shows substantial horizontal variations of crustal structure across the study area. Seismic modelling shows low (<6.1 km/s) velocities suggesting the presence of sedimentary rocks down to a depth of about 20 km in the Polish basin. The shape of the basin in the vicinity of the profile P4 shows significant asymmetry. Three-dimensional modelling also allowed tracing of horizontal irregularities of the basin shape as well as variations of the Moho depth not only along profiles, but also between them. The slice between P2 and P4 profiles shows about 10-km variations of the Moho over a 100-km interval. The crustal thickness varies from about 30 km in SW, beneath the Palaeozoic platform, to about 42 km beneath East European Craton in NE. High seismic velocities of about 6.6 km/s were found in the depth range 2–10 km, which coincides with K trzyn anorthosite massif. The results of this 3-D seismic modelling of the POLONAISE'97 data will ultimately be supplemented by inversion of seismic data from previous experiments.     

河口盆地非线性地震效应及设计地震动参数

针对基岩明显起伏、土层非均匀分布的典型河口盆地场地,考虑土体非线性特征,采用黏弹性人工边界模拟无限域对地震波动的影响,建立大尺度精细化二维有限元模型,分析了盆地地表地震动幅值、频谱、持时、传递函数特征,探讨了基岩起伏土层的地震动聚集效应及盆地边缘效应。结果表明：（1）盆地近地表土层表现出不同程度的地震动放大效应,且随土层深度增加呈非单调递减特征,基岩突变处地震动聚集效应明显,盆地两侧产生较为显著的边缘效应;场地中、长周期地震动的放大作用显著;（2）多遇地震、偶遇地震和罕遇地震水平时,场地卓越周期依次介于0.35～0.65 s、0.40～0.75 s和0.50～1.05 s之间;给出了盆地地表PGA（地表峰值加速度）、卓越周期均值等值线图及地表加速度反应谱放大因子建议值,地表设计地震动参数amax（地震影响系数）与Tg（特征周期）明显大于现行《建筑抗震设计规范》取值;（3）盆地特殊位置地表地震动持时得到不同幅度增长,且与输入地震动特性相关;（4）该盆地对0.5～2.0 Hz频段基岩地震动的放大效应比较显著,对小于0.2 Hz或大于2.5 Hz的基岩地震动,该盆地地震动放大效应不明显;（5）福州城区及其邻近区域地震动放大效应普遍较大。大尺度二维非线性分析一定程度上能合理反映微地形起伏、土层分布及土体非线性对地震波传播过程的影响。     

Neo-deterministic definition of seismic input for residential seismically isolated buildings

This paper deals with the neo-deterministic definition of the seismic input in the municipality of Nimis (Italy), aimed at the design of residential seismically isolated buildings. The seismic input is defined by the computation of realistic synthetic seismograms considering different levels of detail for the earthquake source, both for a bedrock model and taking into account the specific site conditions. The horizontal response spectrum, calculated in the centre of the municipality by modelling the most dangerous source, advises against the construction of a building with a fixed base, but it is compatible with the seismic isolation, and it has been, therefore, used for the design of a residential seismically isolated building. The maximum displacement for the isolation system has been estimated about 17 cm, a value much lower than that provided by the code design response spectrum (28 cm). The importance of the realistic modelling, which permits the generalization of empirical observations by means of, physically sound, theoretical considerations, is evident, as it allows for the optimisation of the structural design with respect to the site of interest.     

有耗媒质中脉冲响应的研究

本文详细论述了一维有耗媒质中脉冲响应的正演理论及其模拟技术。深入研究了导电率、介电常数对脉冲响应的影响。以平面脉冲电磁波在一维介质中的传播理论为基础，编制脉冲响应的仿真模拟软件。最后，用实际模型验证了该模拟技术的正确性和可靠性。     

Stochastic Fracture Propagation Modelling for Enhanced Geothermal Systems

Fractures and fracture networks are the fundamental components of enhanced geothermal systems and determine their technical and economic viability. A realistic fracture model that can adequately describe a fracture-stimulated reservoir is critical for subsequent flow and heat transfer analyses of the system. Fractures in these systems are essentially the product of hydraulic stimulations of the reservoir that, together with ground conditions and the local stress regime, determine how fractures are formed and propagated. This paper describes three methods for generating realistic fracture models for enhanced geothermal systems; two of them incorporate the fracture propagation process in the modelling and hence provide a stochastic fracture propagation model. The methods are: a Bayesian framework in the form of Markov Chain Monte Carlo simulation, an extended Random Sampling Consensus model and a Point and Surface Association Consensus model. The conditioning data used in these methods are seismic events recorded during fracture stimulation. Geodynamics’ Habanero reservoir in the Cooper Basin of South Australia is used as a case study to test these methods.     

与城市地下生命线工程有关的双层地质结构中的几个问题

近代城市地下生命线工程要求对埋藏基岩中的活动断裂作出工程评价,以判断其是否有出现地表破裂的可能。 本文提出了对研究活动断裂的3项内容,即:“活动断裂的活动时间下限、地震复发周期及活动断裂的活动速率;对不同的断裂活动机制,即蠕动和粘滑,从覆盖层厚薄两种情况分析了对地基所产生的影响。 作者通过超声地震模型实验对覆盖层横向非均一性影响地下生命线工程的特点进行了研究,结果表明,横向非均一的覆盖层在地震过程中对地下生命线工程影响显著,不可忽视。     

埋藏基岩地形对地震波动影响的模拟研究

本文以超声地震模型试验为基础,着重从地震波的强度方面来阐明埋藏基岩地形对地震波动的影响,提出了利用振动能量分析地震动问题的重要性,理论分析及实验结果表明,埋藏基岩地形对地震波动的影响表现为:一方面基岩面凹凸不平,致使基岩面上不同部位接收能量大小不同;另一方面,基岩的凹凸面对地震波产生汇聚与发散作用,这两方面共同作用的结果使得地表不同地点表现出不同的地震动强度。     

Topography and Soil Effects in the MS 5.9 Parnitha (Athens) Earthquake: The Case of Adámes

Large concentration of damage to residential and industrial buildings occurred in regions near the banks of the Kifisos river canyon during the 7-September-1999 Parnitha (Athens) Earthquake. One such region, which experienced unexpectedly heavy damage, was the small community of Adámes, which borders the canyon near its deepest point. To explore whether in addition to structural factors the particular topographic relief and/or the actual soil profile contributed to the observed concentration and non-uniform distribution of damage within a 300 m zone from the edge of the canyon cliff, wave propagation analyses are conducted in one and two dimensions. Finite-element and spectral-element formulations are used to this end. To avoid spurious wave reflections at the artificial boundaries, ourtwo-dimensional (2-D) finite-element analyses utilize Bielak's Effective Seismic Excitation method. Soil layering and stiffnesses are determined from 10 SPT-boreholes and 4 crosshole tests. Ricker wavelets and six realistic accelerograms are used as excitation; two of the latter are selected from the literature and four are obtained on the basis of the four strongest motions of the earthquake, recorded in central Athens.The results show that the 2-D topography effects are substantial only within 50 meters from the canyon ridge. These effects materialize only in the presence of the relatively soft soil layers that exist in the profile at a shallow depth. The so-called Topographic Aggravation Factor (TAF), defined as the 2-D over 1-D Fourier spectral ratio, varies around 1.4 over a broad frequency band which covers the significant excitation frequencies. At the location of four collapsed buildings, about 250 m from the edge, 2-D (topography) effects are negligible, but the specific soil profiles amplify one-dimensionally all six ground base excitations to spectral acceleration levels that correlate well with the observed intensity of damage, at least in a qualitative sense.     

Electrical imaging of the groundwater aquifer at Banting,Selangor, Malaysia

A geophysical study was carried out in the Banting area of Malaysia to delineate groundwater aquifer and marine clay layer of the alluvial Quaternary deposits of Beruas and Gula Formations. The Beruas Formation is formed by peat and clayey materials as well as silt and sands, whereas the Gula Formation consists of clay, silt, sand and gravels. Both Formations were deposited on top of the Carboniferous shale of the Kenny Hill Formation. A 2-D geoelectrical resistivity technique was used. Resistivity measurement was carried out using an ABEM SAS 4000 Terrameter. The 2-D resistivity data of subsurface material for each survey line was calculated through inverse modelling and then compared with borehole data. The resistivity images of all the subsurface material below the survey lines show similar pattern of continuous structure of layering or layers with some lenses with resistivity ranging from 0.1 to 50 Ωm. The upper layer shows resistivity values ranging from 0.1 to 10 Ωm, representing a clay horizon with a thickness up to 45 m. The second layer with depth varies from 45 to 70 m below surface and has resistivity values ranging from 10 to 30 Ωm. Borehole data indicate coarse sand with some gravels for this layer, which is also the groundwater aquifer in the study area. The lowermost layer at a depth of 70 m below ground level shows resistivity values ranging from 30–50 Ωm and can be correlated with metasedimentary rocks consisting of shale and metaquartzite.     

青藏工程走廊多年冻土场地地震动加速度峰值特征研究

根据青藏工程走廊北麓河及楚玛尔河场地地震危险性分析结果,合成年超越概率为1.97%、1.00%、0.21%、0.10%、0.04%、0.02%的人造基岩地震波作为输入地震动,结合场地钻孔剖面及波速资料,和已有的冻土动力学研究成果,建立一维模型,通过等效线性化方法进行场地地震反应分析计算,研究了青藏工程走廊多年冻土场地地震动加速度峰值特征及影响因素.研究结果表明,北麓河场地与楚玛尔河场地的人造基岩地震波峰值及持时均存在显著差异,北麓河场地峰值大、持时短,以近震影响为主,楚玛尔河场地峰值小、持时长,以中远震影响为主;多年冻土区场地,夏季场地地震动加速度峰值显著大于冬季,活动层融化对场地地震动加速度峰值有明显的放大效应;冬季场地冻结后,场地地震动加速度峰值随冻土波速增大而减小,最大减小幅度为6.1%,随动剪切模量比减小、阻尼比增大而减小,最大减小幅度为8.9%.活动层的融化有利于放大场地地震动加速度峰值,重大冻土工程抗震设防应予以重视.     

Multichannel analysis of surface waves (MASW) for seismic site characterization using 2D genetic algorithm at Bahrah area,Wadi Fatima,Saudi Arabia

The present study deals with the seismic site classification of Bahrah area, Wadi Fatima, to characterize the local site conditions. The dynamic behavior of sediments was studied by the application of surface wave inversion. The multichannel analysis of surface waves (MASW) shallow geophysical technique was utilized for site classification. MASW survey was carried out at 66 sites along with 13 seismic refraction profiles at suitable localities. MASW and seismic refraction profiles were processed and compared with the available borehole data. The integration of MASW techniques with seismic refraction and borehole data progressively enhanced the subsurface visualization and reliability of the shear wave velocity estimation in the subsurface in the study area. The subsurface shear-wave velocity model was achieved by the solution of an inverse problem-based dispersion of surface waves and propagation in a vertically heterogeneous medium. The 2D genetic algorithm was employed for the inversion of dispersion curves to obtain velocity and thickness of subsurface layers. The depth to engineering bedrock and velocity of shear waves in the first 30 m was deciphered and mapped. The depth of bedrock in study area varies from 4 to 30 m, and V _{S 30} ranges from 320 to 800 m/s. The most of study area falls in B and C class categories in addition to few sites of D class according to the NEHRP guidelines.