基于有限断层动力学模型的逆断层—三维层状沉积盆地地震动模拟

近断层效应使得沉积盆地对地震动放大效应更为复杂。本文针对逆断层发震下三维层状沉积盆地地震反应,基于波动谱元法,采用有限断层动力学模型,模拟断层动力破裂、地壳层地震波传播和层状沉积盆地对地震波散射全过程。在此基础上,对比分析了层状和均质沉积盆地对近断层地震动放大效应的影响,讨论了不同断层倾角下层状沉积盆地地震动加速度特性。结果表明：层状沉积盆地PGA空间分布与均质沉积盆地存在较大差异,由于近断层效应和盆地效应,层状沉积盆地地表局部范围竖向PGA大于水平向PGA;90°断层倾角下层状沉积盆地地表地震动放大范围与60°断层倾角结果明显不同,主要集中在盆地中心区域和断层附近,且幅值远小于60°断层倾角下结果;沿断层走向,盆地内地表地震动加速度峰值对应时刻较盆地外延后。     

Simulation of the Seismic Response of Sedimentary Basins with Vertical Constant-Gradient Velocity for Incident <Emphasis Type="Italic">SH</Emphasis> Waves

The simulation of the seismic response of heterogeneous sedimentary basins under incident plane waves is computed using the Indirect Boundary Element Method (IBEM). To deal with these kinds of basins we used approximate analytical expressions for the two-dimensional Greens functions of a medium with constant-gradient wave propagation velocity. On the other hand, for the homogeneous half space underlying the sedimentary basin, the full space Greens functions were used. The response of semi-circular heterogeneous basins under incident SH waves is explored by means of the displacements in the frequency-space diagrams and synthetic seismograms. Moreover, we compared these results with those obtained for other homogeneous semi-circular models. The principal differences among them are pointed out. This simulation provided interesting results that displayed a complex amplification pattern in a rich spectrum of frequencies and locations. The maximum amplitudes levels were found around the edges of the heterogeneous sedimentary basins. In time domain some features characterize the seismic response of the basin which include enhanced generation and trapping of surface waves inside the sediments, and the reduced emission of seismic energy to the hard rock. In the heterogeneous models the lateral reflections of surface waves greatly influence the total displacements at the free surface in comparison with the homogeneous models where the displacements have a shared influence among both vertical and lateral reflections.     

逆断层地震作用下三维沉积盆地地震动谱元法模拟

近断层地震动对近场区域会造成严重的地震灾害,尤其区域内存在沉积地形时会进一步加重灾害。目前,针对释放能量更大、破坏力更强的逆冲断层地震作用下沉积盆地的地震动响应研究还未见报道。本文即采用谱元法,研究了动力学逆断层地震作用下的三维沉积盆地的动力响应。文中以椭球形沉积盆地为例,对其逆断层地震动响应特性进行了分析,并探讨了改变内外介质波速比和沉积厚度时沉积盆地内部观测点地震动时程和峰值变化规律。研究表明:①沉积内外介质波速比对沉积盆地的地震动影响显著,当沉积内部介质波速比降低时,盆地内部地表的峰值响应增大,地震动持续时间明显延长,尤其是位于盆地中心加速度峰值放大2.08倍,持时延长1.97倍;②沉积盆地厚度同样对其地震动响应产生影响,当沉积厚度增加时盆地中心位置地震动响应减小,加速度峰值缩小约0.64倍,而盆地边缘区域的地震动响应明显增大,峰值放大约1.35倍。      

基于伪谱和有限差分混合方法的兰州盆地强地面运动二维数值模拟

