论震源辐射问题

地震震源辐射过程是现代地震学中一个复杂而关键的课题,而震源谱的研究则是认识震源辐射的重要一环。地震波由震源传播至地球表面,经历了一系列诸如路径衰减、地壳表层放大和场地效应等物理作用。其中震波衰减效应包括几何扩散、非弹性衰减(YQ(f))和近地表高频衰减。地壳放大作用主要发生于表层或浅层波阻抗梯度带。为了更确切地描述和理解震源辐射谱,必须从观测到的地面运动记录中把地震源谱分离出来,从而消除传播路径和地表场地效应。强震运动记录是研究震源谱的基本资料。采用频率域方法,以傅里叶变换为工具,可使时间域的卷积问题简化为频率域的乘积运算。研究区域含日本、墨西哥、土尔其、加利福尼亚、加拿大西部(British Columbia)和北美东部(ENA)等典型构造区。结果表明,在适当消除路径和场地效应之后,震源谱的基本特征只随震级变化,而与研究地区无直接关系,亦即震源谱基本独立于构造区域、震源距和震源深度。这对于未来强震运动预测和地震灾害评估具有十分重要的理论和实际意义。与此同时,高频衰减因子(Kappa)与构造环境有关：低Kappa工资值相应于较稳定的板内构造环境下的硬岩场地(如北美东部),而相对较高的Kappa值则相应于比较活跃的构造环境下的场地条件,如日本、墨西哥、加拿大西海岸、美国西部的加利福尼亚和土耳其的转换构造带。强震运动水平与垂直分量的频谱比RH/V(f)作为频率的函数可近似描述为地壳放大和场地高频衰减的综合效应：即RH/V(f)=A(f)^-πkf。其中,A(f)是地壳表层放大函数,k是Kappa因子。通过震源谱的对比研究。提供了一个新的震级转换公式。     

地震震源与灾害评估分析:中国四川2008年5月12日Ms 8.0汶川地震

汶川5月12日8.0级地震在构造上起因于印度板块与欧亚板块以每年约5 cm的速度聚敛,并因此而引起青藏高原的地壳物质向四川盆地及中国东南大陆运移.主震震源及余震活动集中于以龙门山为中轴的一条长约350 km、宽约100 km的地震活动带.震源深度一般分布丁地壳脆性-韧性转换边界以上约10～20 km区间的地壳震源层之中,属浅源构造地震.主要震源机制与龙门山构造运动方式密切相关,以其地壳厚度向西急剧加厚、重力梯度带、高波速比(Vp/Vs～2.2)等深部异常及逆冲断层兼具走滑性质的地质构造为特征.在震源辐射、路径传播和场地效应研究的基础上,分别计算并比较了岩石和土壤条件下的地震响应谱,特别强调了土壤条件下的场地放大效应;同时对与地震安全性有关的一些问题如地质灾害、地震频谱设计、地震早期预警系统及中、长期至短期地震预报等进行了探讨;特别提供了一个由加权平均计算、以岩石条件下震波衰减模式为基础的地震频谱设计参考实例.地震构造与动力学研究可融人工程地质与环境工程等学科发展.经历汶川地震考验的一些新近设计和建设的工程项目可为今后改进工程建筑规范与标准提供重要而有益的参考.地震预报是当今一大难题,但需探索研究,不可懈怠.地震减灾与预防足目前比较切合实际的安全举措.     

基岩表面地脉动谱比研究

