基岩地震动的一个相干函数模型--倾滑断层情形

本文研究了基岩随机地震动的空间变化规律，考虑了震源破裂速度、子源个数、震源深度和介质传播速度等因素的影响。对应于每个样本，用数值模拟方法计算了采用震源位错模型的弹性半空间近场地震动场，最后通过统计方法给出了一个倾滑断层情形下的近场基岩地震动的相干函数模型。这一方法可以补充常用的统计方法因观测资料有限而导致的欠缺。     

1999年台湾集集大地震强地面运动的模拟

利用有限断层模型对1999年台湾集集大地震所产生的近场强地面运动进行了数值模拟,尝试了将具有明确物理意义的震源动力模型的研究结果应用于描述震源过程。计算结果显示,基岩台或接近工程基岩台的计算模拟记录同观测记录符合得很好。进一步研究表明,当综合考虑场地效应对地震动影响时,模拟计算结果得到了明显的改善,同实际观测记录更为相符。本文的研究结果说明,为了科学合理地预测大地震所产生的强地面运动,应综合考虑影响地震所产生的地面运动的"三要素",即:震源过程、传播路径及场地效应。     

近场强地震动数值模拟的简化计算方法

近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动。在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型。通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性     

断层附近地面地震动空间分布

运用震源位错模型, 分析矩形垂直断层及倾斜断层走向滑动和倾向滑动的近场地震动场，以地表地震动的傅立叶振幅谱比为参量考察断层附近地震动空间分布的特点. 结果表明, 断层附近的地震动强度主要受近旁子断层的控制，高强度的地震动分布在紧靠断层两侧有限的带状区域内，长周期分量受断层破裂传播方向性的影响. 走向滑动的方向性影响主要表现在垂直于断层走向的分量，倾向滑动则表现在平行于断层走向的分量，且深震在地面上引起的地震动强度分布比浅震要平缓，影响范围宽. 倾斜断层产生的地震动有明显的上盘效应，空间分布不对称，与观测结果相符. 最后给出了断层附近近场地震动强度分布拟合函数的表达式，并与美国的1997统一建筑规范规定的近场因子作了比较.      

汶川地震强地面运动模拟

基于确定性震源模型的方法主要用于计算低频(11 Hz)地震动常用经验格林函数法或随机方法,对低频地震动模拟不够准确．本文在确定性震源模型方法基础上,尝试采用分解给定的震源模型的方法来模拟宽频带(0.1——10 Hz)强地面运动,即采用分级离散断层面和分解断层面破裂单元上升时间的方法,增加震源时间函数中的高频信号,从而避免了对地震记录丰富程度和准确性的依赖．文中模拟计算了汶川MS8.0地震在8个地震观测台的地震动,将模拟结果和观测记录进行了加速度时程曲线和傅里叶振幅谱的对比分析．对比结果显示,模拟估计的地震动峰值加速度和持续时间与观测记录的数据基本在plusmn;50%的精度范围内相同,傅里叶振幅谱显示模拟结果有得到10 Hz左右的高频成分. 四川盆地中的台站模拟结果高频衰减比观测记录要快,原因是模拟过程没有考虑场地效应.对强地震动模拟还是要综合考虑震源、传播路径和场地的影响．研究结果表明,此改变震源输入的确定性方法可应用于模拟近断层宽频强地面运动．      

青藏工程走廊地震危险性对比分析及区划