利用伪谱和有限差分混合方法对兰州盆地进行二维强地面运动模拟.假设3个不同的震源深度5、10、20 km的同一地震下,以一个5层的二维剖面作为模型,以此来研究地震波传播过程并分析复杂的二维非均质路径效应,以及局部地下速度构造对强地面运动振幅加乘的影响,从而了解路径效应对强地面运动最大峰值位移的影响程度.研究结果显示:伪谱和有限差分混合方法模拟地震波场结合了有限差分法和伪谱法的优点,弥补了二者的不足,能较好地处理介质不连续面的计算,同时保证了和伪谱法相当的计算精度.通过地震波场模拟可知沉积盆地的强地面运动较基岩相比具有放大效应,震源深度的不同对兰州盆地垂直分量地面最大峰值位移影响较大,在震源深度为5 km时最小,为0.06 cm;震源深度10 km时在水平剖面30～55 km范围内峰值位移最大,达到0.14 cm;但对水平分量的峰值位移影响较小.     

Regional seismic responses of shallow basins incorporating site‐city interaction analyses on high‐rise building clusters

This study assesses the 3D amplification effects in shallow basins and quantifies the effects of site‐city interaction (SCI) on high‐rise buildings. A regional‐scale 3D spectral element simulation is conducted on the Tuen Mun‐Yuen Long basin, which contains multiple subbasins with heterogeneous and nonlinear soil profiles, while 3D city models with various building layouts are fully integrated into the basin model for our SCI study. We found a good correlation between spectral amplification factors and soil depths. Site response is significantly amplified at basin edges and centers due to surface waves generated at basin edges and the focusing effects stemming from 3D basin geometry. Transfer functions of 3D basins can be up to fourfold at fundamental frequencies as compared to 1D response, and further amplifications occur at high frequencies due to surface waves. In the SCI simulations, we observe wave trapping in the open space amid buildings resulting in energy concentration and up to twofold PGA amplifications. The wave trapping effect diminishes as the space between buildings increase beyond their range of influence (～100 m). The SCI analyses show that destructive kinetic energy in superstructures increases 28% in one horizontal direction but decreases 22% in the other. Our study concluded that, 1D site response analysis can significantly underestimate the seismic demand in shallow basins. Site‐city interaction of high‐rise buildings increases the short‐period spectra of ground motions, leading to an increase in their story accelerations by up to 50% and to a substantial decrease in the seismic safety of short structures in their vicinity.     

表层沉积物特性对地震地面运动的影响研究

建立包含震源、沉积盆地和表层低速沉积层的二维模型,采用交错网格有限差分/伪谱混合方法求解地震波传播,讨论沉积层厚度和速度对地震地面运动的作用。结果表明:沉积层内产生的地震波的多重反射以及转换会引起地面运动持续时间的延长,它们的相干叠加会造成地面运动峰值的放大;随着沉积层速度的增加,多重反射与转换波的能量减小,地面运动持续时间减小,但是不同速度或者不同厚度的低速层模型均显示出一致的地面运动峰值放大特征。结果说明,在包含震源、沉积盆地和沉积层的模型中,沉积层对地面运动的作用机理更复杂。在实际应用中有必要同时考虑这些因素的综合作用。     

2015年4月25日尼泊尔地震波场传播及烈度初步模拟分析

2015年4月25日尼泊尔发生Mw7.9地震.本文采用三维曲线坐标网格有限差分方法,依据USGS给出的震源运动学反演结果和震源区域地形数据,模拟了尼泊尔地震波场传播过程,得到震源区域烈度分布模拟结果,与实际观测的地表峰值速度(PGV)大体吻合.结果表明:地震烈度的空间分布整体上受控于震源的单边破裂特征,高烈度区域主要分布在破裂传播方向上,即震源东半部.震源南侧到东南侧近场,由于受到震源和地形的双重影响,形成最大烈度分布区域,最大烈度约为IX.南侧平原受低速沉积层影响形成高烈度区域.震源西侧及盆地内烈度相对较低.     

SH波入射下圆弧状沉积盆地地震反应瞬态解析解计算

地震波散射问题的解析解是研究局部场地、地形、盆地等不规则地层结构对地震动参数放大效应影响的重要理论工具.现有解析解大部分在频域内给出,无法直接用于研究不规则地层结构对地震动峰值、反应谱等参数的放大效应.本文基于平面SH波入射下圆弧状沉积盆地动力响应宽频带稳态解析解,通过Fourier变换,获取瞬态响应解析解.基于此,研...     