目前,场地影响评价方法可以分为理论法和经验法。地脉动水平分量与竖向分量谱比法属经验法。这种方法将地脉动水平分量与竖向分量谱比的最大值对应的频率和幅值分别看作是观测场地的卓越频率和放大因子。虽然该法在国内外很多工程实例中得到了理想的结果。但由于这种方法建立过程中应用了许多假设,所以很多地震工程学家对这种方法持怀疑态度。本文讨论了基岩地脉动水平分量与竖向分量谱比为1这一假定的合理性。并采用泊松弹性半空间模型,应用概率分析方法考虑地脉动的随机性,推导了地下体波斜入射弹性半空间自由表面形成地脉动的水平分量与竖向分量谱比表达式。通过计算论证体波斜入射泊松弹性半空间形成的地脉动水平分量与竖向分量谱比特性。结果表明：当地脉动主要成分为体波时,基岩表面地脉动水平分量与竖向分量谱比接近于1;当地脉动主要成分为面波时,基岩表面地脉动水平分量与竖向分量谱比介于0.54～0.79之间。对泊松基岩,该比值是0.68。     

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.     

Assessment of the stability of H/V spectral ratios from ambient noise and comparison with earthquake data in the Cologne area (Germany)

Situated in an active tectonic region the highly industrialised Cologne area (Germany) suffers from moderate-sized earthquakes. Our contribution to the mitigation of earthquake risk includes a microzonation study using ambient seismic noise and earthquake recordings from two field campaigns. An analysis of these data using the horizontal-to-vertical (H/V) as well as the classical spectral ratio (CSR) technique led to the following observations and conclusions: (1) The spatial variation in the thickness of the sedimentary cover is reliably retrieved using the fundamental resonance frequency estimated from the peak in the (H/V) ratio of ambient noise. (2) This fundamental resonance frequency is stable in time, but the amplitude of the peak (the amplification factor) is not. (3) The relative amplification variation of the H/V ratio in the area should therefore be checked systematically with repeated measurements before drawing conclusions about its significance. (4) The analysis of the H/V ratio of ambient noise provides the lower frequency bound, above which amplification may occur. (5) The shapes of the spectral ratios obtained by the different methods are generally in good agreement. However, the analysis of earthquake data shows that significant amplification of ground motion may also occur at frequencies higher than the fundamental one even when thick sediments are present.     

Estimation of Response Spectra in the Severely Damaged Area During the 1995 Hyogo-Ken Nambu Earthquake

We calculated theoretical acceleration response spectra at the ground surface in the region near the 1995 Hyogo-ken Nambu earthquake source based on the spectral moment method. To estimate earthquake motion on the ground surface, a formula of earthquake motion at base rock level was derived. The amplification effect of the ground was introduced by using multiple reflection theory. Theoretically estimated response spectrum were modified by the response spectra calculated using observed earthquake motions.     

Effect of surface wave inversion non-uniqueness on 1D seismic ground response analysis

Surface wave methods are increasingly being used for geotechnical site characterization. The methodology is based on the dispersive characteristic of Rayleigh waves in vertically heterogeneous medium. Experimental dispersion curve is inverted to obtain one-dimensional shear-wave velocity profile by inverse problem solution. Uncertainty associated with this surface wave inversion has drawn much attention. Inverse problem solution can provide different equivalent shear-wave velocity profiles, which may lead to different seismic site response analysis. In this study, a neighborhood algorithm has been used for inversion of dispersion curve to get a set of equivalent shear-wave velocity profiles. These equivalent velocity profiles are then used for 1D ground response analysis for different input motion record of the same earthquake at different epicentral distances. Results show significant variation in amplification spectrum in terms of maximum amplification as well as peak frequency. The extent of this uncertainty largely depends on the characteristics of the ground motion records at different epicentral distances. A linear variation is observed between mean coefficients of variation of amplification spectrum and epicentral distance of ground motion records. A gradual increase in mean value of peak frequency and peak amplification with the epicentral distance is also observed.     

断层穿越斜坡场地效应分析:以绵竹九龙镇山前斜坡为例