潜在震源区地震活动性参数、地震动衰减关系对地震危险性分析结果至关重要。以中国第五代地震动区划图潜在震源区划分方案为基础,采用2类震级分档分别建立自编及五代图潜在震源区空间分布函数,收集4组青藏高原及周缘地区地震动衰减关系,采用不同组合对青藏工程走廊沿线的81个场点进行概率地震危险性分析计算,得到50年超越概率10%(地震重现期475年)的各场点基岩地震动峰值加速度(PGA),并转换为一般场地(Ⅱ类)PGA,对计算结果进行对比分析,并与第五代地震动区划图归档上下限值进行比较。结果显示:采用我国西部地区地震动衰减关系计算得到的PGA最大,采用云南地区地震动衰减关系得到的PGA最小,采用川藏地区及青藏高原东北缘地震动衰减关系时居中;在同一地震动衰减关系下,采用自编空间分布函数计算得到的PGA普遍略大于采用第五代图空间分布函数时;在震级上限为8.5的潜在震源区及附近地区,潜在震源区空间分布函数震级分档对计算结果有显著影响。综合分析表明,采用自编Ⅱ型震级分档空间分布函数方案与川藏地区地震动衰减关系组合方案的计算结果最为理想。最后,采用该组合方案对青藏工程走廊50年超越概率10%的基岩场地PGA及一般场地PGA进行了区划。     

场地条件对地震动相干函数的影响

本通过弹性半空间内位错源的数值解法研究了曲岩地震动相干函数，采用有限元方法分析了一些典型场地的地表地震动相干函数，两的对比结果表明：复杂场地对地震动相干函数的影响强烈。     

近场强地震动模拟中对破裂的方向性效应和上盘效应的表达.

为了模拟近场强地震动，采用了基于有限断层模型的一种随机合成方法．震源破裂面被剖分成一定数量的子源，总的地震矩分成数量更多的子震的矩．一给定子源中一次子震引起的场地地震动的傅氏谱，通过考虑点源的谱、随着距离的衰减、能量耗散及近地表效应等导出．据此幅值谱和一个随机相位谱，并与一个时程包络函数结合，合成一个子震时程．将各子源中各个子震引起的所有时程叠加，得出场地的地震动时程．叠加中，各子震时程之间的时滞，据子震发震时差和子源至场地的距离差别引起的时差确定．对一个设定地震，选用4个断层面倾角，计算了近场21个地点的地震动．结果表明，本文方法可以很好地表达破裂的方向性效应和上盘效应．为了验证方法的可靠性，对1994年美国北岭地震中3个近断层台站MCN，LV3和PCD模拟的地震动与实际记录的加速度反应谱和时程作了比较．      

新疆北天山地区宽频带地震动模拟关键技术研究

近年来,我国强震动观测记录数据库有了大幅扩充,一定程度上满足了重大工程抗震设防对宽频带地震动输入的迫切需求。然而,对于地震频发、强震动观测记录少的地区,宽频带地震动的数值模拟仍是确定合理地震动输入的重要途径。北天山地区是新疆经济社会发展的核心区域,现有强震动观测记录难以满足该区重大建设工程的抗震设防需求。研究构建考虑北天山地区震源机制、地质构造以及场地特征的宽频带地震动模拟方法对于高地震风险的北天山地区具有重要工程应用价值。宽频带地震动模拟方法由具有强区域特征的高频地震动随机模拟方法及基于求解无限域波动问题的低频地震动模拟方法组成。本文围绕北天山宽频带地震动模拟方法构建中的关键问题展开研究:针对高频地震动模拟,主要研究中小震震源模型的构建及非平稳地震动合成所需的群速度模型;针对低频地震动模拟,主要研究无限域波动数值模拟中高精度人工边界的构建,以尽可能缩小模拟工作区和提高计算效率。取得的研究成果如下:(1)北天山地区震源和路径传播模型构建。利用该地区2009—2015年843个M_L≥3.0地震记录,采用多台多地震联合反演方法建立了该地区的路径衰减、台站场地响应以及震源谱模型。结果表明,北天山地区中小震震源谱具有分段特性:4级以下地震的震源谱与引入震源谱高频衰减因子(κ=0.003)的Brune模型在统计意义上基本一致,4级以上地震的震源谱逐渐偏离Brune模型,呈现较为明显的双拐角频率特性;另外,震源谱存在分区差异,库拜盆地、塔里木盆地内地震的震源谱在高频段衰减更为显著。(2)北天山地区等效群速度模型构建。在已有基于等效群速度模型合成地震动的工作基础上,引入震相微分计算方法改进了等效群速度的计算方法;提出应用测震记录P波到时构建最快波群速度模型,分析了最快波群速度选取对等效群速度模型构建的影响;基于北天山地区的强震动观测记录,统计得到该地区的等效群速度模型。等效群速度模型引入到高频随机模拟方法的算例结果表明,基于等效群速度模型考虑地震动非平稳特征可给出更为合理的地震动模拟结果,对于面波发育的深厚覆盖层场地影响尤为突出。(3)高精度人工边界条件完美匹配层构建。完美匹配层内场方程和界面条件通常分别采用复坐标延伸技术变换强形式无限域内波动方程和界面条件得到,其构建过程相互独立,可能出现匹配不合理而引发的数值失稳、计算精度低下等问题。为此提出采用复坐标延伸技术变换弱形式无限域波动方程以构建完美匹配层的方法。弱形式波动方程耦合了波动方程及界面条件,规避了变换后所得场方程与界面条件的匹配不合理问题。新方法可直接给出弱形式匹配层,结合勒让德谱元建立了弹性介质近场波动谱元模拟方案。利用算例验证了新方案的精度及数值稳定性。(4)新源—和静6.6级地震地震动场的数值模拟。2012年6月30日新源—和静地震是北天山地区近50年来发生的最大地震,测震及强震动观测记录均较为丰富。本文应用反演所得滑动分布和随机滑动分布分别拟合了强震动台站的地震动时程,对比分析了拟合反应谱与观测记录反应谱的差异性;模拟了该地震的烈度分布图,分析其与实际调查烈度图和基于新疆常用地震烈度衰减关系计算的烈度图之间的差异,探讨了破坏性地震发生后应用高频随机有限断层方法生成地震烈度分布图的可行性。     