Effects of strong lateral discontinuity on ground motion characteristics and aggravation factor

This paper presents the effects of strong lateral discontinuity (SLD) in basins on the ground motion characteristics, differential ground motion (DGM) and aggravation factor (aggravation factor is simply the extra spectral amplification due to complex 2D site effects over the 1D response of the soil column). The seismic responses of open- and closed-basin models with SLD were simulated using SH-wave finite difference algorithm with fourth-order spatial accuracy. Simulated seismic responses, DGM and its spectra revealed that SLD induces Love waves, and the lower cutoff frequency for the same is equal to the fundamental frequency (f₀) of the soil layer. The maximum average aggravation factor (AAF) and DGM were obtained near the edge of open basin and at the centre of closed basin. A decrease of amplitude of Love wave, DGM level and AAF with offset from SLD was observed in an open basin. On the other hand, in closed basin, spatial variation of AAF and DGM level was highly variable. Duration of shaking, AAF and DGM level was more in the closed basin than in an open basin. Increase of DGM level and AAF with decrease of the width of basin was observed.     

面波频散与接收函数联合反演南北地震带北段壳幔速度结构

通过对南北地震带北段区域所布设的676个流动地震台站观测资料进行处理,联合反演面波频散与接收函数数据,获得了研究区内地壳厚度、沉积层厚度的分布情况以及地壳上地幔高分辨率S波速度结构成像结果.反演结果显示研究区地壳厚度从青藏高原东北缘向外总体逐渐变薄,秦岭造山带地壳厚度较同属青藏高原东北缘的北祁连块体明显减薄;鄂尔多斯盆地及河套盆地分布有非常厚的沉积层,阿拉善块体部分区域也有一定沉积层分布,沉积层与研究区内盆地位置较为一致;松潘—甘孜块体、北祁连造山带等青藏高原东北缘总体表现为S波低速异常;在中下地壳,松潘—甘孜块体下方的低速体比北祁连造山带下方的低速体S波速度值更小、分布深度更浅,更有可能对应于部分熔融的地壳;鄂尔多斯盆地在中下地壳以及上地幔内有着较大范围的高速异常一直延伸到120 km以下,而河套盆地地幔只在80 km以上部分有着高速异常的分布,此深度可能代表了河套盆地的岩石圈厚度,来自深部地幔的热物质上涌造成了该区域的岩石圈减薄;阿拉善块体在地壳和上地幔都表现出高低速共存的分布特征,暗示阿拉善块体西部岩石圈可能受青藏高原东北缘的挤压作用发生改造.     

Non‐uniqueness with refraction inversion – a syncline model study

Non‐uniqueness occurs with the 1D parametrization of refraction traveltime graphs in the vertical dimension and with the 2D lateral resolution of individual layers in the horizontal dimension. The most common source of non‐uniqueness is the inversion algorithm used to generate the starting model. This study applies 1D, 1.5D and 2D inversion algorithms to traveltime data for a syncline (2D) model, in order to generate starting models for wave path eikonal traveltime tomography. The 1D tau‐p algorithm produced a tomogram with an anticline rather than a syncline and an artefact with a high seismic velocity. The 2D generalized reciprocal method generated tomograms that accurately reproduced the syncline, together with narrow regions at the thalweg with seismic velocities that are less than and greater than the true seismic velocities as well as the true values. It is concluded that 2D inversion algorithms, which explicitly identify forward and reverse traveltime data, are required to generate useful starting models in the near‐surface where irregular refractors are common. The most likely tomogram can be selected as either the simplest model or with a priori information, such as head wave amplitudes. The determination of vertical velocity functions within individual layers is also subject to non‐uniqueness. Depths computed with vertical velocity gradients, which are the default with many tomography programs, are generally 50% greater than those computed with constant velocities for the same traveltime data. The average vertical velocity provides a more accurate measure of depth estimates, where it can be derived. Non‐uniqueness is a fundamental reality with the inversion of all near‐surface seismic refraction data. Unless specific measures are taken to explicitly address non‐uniqueness, then the production of a single refraction tomogram, which fits the traveltime data to sufficient accuracy, does not necessarily demonstrate that the result is either ‘correct’ or the most probable.     