2008年汶川MS8.0强震中频繁出现断层穿越的地震滑坡,除发震断层的地震特性外,其自身场地效应也会影响斜坡动力响应,甚至加剧斜坡失稳。本文以汶川极震区绵竹九龙镇山前斜坡为典型实例,根据余震作用下斜坡不同高程实测地震记录及地脉动测试结果,通过单点谱比法(H/V),获得斜坡地震动加速度随高程的放大系数和地脉动的频谱特征曲线,通过曲线对比分析发现: 1)斜坡两次典型余震PGA放大系数随高程先减小后增大,呈明显的凹形特征,凹形部位位于断层位置,其PGA放大系数约为斜坡底部测点0.4~1.0,坡顶测点的PGA放大系数则达到1.0~2.0倍。由于PGA放大系数是在断层位置出现的明显拐点,从地质上表明了断层场地效应明显; (2)各点NS/UD谱比普遍大于EW/UD谱比。断裂区域卓越频率为低频1Hz,小于其他测点的频率,对应谱比最大值高达3.0~4.0,高于其他测点谱比最大值。说明地表破裂处岩土体松散破碎,导致断层处的卓越频率较低,近场余震传播过来的高频地震波被断层隔断,地震加速度放大系数在该处发生了衰减,场地效应显著。本研究有助于增强断层场地对斜坡动力响应影响的认识。     

青藏高原东缘—扬子特提斯构造域深部结构与地壳形变研究

青藏高原东缘和扬子西缘的构造带是中国特提斯构造域的重要组成部分,该构造域受欧亚板块与印度板块陆—陆碰撞、高原隆升、块体裂解或拼接挤压等强烈构造活动的影响,记录和保存了多期次的特提斯构造演化历史痕迹。同时,该研究区域也是中国西部地区地壳形变最强烈的地区之一,其浅表形变特征与深部构造之间存在怎样的关联和制约机制是目前国际地球科学的一个研究热点。本研究依据作者十多年来持续在该区域开展的地质—地球物理研究,通过深部地球物理多参数结构成像、沉积盆地分析、地壳形变和强震孕育机制等综合对比分析,发现在青藏高原东缘的下地壳存在低速和高泊松比异常带,该异常体与来自青藏高原上涌的软流圈热物质汇聚,导致从扬子西缘到青藏高原的下地壳和上地幔的深部结构发生显著变化。沿着龙门山断裂带,中、下地壳存在交叠相间的低速(高泊松比)和高速(低泊松比)区域,这些深部结构分布特征与地表形变及前陆盆地隆坳格局具有较好的一致性。基于上述认识,提出了青藏高原东缘—扬子板块的深部接触模式及其相应的盆山耦合关系,阐明了板块碰撞—耦合的深部动力学过程对剧烈地壳形变、盆地隆坳格局和强震诱发的制约关系。本研究成果将为深入认识青藏高原东缘高原急剧隆升、盆地基底结构与隆拗格局,以及强烈地壳形变的深部动力学机制提供参考信息。     

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.     

Estimation of high-frequency seismic wave radiation on fault plane of 2008 Wenchuan earthquake based on improved empirical Green’s function

A new strategy for inversion of high-frequency wave radiation condition on the fault plan is exhibited. One-dimensional source model of large earthquake was divided into subfaults, each subfault contains a series of subsources to express high-frequency wave radiation. Envelope of large earthquakes can be expressed as a root-mean-squared with combination of envelope attenuation relationship from all subsources. The envelope attenuation relationship is considered as the empirical Green’s function. Distribution of subsources is estimated by envelope inversion. According to this method, the high-frequency (>?1 Hz) wave radiation areas of 2008 Wenchuan earthquake are generally inverted by the differential evolution using acceleration data from 27 near-field stations, while acceleration waveforms of the Lushan earthquake from 43 near-field stations were utilized to create attenuation envelope. High-frequency waves radiated in: (1) surface rupture areas, including Yingxiu and Beichuan areas; (2) close to the boundaries of asperities near Yingxiu, Yuejiashan, Beichuan and Nanba areas; (3) within 30 km length near the fault northeastern tip; (4) around Qingchuan area.

     

基于包络的汶川大地震高频地震波辐射区域反演及近场加速度合成