1604年福建泉州海外地震震级和震源参数的探讨——基于强地震动数值模拟方法

根据历史资料和已有的研究成果,利用两种地震动数值模拟方法对1604年福建泉州海外地震震级和震源参数进行了研究.首先利用经验公式确定了该次地震震级为71/2级和8级时的震源参数,然后基于离散波数法和有限断层地震动随机合成法,并利用基岩与一般场地水平峰值加速度的对应关系,计算了福建沿海8个场点的加速度、速度时程及其峰值,发...     

A coherency function model of ground motion at base rock corresponding to strike-slip fault

IntroductionEarthquakedamagesurveyandresearchresultshavedemonstratedthatspatialdistributiondifferenceofgroundmotionisoneoftheimportantreasonswhichcausedlongstructure(eglongspanbridge,undergroundpipe)destroy.Thathowtoprovideareasonableinputofgroundmotionfieldforaseismicdesignoflongstructureisaurgentprobleminearthquakeengineeringfield.Atpresent,themethodtostudyspatialvariationofgroundmotionsisadoptingstatisticanalysisbasedondensearrayrecordssuchasSMART-1array,etc,togetcoherencyfunctionofground…     

Strong-motion accelerograph array in Santiago,Chile, and preliminary evaluation of site effects

A strong-motion accelerograph array in Santiago, Chile has been installed. One of the sites is located on rock and the other six sites are on soil ground with different surface geology, so that local site effects on ground motions can be studied. As a preliminary evaluation of the site effects, the spectral ratios of weak-motion records at soil sites with respect to the rock site are calculated. The spectral ratios show that the amplification of ground motions with respect to the rock site is approximately 1.25 on dense gravel deposits, 2.5 on stiff pumice ground and 3.5 on soft silt ground.     

Ground Motion Modeling for Seismic Hazard Analysis in the Near-source Regime: An Asperity Model

—A new, yet simple, method using the asperity model to estimate ground motion in the near-source regime for probabilistic seismic hazard analyses is proposed in this study. This near-source model differs from conventional empirical attenuation equations. It correlates peak ground motions with the local contributing source in terms of the static stress drop released non-uniformly on the causative fault plane rather than with the whole seismic source in terms of magnitude. Here the model is simplified such that ground motions at a rock or firm soil site near extended vertical strike-slip faults are dominated by direct shear waves. The proposed model is tested by comparing its predictions with strong ground motion observations from the 1979 Imperial Valley and the 1984 Morgan Hill earthquakes. The results have revealed that ground motions in the near-source region can be adequately predicted using the asperity model with appropriate calibration factors. The directivity effect of ground motion in the near-source region is negligible for high-frequency accelerations. The cut-off frequency (?_max?) at a site is an important parameter in the near-source region. Higher values of ?_max yield higher estimates of peak ground accelerations. For high-frequency structures, ?_max should be carefully estimated. In the near- source region both non-uniform and uniform source models can produce non-stationary high-frequency ground motions. Peak motions may not be caused by the nearest sections of the fault (even if the uniform source model is considered).     

CHARACTERISTICS OF TORSIONAL GROUND MOTIONS

Due to the inherent difficulty in directly recording the rotational ground motions, torsional ground motions have to be estimated from the recorded spatially varying translational motions. In this paper, an empirical coherency function, which is based on the recorded motions at the SMART-1 array, is suggested to model the spatial variation of translational motions. Then, the torsional ground motion power spectral density function is derived. It depends on the translational motion power spectral density function and the coherency function. Both the empirical coherency function and the torsional motion power spectral density function are verified by the recorded motions at the SMART-1 array. The response spectra of the torsional motions are also estimated. Discussion on the relations between the torsional motion response spectrum and the corresponding translational motion response spectrum is made. Numerical results presented can be used to estimate the torsional ground motion power spectral density function and response spectrum.     

Testing the Validity of Simulated Strong Ground Motion from the Dynamic Rupture of a Finite Fault, by Using Empirical Equations

This paper is concerned with testing the validity of the ground motions estimated by combining a boundary integral equation method to simulate dynamic rupture along finite faults with a finite difference method to compute the subsequent wave propagation. The validation exercise is conducted by comparing the calculated ground motions at about 100 hypothetical stations surrounding the pure strike-slip and pure reverse faults with those estimated by recent ground motion estimation equations derived by regression analysis of observed strong-motion data. The validity of the ground motions with respect to their amplitude, frequency content and duration is examined. It is found that the numerical simulation method adopted leads to ground motions that are mainly compatible with the magnitude and distance dependence modelled by empirical equations but that the choice of a low stress drop leads to ground motions that are smaller than generally observed. In addition, the scatter in the simulated ground motions, for which a laterally homogeneous crust and standard rock site were used, is of the same order as the scatter in observed motions therefore, close to the fault, variations in source propagation likely contribute a significant proportion of the scatter in observed motions in comparison with travel-path and site effects.     

Stochastic model for simulation of near‐fault ground motions

A parameterized stochastic model of near‐fault ground motion in two orthogonal horizontal directions is developed. The major characteristics of recorded near‐fault ground motions are represented. These include near‐fault effects of directivity and fling step; temporal and spectral non‐stationarity; intensity, duration, and frequency content characteristics; directionality of components; and the natural variability of ground motions. Not all near‐fault ground motions contain a forward directivity pulse, even when the conditions for such a pulse are favorable. The proposed model accounts for both pulse‐like and non‐pulse‐like cases. The model is fitted to recorded near‐fault ground motions by matching important characteristics, thus generating an ‘observed’ set of model parameters for different earthquake source and site characteristics. A method to generate and post‐process synthetic motions for specified model parameters is also presented. Synthetic ground motion time series are generated using fitted parameter values. They are compared with corresponding recorded motions to validate the proposed model and simulation procedure. The use of synthetic motions in addition to or in place of recorded motions is desirable in performance‐based earthquake engineering applications, particularly when recorded motions are scarce or when they are unavailable for a specified design scenario. Copyright © 2016 John Wiley & Sons, Ltd.     

Modeling of spatially correlated,site-reflected,and nonstationary ground motions compatible with response spectrum

Seismic risk analysis and mitigation of spatially extended structures require the synthesis of spatially varying ground motions in the response history analysis of these structures. These synthetic motions are usually desired to be spatially correlated, site reflected, nonstationary, and compatible with target design response spectra. In this paper, a method is presented for simulating spatially varying ground motions considering the nonstationarity, local site effects, and compatibility of response spectra. The scheme for generating spatially varying and response spectra compatible ground motions is first established for spatial locations on the ground surface with varying site conditions. The design response spectrum is introduced as the “power” spectrum at the base rock. The site amplification approach is then derived based on the deterministic wave propagation theory, by assuming that the base rock motions consist of out-of-plane SH wave or in-plane combined P and SV waves propagating into the site with assumed incident angles, from which tri-directional spatial ground motions can be generated. The phase difference spectrum is employed to model ground motions exhibiting nonstationarity in both frequency and time domains with different site conditions. The proposed scheme is demonstrated with numerical examples.     

复杂场地条件下震源参数对断层附近长周期地震动的影响

基于显式有限元方法和运动学震源模型并利用昆明盆地三维地下构造模型,本文研究了震源参数对断层附近长周期地震动的影响.结果表明,断层的破裂方式、埋藏深度、破裂速度以及断层面上位错的不均匀分布对断层附近长周期地震动有重要影响.不同破裂方式下,破裂的方向性强的区域分布不同,由于破裂的方向性效应和复杂场地条件的共同作用,导致不同破裂方式的断层附近地震动分布差别很大.随着破裂速度的增加,方向性效应更加明显,断层附近的长周期地震动也随之增大;对于浅源地震,随着断层埋深的增加,地震动明显下降.对于埋藏深度很浅的断层,当Asperity靠近断层上沿时,会显著增大其在地表投影附近的长周期地震动.能否合理地估计这些基本震源参数,是预测未来发震断层周围地震动场的关键.     

Resonance and criticality measure of ground motions via probabilistic critical excitation method

Resonant characteristics of recorded ground motions are investigated and a new measure of criticality of ground motions is proposed. Four classes of recorded ground motions, i.e., (i) near fault motions (rock records), (ii) near fault motions (soil records), (iii) long duration motions (rock records) and (iv) long duration motions (soil records), are taken from Abrahamson N, Ashford S, Elgamal A, Kramer S, Seible F, Somerville P, Proc of First PEER Workshop on Characterization of Special Source Effects, 1998. It is shown that resonant characteristics of recorded ground motions can be captured appropriately by means of the probabilistic critical excitation method due to the present author regardless of the type of ground motions and the distance between the critical response and the actual one can be a new measure of criticality of ground motions. The time-averaged approximate treatment of nonstationary ground motions as stationary ones is shown to be adequate for structures with shorter natural periods subjected to long duration ground motions.     

The effects of torsion and motion coupling in site response estimation

Soil amplification characteristics are investigated using data from the Chibaken‐Toho‐Oki earthquake and its aftershocks recorded at Chiba dense array in Japan. The frequency‐dependent amplification function of soil is calculated using uphole‐to‐downhole spectral ratio analysis, considering the horizontal components of shear wave. The identified spectral ratios consistently demonstrate the splitting of peaks in their resonance frequencies and low amplification values in comparison with a 1D model. The torsional behaviour and horizontal ground motion coupling are clarified as the reasons for these phenomena at the site. To prove the hypothesis, the torsional motion is directly evaluated using the data of the horizontal dense array in different depths at the site. The comparison between Fourier spectra of torsional motion and identified transfer functions reveals the peaks at the same frequencies. The wave equation including torsion and horizontal motion coupling is introduced and solved for the layered media by applying wave propagation theory. Using the developed model, the effects of torsional motion with horizontal motion coupling on soil transfer function are numerically examined. Splitting and low amplification at resonance frequencies are confirmed by the results of numerical analysis. Furthermore, the ground motion in two horizontal directions at the site is simulated using site geotechnical specification and optimizing the model parameters. The simulated and recorded motions demonstrate good agreement that is used to validate the hypothesis. In addition, the spectral density of torsional ground motions are compared with the calculated one and found to be well predicted by the model. Finally, the results are used to explain the overestimation of damping in back‐calculation of dynamic soil properties using vertical array data in small strain level. Copyright © 2003 John Wiley & Sons, Ltd.