浅硬场地剪切波速测试标准差对地震动影响初探

采用Shake 2000程序,以Turkey Flat试验场地为模型,通过输入不同类型下多种强度的地震波,计算研究多工况下剪切波速测试标准差对地表加速度反应谱和峰值加速度的影响。结论为:(1)浅硬场地上剪切波速测试标准差对地震动的影响很大,影响程度与输入地震波的强度和频率以及场地剪切波速计算值有关;(2)如果将反应谱残差大于20%或加速度峰值差别大于20%定义为统计意义上的不可忽略,那么剪切波速测试标准差对计算结果的影响在大多数情况下均不可忽略;(3)当输入波的卓越周期与场地特征周期接近时,浅硬场地上剪切波速测试标准差引起的反应谱变化非常显著;(4)只有当输入波的卓越周期与场地特征周期相差较大且输入波强度偏小时,剪切波速测试标准差引起的反应谱变化才可略去;(5)当浅硬场地上剪切波速实测结果低于统计均值时,地震动计算结果的偏差一般明显大于剪切波速实测结果,高于统计均值时引起的偏差,且地震输入越强表现越明显。     

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.     

OnP F and other phases in the early part of seismogram

Summary It has been found that when seismic energy propagates along the surface of the homogeneous crust beside usual Rayleigh waves, it produces certain instability in layers through which it propagates. In the light of this instability, a type of motion corresponding to longitudinal wave will be prominent in horizontal component compared to the vertical component; while transverse wave will be prominent in the vertical component but weak in the horizontal component, a contradiction with the existing knowledge. This has been identified withP _F phase. On taking the medium of propagation as slightly heterogeneous which allows existence of low velocity layer, a few larger number of such instabilities have been found. Velocity equation for Rayleigh waves for such media reveals existence of different velocities corresponding to vertical and horizontal components. Table for these velocities has been furnished.     

回填土不确定性对场地地震反应的影响

选取某核电场地控制性钻孔的厚度、剪切波速、密度等实际勘探数据,通过改变回填土剪切波速,分析了回填土不确定性对场地地震动参数的影响。研究结果表明:在回填土层厚度不变和模型总厚度不变的情况下,地表的水平向峰值加速度随着回填土剪切波速的增大而减小,但水平向峰值加速度增幅逐渐减小;回填土剪切波速到达一定的波速就不再影响地表水平峰值加速度;随着回填土剪切波速的增加,整个反应谱的谱值都普遍减小。     

Validation of 3-D basin structure models for long-period ground motion simulation in the Osaka basin, western Japan

We studied the applicability of two types of existing three-dimensional (3-D) basin velocity structure models of the Osaka basin, western Japan for long-period ground motion simulations. We synthesized long-period (3–20 s) ground motions in the Osaka basin during a M6.5 earthquake that occurred near the hypothetical Tonankai earthquake source area, approximately 200 km from Osaka. The simulations were performed using a 3-D finite-difference method with nonuniform staggered grids using the two basin velocity structure models. To study the ground motion characteristics inside the basin, we evaluated the wave field inside the basin using the transfer functions derived from the synthetics at the basin and a reference rock site outside the basin. The synthetic waveforms at the basin site were obtained by a convolution of the calculated transfer function and the observed waveform at the reference rock site. First, we estimated the appropriate Q values for the sediment layers. Assuming that the Q value depends on the S wave velocity V _S and period T, it was set to Q = (1/3V _S)(T ₀/T) where V _S is in m/s and the reference period T ₀ is 3.0 s. Second, we compared the synthetics and the observations using waveforms and pseudovelocity response spectra, together with a comparison of the velocity structures of the two basin models. We also introduced a goodness-of-fit factor to the pseudovelocity response spectra as an objective index. The synthetics of both the models reproduced the observations reasonably well at most of the stations in the central part the basin. At some stations, however, especially where the bedrock depth varies sharply, there were noticeable discrepancies in the simulation results of the models, and the synthetics did not accurately reproduce the observation. Our results indicate that the superiority of one model over the other cannot be determined and that an improvement in the basin velocity structure models based on simulation studies is required, especially along the basin edges. We also conclude that our transfer function method can be used to examine the applicability of the basin velocity structure models for long-period ground motion simulations.     

