近断层地震动对大型地下洞室群地震响应的影响研究

针对现行水工抗震规范中对近断层地震缺乏考虑这一实际问题,首先讨论了近断层、远场地震动的差别,对比了近断层、远场地震动对地下洞室地震响应的影响,继而提出一种新的人工合成近断层地震动时程的方法,并采用合成的近断层地震动时程对某水电站地下洞室群进行了近断层地震稳定性专门研究。结果表明：与常规地震分析中考虑的远场地震动相比,近断层脉冲型地震动的主要特征为Vmax/Amax和Dmax/Amax（Vmax、Dmax、Amax分别为速度峰值、位移峰值、加速度峰值）指标更大,且近断层脉冲型地震动的主要能量集中在1 Hz以下的频段;在相同峰值加速度幅值、相同反应谱的前提下,近断层地震动对高边墙、大跨度地下工程的破坏程度远比远场地震动大;较之前人的研究方法,文中提出的近断层脉冲型地震动合成方法考虑了[1/Tp, 1] Hz这一频段内的地震波信息,更具合理性;在设定的计算条件下,地下厂房洞室群有在近断层地震动作用下发生失稳的危险性,需要进一步针对近断层地震动进行专项风险性评价,以讨论是否需要增设针对性的抗震支护措施。     

基于间接边界元法的近断层沉积谷地地震动模拟

将间接边界元方法拓展到近断层复杂场地地震动模拟求解,通过建立有限断层二维运动学模型,量化分析了走滑型断层逐次破裂错动下近场沉积谷地地震动放大效应。首先,将断层上、下盘以及沉积谷地分为不同的计算域;其次,在频域内通过利用不同交界面的应力、位移边界条件来建立边界积分方程,离散求解得到虚拟荷载密度,进而得到频域内场地的地震反应,时域结果由傅里叶变换得到;然后,通过与解析解对比验证了该方法的准确性;最后,分析了近断层沉积谷地地震动放大规律;研究了断层上界埋深、断层倾角、沉积谷地断层距以及断层单元的破裂速度等参数变化对场地地震反应的影响。结果表明:沉积谷地对近断层地震动幅值具有明显的放大效应,所分析模型加速度反应谱峰值可放大4.64倍;在沉积谷地内部,地震动持时明显延长,同时会出现幅值更大的长周期速度大脉冲;近断层地震动震害具有典型的集中效应,断层距增大10 km,半空间地表位移衰减大于50%。为近断层复杂场地地震动模拟提供了一种新的有效方法,对近断层复杂场地地震区划和工程结构抗震设计具有重要意义。     

近断层地震动模拟现状

地震动是由3个物理过程（震源破裂过程、波传播过程、场地反应）组成的一种复杂系统的产物,地震动模拟均是围绕这3个物理过程的建模开展的。地震动模拟目前仍然是一门相对较新的科学,强震观测中不断发现的新情况、新问题及其深入研究进一步推动近断层地震学理论和实践的发展。减少建模中的不确定性,用基于观测物理学的统计特征逐渐取代基于现象的假设描述,以改善地震动模拟的精度。基于大量地震动模拟的研究文献和资料,归纳、评述了近断层地震动模拟方法的现状、3个物理过程的建模方法及其发展趋势。     

大连开发区场地地震设计参数

根据工作区潜在震源区划分、潜在震源区地震活动性参数和地震动衰减关系, 进行场地地震危险性分析, 得到不同概率水平下场区相应地震烈度和基岩水平加速度峰值及其反应谱。根据场地工程地震条件划分不同地质单元及相应的场地类别, 进行不同概率水平的地震反应分析计算, 确定地震动设计参数。     

场地条件对地震动和震害影响的研究进展与建议

场地条件对地震动和震害有重要影响。本文简要回顾了场地条件对地震动和震害影响的研究历史,总结和评述了这一领域的研究进展。重点对地形地貌、岩土类型、覆盖层厚度、土层结构、地下水、岩土动力性质及物理地质现象等对地震动和震害的影响进行了总结和评述。提炼和概括了目前场地条件对地震动和震害影响研究中比较一致的几点认识。在此基础上,提出了这一领域当前应进一步研究的问题,主要包括地震现场资料调查、局部地形几何尺度、场地分类、覆盖层安全厚度界限、近断层地震动参数分布规律、活断层避让距离,以及强震资料积累等关键问题。     

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

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

集集地震的综合研究

本文收集、整理了台湾集集大地震前后福建、台湾地区数字地震观测、地震地质、大地形变测量以及板块运动的大量资料。对这些资料进行了系统的分析和研究 ,内容涉及集集大地震的序列特征、震源参数、震源机制、破裂过程、地壳形变、断层动力学、板块运动、地震成因以及前兆现象等。可供深入研究集集地震参考。     

龙门山断层破裂的频率非平稳特性

为研究地震过程中的频率非平稳特性,对近年来龙门山断层发生的两次大地震:汶川大地震和芦山大地震的近断层地震记录进行频谱分析。结果表明:相对于芦山地震有较大走向滑动分量的汶川地震,大多数位于汶川地震断层滑动前方的台站接收到更高的频率成分,位于断层滑动后方的台站接收到的地震波频率较低;尽管芦山地震断层相对汶川地震有较小的走向滑动分量,但仍然可以得出与汶川地震相同的结论,不同的是虽然芦山地震沿断层面向上方向分量大,但是其同一台站东西、南北、竖直三方向分量记录幅值相当。把芦山地震三分量记录变换到走向和沿断层面向上方向,证实了沿断层面向上方向高频成分更丰富。在断层滑动前方接收到的地震波频率较高,在断层滑动后方接收到的频率较低,这正是多普勒效应影响的结果。由于多普勒效应的客观存在,其对频率非平稳特性的影响与震源、传播路径和场地效应一样具有普遍性;所以,工程场地接收到的地震波的频率不仅取决于震源、传播路径、场地效应,还取决于断层滑动速度(多普勒效应)。     

近地表隐伏断层地震波场响应特征与初至波成像应用研究

近地表隐伏断层直接影响工程安全施工。通过构建浅覆盖层下方陡立界面断层数值模型,分析近地表隐伏断层地震波场的传播与转换规律。通过工程实测数据验证了数值模拟波场响应特征的可靠性,并利用初至波走时层析成像技术对地下介质波速进行成像。研究结果表明：(1)初至波转换形成的断点绕射波具有1频率升高、速度突变的特点,可清晰指示测区内断裂带的存在,可作为现场数据采集时断层的直观反映;(2)工程实测中通过原始波形到时、频率信息可初步判断断层的存在,基于射线追踪的初至波层析成像结果与后期钻孔揭露情况基本吻合,说明初至波层析成像在浅层地震勘探中的实用性;(3)接收测线排列内布设震源的方式,两侧震源延拓距离短,无需额外为震源激发提供场地,在工程施工中具有经济、快速、高效的特点,可作为野外复杂地表条件地下浅层结构探测的有效观测系统使用。     

考虑土体动力特征参数相关性的工程场地随机地震反应分析

研究等效线性化方法中土体动力参数不确定性对场地地震反应的影响,提出了可考虑土体动力特征参数间相关性的土体动剪切模量和动阻尼曲线随机样本生成方法,基于Matlab开发一维等效线性化场地地震反应分析计算程序,用于开展场地地震反应随机动力分析。以Ⅱ类成层场地为例建立一维自由场分析模型,采用不同地震设防水准下露头基岩场地加速度反应谱为目标反应谱,合成人工地震记录折半后作为场地下卧基岩处的输入。计算结果表明,考虑土体动力特征参数不确定性对场地地震反应具有较大的影响,且影响程度与地震动强度、频谱成分及场地基本周期均密切相关;场地最大峰值加速度和地表加速度反应谱的波动范围随地震动强度等级的增加而变大,且最大峰值应变和最大峰值加速度的波动范围可达10%和14%;目标反应谱平台段及场地基本周期处对应的场地地表加速度反应谱的波动范围超过20%。     

Seismic performance evaluation of dam-reservoir-foundation systems to near-fault ground motions

Ground motion records obtained in recent major strong earthquakes have provided evidence that ground motions recorded near the near-fault regions differ in many cases from those observed further away from the seismic source. As the forward directivity and fling effect characteristics of the near-fault ground motions, they have the potential to cause more considerable damage to structures during an earthquake. Therefore, understanding the influence of near-fault ground motions on the performance of structures is critical to mitigate damage and perform effective response. This paper presents results of a study aimed at evaluating the effects of near-fault and far-fault ground motions on seismic performance of concrete gravity dams including dam-reservoir-foundation interaction. Koyna gravity dam is selected as a numerical application. Four different near-fault ground motion records with an apparent velocity pulse are used in the analyses. The earthquake ground motions recorded at the same site from other events that the epicenter far away from the site are employed as the far-fault ground motions. The seismic performance evaluation method based on the demand-capacity ratio, the cumulative overstress duration and the spatial extent of overstressed regions is presented. The concrete damaged plasticity model including the strain hardening or softening behavior is employed in nonlinear analyses. Nonlinear seismic damage analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. The results obtained from the analyses show the effects of near-fault ground motions on seismic performance of concrete gravity dams and demonstrate the importance of considering the near-fault ground excitations.     

Selection and quantification of near-fault velocity pulses owing to source directivity

Design ground motions are typically governed by large earthquakes at close distances, but the increasing number of near-field recordings manifest the large variability in near-source ground-motion amplitudes which result in significant differences in the building response. In the near-field, this variability arises mainly from source-directivity effects, which generate strong near-fault velocity pulses. We investigate the effect of source complexity on the generation of velocity pulses using a geometrical approach to quantify directivity at near-fault sites based on kinematic rupture models. We propose selection algorithms that can be implemented as search strategies for directivity-related strong ground motion records which may serve as ground motion selection tools in large databases for structural analysis and liquefaction studies. We find that the existence of directivity pulses is strongly related to slip heterogeneity on the fault plane, i.e. that the location and size of asperities (large slip areas) determine directivity pulse generation. In this context we quantify several pulse properties, testing a variety of approaches, and develop predictive relationships between a number of source parameters and pulse properties. We find strong dependence of pulse period on total area of asperities, as well as on a geometrical directivity parameter. The empirical observations on velocity pulse generation determined by the proposed selection procedure are compared with the predictions using the geometrical directivity model. The results are important for determining the probability of observing a pulse at a site as a function of magnitude, distance, and slip heterogeneity on the fault plane.     

近断层速度脉冲地震动引起的边坡永久位移预测模型

基于小波分析方法,从NGA数据库的3 551条地震记录中选取189条速度脉冲地震动,地震动均转换成发生最强脉冲的方向。基于Newmark方法,分析了近断层速度脉冲地震动作用引起的边坡永久位移值。结果表明：近断层速度脉冲地震动对边坡产生特殊的破坏作用,表现在滑动位移值大、滑动体破坏力强等方面;边坡永久位移值与速度脉冲地震动的峰值速度具有高度相关性,位移值较大时尤为明显。建立了基于单变量形式的峰值速度及双变量形式的峰值速度、峰值加速度两种边坡永久位移预测模型,模型简单实用,与回归数据具有很好的相关性,前者更适用于预测对实际工程影响较大的永久位移值,且离散性较小。提出的预测模型为考虑近断层地震动速度脉冲特性影响的边坡永久位移值的概率地震灾害分析提供了基础。     

复杂艰险构造区铁路桥梁系统地震响应分析Part Ⅰ: 断层效应影响

近断层地震对桥梁的影响日益引起关注.本研究提出了桥梁?土?桩基全局建模方法,强调了更详细的桥墩及土壤非线性的真正好处,它可以比一系列轴载更真实地描述物理现象.协同SHAKE91程序并利用p-y曲线、t-z曲线和q-z曲线建立土?桩基非线性模型,采用双线性模型模拟桥墩及桩基础的滞回特性,建立不良地质发育区铁路桥梁?土?桩基多跨简支梁桥体系模型,计算其弹塑性地震响应,分析A_p/v_p等对桥梁的弹塑性地震响应的影响.研究结果表明:桥梁横竖向响应受A_p/v_p影响特点不同,相比墩底固结工况,考虑桩基后桥梁横向地震响应减小;对于竖向响应,在A_p/v_p > 10时桥梁竖向地震响应降低,说明竖向地震在较高频率影响桥梁地震响应.      

Comparison of near-fault and far-fault ground motion effects on geometrically nonlinear earthquake behavior of suspension bridges

This paper presents a comparison of near-fault and far-fault ground motion effects on geometrically nonlinear earthquake behavior of suspension bridges. Bo?azi?i (The First Bosporus) and Fatih Sultan Mehmet (Second Bosporus) suspension bridges built in Istanbul, Turkey, are selected as numerical examples. Both bridges have almost the same span. While Bo?azi?i Suspension Bridge has inclined hangers, Fatih Sultan Mehmet Suspension Bridge has vertical hangers. Geometric nonlinearity including P-delta effects from self-weight of the bridges is taken into account in the determination of the dynamic behavior of the suspension bridges for near-fault and far-fault ground motions. Near-fault and far-fault strong ground motion records, which have approximately identical peak ground accelerations, of 1999 Chi-Chi, 1999 Kocaeli, and 1979 Imperial Valley earthquakes are selected for the analyses. Displacements and internal forces of the bridges are determined using the finite element method including geometric nonlinearity. The displacements and internal forces obtained from the dynamic analyses of suspension bridges subjected to each fault effect are compared with each other. It is clearly seen that near-fault ground motions are more effective than far-fault ground motion on the displacements and internal forces such as bending moment, shear force and axial forces of the suspension bridges.     

Seismic demands and analysis of site effects in the Marmara region during the 1999 Kocaeli earthquake

The characteristics of the strong ground motion accelerograms from the 1999 Kocaeli earthquake are investigated in detail in this study. The emphasis is on the comparison of the response spectra for the fault normal (FN) and fault parallel (FP) components of the ground motions. The results show that the near-fault records with directivity effects characterize themselves with increased base shear demands rather than increased displacement demands for both the FN and FP components and a narrower velocity sensitive region for the FN component. This study also shows that the effectiveness of base isolation may vary from site to site and for a given site, from component to component. The site effects in the Marmara region during the 1999 Kocaeli earthquake are examined. Site amplifications are predicted by the classical spectral ratio (CSR) and the receiver function (RF) methods. The CSR method gives higher estimates for the site amplifications compared to the RF method and is in better conformity with the observed damage during the Kocaeli earthquake. The districts of Istanbul that are especially susceptible to site amplification hazard are determined. It is apparent from the results that the site amplification hazard risk is the highest for Avcilar and Bakirkoy districts. This study also shows that for sites which have the risk of soil amplification for long-period structures, liquefaction may not be beneficial as a natural base isolator, and may result in shifting the eigenperiod of the low- and mid-rise structures to the critical periods with high site amplifications. This may be especially the case for Avcilar and Bakirkoy districts. In Fatih, Bakirkoy, and Cekmece districts, the predominant period of the ground motion is calculated to be very close to the eigenperiods of the typical residential buildings. Therefore, these three districts are expected to experience heavier damages in future earthquakes due to resonance effects.     

Near-fault strong motion complexity of the 2000 Tottori earthquake (Japan) from a broadband source asperity model

The October 6/2000 Tottori earthquake that occurred in central Japan was an intermediate size strike-slip event that produced a very large number of near field strong motion recordings. The large amount of recorded data provides a unique opportunity for investigating a source asperity model of the Tottori earthquake that, combined with a hybrid strong motion simulation technique, is able to reproduce the observed broadband frequency near-fault ground motion.

We investigated the optimum source asperity parameters of the Tottori earthquake, by applying a Genetic Algorithm (GA) inversion scheme to optimise the fitting between simulated and observed response spectra and Peak Ground Acceleration (PGA) values. We constrained the initial model of our inversion by using the heterogeneous slip distribution obtained from a kinematic inversion of the source of previous studies. We used all the observed near-fault ground motions (−100 m) from the borehole strong motion network of Japan (KiK-Net), which are little affected by surficial geology (site effects).

The calculation of broadband frequency strong ground motion (0.1–10 Hz) is achieved by applying a hybrid technique that combines a deterministic simulation of the wave propagation for the low frequencies and a semi-stochastic modelling approach for the high frequencies. For the simulation of the high frequencies, we introduce a frequency-dependent radiation pattern model that efficiently removes the dependence of the pattern coefficient on the azimuth and take-off angle as the frequency increases. The good agreement between the observed and simulated broadband ground motions shows that our inversion procedure is successful in estimating the optimum asperity parameters of the Tottori earthquake and provides a good test for the strong ground motion simulation technique.

The ratio of background stress drop to average asperity stress drop from our inversion is nearly 50%, in agreement with the theoretical asperity model of Das and Kostrov [Das, S., Kostrov, B.V., 1986. Fracture of a single asperity on a finite fault: a model for weak earthquakes? Earthquake Source Mechanics, AGU, pp. 91–96.], and an empirical ratio of asperities to rupture area [Seismol. Res. Lett. 70 (1999) 59–80.].

The simulated radiation pattern is very complex for epicentral distances within half the fault length, but it approaches the radiation of a double-couple point source for larger distances.

The rupture velocity and rise time have a significant influence on the Peak Ground Velocity (PGV) distribution around the fault. An increase in rupture velocity produces a similar effect on the ground motion as a reduction in rise time.     

On the problem of the horizontal heterogeneity of the earth's crust and uppermost mantle in southern Eurasia

An analysis of the geological structure and history of development of the Himalayas points to the conclusion that prior to the Neogene most of that territory was part of the epi-Proterozoic Indian Platform; it was reworked only in the epoch of Neogene-Quaternary activation of tectonic movements. It was precisely during this period that the mountain ridge of the Himalayas formed. Only the narrow northern zone — the Himalayas of Tibet — which is composed of marine deposits of all systems of the Paleozoic and Mesozoic groups possibly did not develop on the platform with the Precambrian basement. Most likely, it is a near-fault folded zone that arose on the site of a near-fault trough, which already in the Cambrian Period had formed along the northern boundary of the Indian Platform.