合理地估计汶川破坏区域的地震动有助于地震灾害的研究.通过利用芦山地震记录建立的加速度包络衰减关系和汶川地震近场30个台站的加速度包络,基于线源模型,采用差分进化方法反演了汶川地震断层面上高频 (＞1 Hz) 辐射区域分布.结果表明:断层面上高频辐射分布很不均匀,辐射较强的区域主要位于:(1) 产生较大地表破裂的映秀、北川和南坝区域;(2) 映秀和北川等凹凸体的周边区域,包括震中东北侧60~90 km区域、北川和南坝东北侧30 km处;(3) 断层破裂停止的东北端约30 km长的区域.其中,破裂贯穿到地表的映秀、北川和南坝是低频和高频辐射都很强的区域.对于无观测记录场点,选择其临近且场地条件类似的台站加速度提取平稳随机过程,结合高频辐射分布和衰减关系得到的包络,合成了加速度时程,可为汶川地震结构震害分析提供地震动输入.      

汶川Ms8.0地震后龙门山断裂带地壳应力场及其构造意义

2008年5月12日在青藏高原东缘龙门山断裂带中段发生汶川8.0级特大地震。大震发生时释放应力并对震源区及外围构造应力场产生影响,受汶川地震断层破裂方式和强度空间差异性的影响,震后龙门山断裂带地壳应力场也应表现差异特征,至今鲜有针对该科学问题深入的分析和讨论。经过系统收集、梳理汶川地震后沿龙门山断裂带水压致裂地应力测量数据与2008年汶川地震中强余震序列震源机制解资料,对汶川地震后龙门山断裂带中上地壳构造应力场进行厘定,通过与震前构造应力场对比,深入探讨了汶川8.0级地震对龙门山断裂带地壳应力场的影响,进而对汶川震后应力调整过程及青藏高原东缘龙门山地区深部构造变形模式进行研究,研究结果表明:受汶川8.0级地震的影响,震后龙门山断裂带地壳构造应力场空间分布具有差异性,近地表至上地壳15 km深度范围,映秀—青川段最大主应力方向为北西西向、地应力状态为逆走滑型,青川东北部最大主应力方向偏转至北东东向、应力状态转变为走滑型;15~25km深度范围,龙门山断裂带最大主应力方向仍为北西—北西西向、应力状态以逆冲型为主。汶川8.0级地震后,龙门山断裂带中地壳北西西向逆冲挤压的构造应力特征进一步支持了青藏高原东缘龙门山地区东西两侧刚性块体碰撞挤压、逆冲推覆的动力学模式。     

龙门山多层分层拆离地壳结构:新构造变形与深部构造证据

2008年5.12汶川地震发生后,对于龙门山地壳结构及其与汶川地震的成生联系成为构造地质学研究极为关注的科学问题。然而,现有的多种龙门山地壳结构模式在综合解释表层构造变形及深部构造时与调查和探测资料均有不符。利用前人地球物理探测成果,结合穿过龙门山主要发震构造单元彭灌杂岩及雪隆包岩体的综合构造剖面,将地表构造与深部地壳探测资料结合进行了综合解释,认为龙门山逆冲构造带中的多重冲断推覆构造由约10 km深处的拆离断层分隔,因而应该只是浅层次构造变形的组合样式;在中、下地壳韧性流壳层的主导下,扬子地块基底被动?入并形成多层拆离的韧性流变构造组合。5.12汶川地震及余震是由于以彭灌杂岩和雪隆包岩体为代表的刚性体,在上部韧性流壳层前端的持续推挤作用下,发生破裂而形成的。     

龙门山构造带南段向西南延伸的遥感影像证据及地震地质意义

在系统总结活动断裂遥感影像解译标志的基础上,利用Landsat ETM、Google Earth及ASTER GDEM等影像资料,结合前人研究成果,重点分析了龙门山构造带南段主要活动断裂的空间展布及几何学与运动学特征。研究结果表明,龙门山构造带在向南延伸过程中发生了明显的断裂分散现象,整个断裂带逐渐展宽,主要包括5条断裂带且其中包含多条次级断裂,至最南端被北西向鲜水河左旋走滑断裂带阻挡。其中活动性较为明显的断裂自西向东主要有4条:泸定断裂、天全断裂、芦山断裂和大邑-名山断裂,前两者是北川-映秀断裂的南延分支,而后两者是安县-灌县断裂的南延部分。由于龙门山构造带南段的构造变形被分解至多条次级断裂上,导致单条断裂错断地表的活动迹象明显变弱,因此单条断裂的潜在发震频率和强度也将相应变小,但潜在震源区会更为分散。结合已有的地震地质资料认为,未来应注意泸定和雅安2个地区的地壳稳定性及未来强震危险性问题。      

华北盆地Lg波衰减及台站场地响应特征

华北盆地是中国大陆地震活跃区之一,通过地震波衰减及场地响应参数研究该区构造介质属性及台基属性对地震预测预报、灾害评估具有重要意义.基于Lg波谱比的联合反演方法是获得地震波衰减参数及场地响应的有效方法,通过随机重采样方法可以检验解的稳定性.使用华北盆地68个台站记录的2004—2008年的149次地震,震级M_L为1.7~5.3的震中距为100~600 km,按信噪比大于2的标准挑选有效垂向记录1 000多条,地震射线较好地覆盖了华北盆地38°N~41°N、114°E~120°E区域.采用2.60~3.65 km/s的速度窗截取Lg波形并转化为频谱,研究频率范围为1~7 Hz,频率间隔0.2 Hz.计算得到的地震波衰减品质因子Q(f)与频率f的关系可表示为Q(f)=125±4.4f0.86±0.03,研究区为低Q₀(对应频率1 Hz),高频率依赖性的构造活跃区.基岩台站对地震波没有表现出明显放大作用,黄土沉积台站低频端比高频端明显放大;场地响应波动较大台站其解的标准偏差也大,说明场地响应的不稳定性体现了台基属性的非稳定性特征.      

Deep seismic reflection profiling across the Ou Backbone range, northern Honshu Island, Japan

Knowledge of the crustal structure, especially the geometry of seismogenic faults, is key to understanding active tectonic processes and assessing the size and frequency of future earthquakes. To reveal the relationship between crustal structure and earthquake activity in northern Honshu Island, common midpoint (CMP) deep reflection profiling and earthquake observations by densely deployed seismic stations were carried out across the active reverse faults that bound the Ou Backbone range. The 40-km-long CMP profiles portray a relatively simple fault geometry within the seismogenic layer. The reverse faults merge at a midcrustal detachment just below the base of the seismogenic layer, producing a pop-up structure that forms the Ou Backbone range. The top of the reflective middle to lower crust (4.5 s in travel time (TWT)) nearly coincides with the bottom of seismogenic layer. The P-wave velocity structure and surface geology suggest that the bounding faults are Miocene normal faults that have been reactivated as reverse faults.     

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

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

Lateral variations in shear velocity and attenuation in the upper mantle

We have used digitally recorded long period seismic data from the HGLP and SRO networks to obtain estimates of multiple ScS travel times and spectra and have used these to place constraints on the lateral variations in upper-mantle shear velocity and attenuation structure. The travel time results are consistent with our previous findings: to a first approximation, upper-mantle heterogeneity as expressed in vertical travel time differences can be correlated with surface tectonic features. A spectral stacking procedure for retrieving the attenuation operator for SH polarized multiple ScS waves has been applied to seismograms representing a variety of tectonic provinces. The spectral analysis has yielded stable estimates of the phase and amplitude response of the attenuation operator in the frequency interval 6–60 mHz. These estimates indicate that the apparent Q can also be correlated with surface tectonic regimes.     

PRELIMINARY RESULTS OF EARTHQUAKE MOMENT TENSOR INVERSION IN WESTERN CHINA

PRELIMINARY RESULTS OF EARTHQUAKE MOMENT TENSOR INVERSION IN WESTERN CHINA