三维沉积盆地对地震动的放大效应间接边界元法模拟

基于一种高精度间接边界元法(IBEM), 实现了沉积盆地三维地震响应的频域、 时域精细求解, 并以半空间中椭球形沉积盆地对平面P波和SV波的散射为例, 着重探讨了入射角度、 入射波型、 入射频率、 盆地长宽比和深宽比对沉积盆地地震动放大效应的影响规律． 结果表明: 盆地形状对地震波的放大效应和空间分布状态具有显著影响, 且具体规律受控于入射波频段． ① 随着盆地深度增大, 盆地边缘面波发育更为充分, 在较宽频段内均会出现显著的地震动放大效应, 且深盆地的放大区域集中于盆地中部． ② 圆形盆地对地震波的汇聚效应最为显著, 而狭长盆地对地震波的汇聚作用相对较弱, 高频情况下可在盆地内部形成多个聚焦区域． ③ 不同波型入射下, 盆地对地震动放大效应的机制有所差异: P波入射下, 竖向位移放大主要是由于盆地边缘面波由四周向中部汇聚所致; SV波入射下, 边缘面波汇聚效应相对较弱, 而当盆地较深时, 底部透射体波和边缘面波易形成同相干涉从而显著放大地震动． 按盆地内外介质波速比为1/2, P波和SV波垂直入射下频域最大放大倍数分别为25和15, 时域放大倍数约为4.0和3.7(雷克子波)． ④ 低频波入射下, 位移从盆地中部向边缘逐渐减小, 且浅层沉积盆地对地表位移幅值的放大作用不明显． ⑤ P波和SV波的入射角度对盆地地震动放大幅值及空间分布特征也具有显著影响．      

岩土层速度结构回折波探测原理与方法

岩上层波速有随深度增加的总趋势。用浅层折射法记录到的折射波均属回折波。本文采用垂向速度梯度为不同正的常数多层水平层状介质模型，推导出回析波正反演公式，有简单解析表达式，据此解释的实测剖面，与验证孔波速测井资料相当吻合。表明在岩土层保持天然结构和横向比较均匀地区，回折波探测方法有良好应用前景。     

Basin and crustal velocity structure models for the simulation of strong ground motions in the Kinki area,Japan

We constructed a prototype of the basin and crustal structure model for the Kinki area, southwest of Japan, for the simulation of strong ground motions of hypothetical crustal and subduction earthquakes. We collected results of the deep seismic velocity profiles obtained by the reflection experiments and seismic imaging results, which were conducted in the Kinki area. The obtained profiles give underground velocity structures of the crust, from the surface to the subducting slab. We also gather the basin velocity structure information of the Osaka, Kyoto, Nara, and Ohmi basins. To examine the applicability of the constructed velocity structure model to the ground motion simulation, we simulated waveforms of an intermediate size event occurred near the source area of the hypothetical subduction earthquakes. Simulated ground motions using the basin and crustal velocity structure model are fairly well reproducing the observations at most of stations, and the constructed basin and crustal velocity structure model is applicable for the long-period ground motion simulations.     

Seismic dispersion and attenuation in layered porous rocks with fractures of varying orientations

Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C₃₃ (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C₃₃. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures.