用集集主震记录研究近断层强震记录的基线校正方法

从集集主震记录中选取有代表性的246个台站的记录作为数据集。在总结国内外基线校正方法的基础上,提出以下校正步骤:(1)初步判定记录中是否包含永久位移信息以便选用不同的校正方法。(2)以积分位移时程的末端平稳为标准,来确定积分速度时程末端拟合出的最优直线斜率作为校正量,去校正原始加速度时程,从而保留永久位移信息。永久位移等值线结果与GPS测量结果大体相符。根据校正结果推测认为:瞬变后保持或扩张或部分回弹,可能是决定永久位移形态的主要物理机制。     

近场数字强震仪记录误差分析与零线校正方法

本文对数字强震仪记录进行了误差分析,并对国家强震动台网入网的5种型号数字强震仪系统作了振动台对比试验,分析了该系统加速度记录积分后速度和位移时程零线漂移的原因.本文提出了加速度记录的零线漂移校正方法和校正准则.为了印证零线校正方法的可靠性,对振动台试验中强震仪记录到的加速度两次积分得出位移时程与试验时记录到的绝对位移进行比较,计算位移和振动台绝对位移完全一致;对2008年5.12汶川8.0级大地震和1999年台湾9.21集集7.6级地震现场加速度记录两次积分后得出永久位移与两次大地震的GPS同震位移进行比较.结果表明,该方法对大地震近场仪器墩会发生倾斜或产生永久位移时加速度记录的零线校正有明显效果,可以给出加速度积分后的速度和位移并符合校正准则.本文方法解决了对大地震近场地面运动的研究停留在对峰值加速度和反应谱的研究阶段的困惑,满足了结构抗震对地面永久位移的需求.     

基于近场强震记录和同震位移的汶川地震破裂过程

基于方位角覆盖较均匀的近场宽频带强震记录和同震位移资料,采用非负最小二乘法和多时间视窗技术,反演了2008年汶川地震破裂过程。从运动学的角度,对北川断层南段、彭灌断层和小鱼洞断层可能的破裂顺序进行了分析,同时探讨了北川断层南段浅部高倾角部分起始破裂的时间。得到如下结论:1)通过近场记录波形拟合残差结合同震位移资料得到,只有北川断层南段在与小鱼洞断层相交处发生双侧破裂,断层西南侧台站合成记录才符合观测记录,同时不会在彭灌断层南段近地表处产生远大于观测值的位错。得到1种可能的破裂方式为:主震从北川断层南段深部低倾角部分开始,造成浅部的彭灌断层破裂,彭灌断层往北侧破裂传播触发小鱼洞断层,进而引起北川断层南段浅部高倾角部分与小鱼洞断层相交处发生双侧破裂。2)通过分析断层西南侧台站合成记录的第2个波包波形拟合的结果,发现北川断层南段高倾角部分起始破裂可能有8s的"停滞"。从时间上来说,相互平行的北川断层南段和彭灌断层在破裂时可能存在先后顺序,北川断层南段浅部区域的破裂可能滞后于彭灌断层。同时,北川断层南段可能存在多点破裂的情况。3)断层面上滑动速度较大的区域,滑动速率的值与周边台站有较好的对应关系。断层面上产生较大滑动速率的区域,其周边台站往往有较大的PGV产生。     

考虑动态临界加速度的地震边坡永久位移预测模型研究

在区域边坡地震危险性评价中主要采用永久位移预测模型进行地震边坡永久位移计算.永久位移预测模型以Newmark滑块理论为基础,通过大量实测地震时程记录统计拟合得出.针对Newmark理论中滑动面抗剪强度参数保持不变和已有位移预测模型的计算位移小于实测位移的问题,利用动态临界加速度理论,分别构建含有峰值加速度和阿里亚斯强度...     

断层地表同震位错的平均位移估计方法

根据管道性能设计的理念,现行管道抗震设计规范规定现行输油气管道重要区段和一般区段的设防断层位移分别为预测的最大位移值和平均位移值,目前的断层位移估计方法一般给出的只是断层未来一百年的最大位移值,而不能给出平均位移值。针对此问题,本文综合国内外发震断层在地表的同震位错分布调查结果,总结分析了断层同震位错沿断层地表破裂走向的分布形态特点,并通过对实际震例中不同计算方法得到的平均位移进行对比分析,得到了断层的平均同震位错量与断层最大同震位错量的比值关系。基于国内外地震断层位移的大量数据,考虑一定的安全系数,本文建议采用预测断层最大位移值的2/3作为一般区段输油气管道的设防断层位移。      

基于第四代地震区划的跨越发震断层永久位移概率分析方法

本文主要借鉴了Lee(2007)和Youngs(2003)提出的地表永久位移估计方法,联合中国第四代地震区划图的潜在震源区参数,给出了适合中国跨越发震断层永久位移概率分析方法。震害调查表明,近场的结构破坏不仅是地震动引起的,地震地表的破裂或永久位移引起的破坏也占很大比重。应用与概率估计地震动相似的方法,给出不同概率水准的地震地表永久位移,使之服务于穿越地震发震断层管线、桥梁工程的抗震设防,采取有效的抗御永久地表位移措施,减小地震灾害带来的损失。     

GPS测得的汶川大地震同震位移

2008年汶川8.0级地震是建国以来破坏最惨重的大地震,是发生在地震前后GPS观测资料最多,震中在连续观测网中部及非连续观测站相当密集的地区.许多GPS连续观测站测得汶川地震同震位移.同震位移具有突变性和时间上的同步性以及变化幅度的显著性,是最确凿的与地震直接相关的地壳运动现象.震区外GPS连续观测站位移时间序列表明,此次大地震同震水平位移范围大,远离震中的许多地区观测到同震水平位移,与震前位移对比表现为弹性回跳.GPS连续观测站震前水平位移和同震水平位移过程揭示了此次地震震前大范围有显著变化的观测站水平位移与地震孕育过程的联系.但华北及邻近地区无明显同震水平位移.震中区外均未观测到明显的同震垂直位移.震区GPS站观测到的同震位移则主要为永久形变,不仅有大幅度的水平位移,也有幅度相对较小的垂直位移.本文研究了汶川地震同震位移的特点与机制,并由此进一步讨论此次大地震的成因、前兆地壳运动特征及其复杂性.     

基于同震位移的隧道抗震设计及其工程应用

首次在隧道抗震设计中考虑地震同震位移场,通过搜集历史地震破裂信息估计工程区地震位错,采用"分裂节点法"获得隧道周边的同震位移场,并联合拟静力隧道抗震设计方法,成功应用于某公路隧道。研究表明:基于同震位移的隧道抗震设计方法能够为我们估算未来可能的地震灾害提供有力的定量分析工具;当地层自振周期和基岩反应谱未知,可以利用同震位移场估算远场剪应变开展隧道拟静力抗震设计;例如对于云南某隧道,如果龙陵—澜沧断裂带发生地震,隧道南北两端强制位移约为6 cm。如果龙陵—瑞丽断裂发生地震,由于强制位移导致隧道南北部分出现反向运动,将会对隧道中部产生剪切作用,且隧道出现明显的横向和垂向差异运动。     

2014年鲁甸6.5级地震GPS同震位移及反演分析

本文综合GPS流动和连续观测结果,并利用最速下降法（SDM）反演方法,给出并分析了2014年鲁甸6.5级地震同震位移、断层面滑动位移分布特征.GPS同震位移结果表明：此次地震沿北西方向表现出左旋应变释放特征、沿北东向表现出拉张应变释放的同震特征,并且随着离开断裂带距离的增加,拉张变形衰减;受包谷垴—小河断裂控制的左旋剪切应变释放的位移在莲峰、昭通—鲁甸断裂附近较弱,说明该断裂可能没有完全切割昭通—鲁甸断裂,不属于该区域主干断裂;昭通—鲁甸断裂带有一定的右旋应变释放,而逆冲应变释放不明显,表明该断裂带处于受南东向挤压的强闭锁状态.SDM反演结果表明,鲁甸地震以左旋走滑为主并兼有拉张性质,地震矩震级为M_W6.3左右.综上所述,并结合其他研究成果,我们认为莲峰、昭通—鲁甸断裂带仍存在强震危险性.

     

Approximate recovery of coseismic deformation from Taiwan strong-motion records

Since 1990, digital strong-motion accelerographs and global positioning system (GPS) instruments have been widely deployed in the Taiwan region (Shin et al. 2003; Yu et al. 2001). The 1999 Chi-Chi, Mw 7.6 earthquake and the 2003 Chengkung, Mw 6.8 earthquake were well recorded by both digital accelerographs and GPS instruments. These data offer a good opportunity to determine coseismic displacements from strong-motion records and to compare the results with those derived from GPS measurements. As noted by Boore (2001), a double integration of the acceleration data often leads to unreasonable results, and baseline corrections are therefore required in most cases before the integration. Based on the works of Iwan et al. (1985) and Boore (2001), we developed an improved method for baseline correction and validated it using an extensive set of data from shake-table tests of a known “step” displacement on 249 accelerographs. Our baseline correction method recovered about 97% of the actual displacement from the shake-table data. We then applied this baseline correction method to compute coseismic displacements from the strong-motion data of the Chi-Chi and Chengkung earthquakes. Our results agree favorably with the coseismic displacements determined by the GPS measurements at nearby sites. The ratio of seismic to geodetic displacement varies from 0.78 to 1.41, with an average of about 1.05.     

Modeling coseismic displacements resulting from SAR interferometry and GPS measurements during the 1997 Umbria-Marche seismic sequence

In this study we analyse coseismic GPS displacements and DInSAR data to constrain a dislocation model for the three largest earthquakes of the 1997 Umbria-Marche seismic sequence. The first two events, which occurred on September 26 at 00:33 GMT (Mw 5.7) and 09:40 GMT (Mw 6.0) respectively, are investigated using both GPS displacements and DInSAR interferograms. We discuss and compare the results of previous studies which separately modeled a smaller subset of geodetic data. We provide a dislocation model for these two earthquakes which fits well both GPS and DInSAR data and agrees with the results of seismological and geological investigations. The first event consists of a unilateral rupture towards the southeast with a uniform dislocation. The strike, rake and dip angles are those resulting from the CMT solution. The second event consists of an unilateral rupture towards the northwest and a variable slip distribution on the fault plane. The strike and the rake are consistent with the CMT solution, but the dip angle has been slightly modified to improve the simultaneous fit of GPS and DInSAR data. While the second rupture (09:40 GMT) arrived very close to the surface, the fit to geodetic data shows that the first rupture (00:33 GMT) is deeper (2 km), despite the more evident surface geological effects. The analysis of new SAR interferograms allows the identification of a 5–6 cm additional displacement caused by the October 3 (Mw 5.2) and 6 (Mw 5.4) seismic events.We use data from a new DInSAR interferogram to model the displacement field of the Sellano earthquake of October 14, 1997. For this event significant GPS measurements were not available. We tested two different fault plane geometries: a blind, planar fault (top depth = 2.4 km), and a curved (listric) fault reaching the surface. The two models provide a generally similar fit to the data, and show that most of the slip was released at depths greater than 2.4 km along a gently dipping (40^°–45^°) fault surface. They also show that a unilateral rupture does not allow fitting the interferometric fringes since there is evident surface deformation to the northwest of the hypocenter. Moreover, we suggest that the concentration of high residuals in the southern part of our uniform slip model may in fact indicate a certain slip variability in this area.We conclude that, despite the moderate magnitudes and the lack of significant surface faulting, the space geodetic data allowed to constrain dislocation models giving new insights in the rupture process of the three largest events of the sequence.     

A note on the L’Aquila earthquake of 6 April 2009: Permanent ground displacements obtained from strong-motion accelerograms

On 6 April 2009, 01:32 GMT, an M_w 6.3 earthquake hit the Abruzzi region of central Italy causing widespread damage in the City of L’Aquila and its nearby villages. The mainshock of this earthquake was recorded by 57 digital strong-motion instruments, four of which are located on the hanging wall of the Paganica Fault near L’Aquila. These stations are no more than 6 km from the epicentre. We use accelerometric data from these four stations to estimate permanent ground displacements caused by the mainshock. Our numerical results reveal south-east and downwards directed permanent co-seismic displacements which are in fair agreement with the outcomes of GPS and InSAR measurements reported in preliminary Istituto Nazionale di Geofisica e Vulcanologia (INGV) reports.     

2004年苏门答腊地震的几个断层滑动模型的全球同震位移对比

2004年苏门答腊大地震后,不同作者根据地震波和/或GPS观测,提出了不同的断层错动模型.在利用同震位移观测资料反演断层滑动模型时,由于使用半无限空间均匀介质模型或半无限空间分层介质模型,一般只能利用近场位移GPS观测约束,无法利用远场资料,这些模型有时差别颇大,如何区别这些模型的优劣是一个仍尚未解决的问题.本文采用等效体力有限元方法,在考虑地球球形和分层的条件下,对四个不同作者提供的2004年苏门答腊地震的断层滑动模型计算全球同震位移.由于采用了球形模型,所以不仅可以把四个模型的近场位移计算结果与GPS数据进行对比,而且可以把远场位移计算结果与GPS数据进行对比.我们发现,垂直位移对断层滑动模型的依赖性小于水平位移.四个模型计算的近场位移与GPS位移符合程度均较好,但是四个模型计算的远场位移与GPS位移符合情况有很大不同,其中Chlieh等（2007）模型在近场与远场符合程度均很好,是四个模型中最好的.另外还探讨了断层反演数据资料、断层几何模型以及地球模型对计算结果的影响.对于特大地震,全球同震位移观测与计算值吻合程度的好坏是衡量断层滑动模型的合理性的一个重要依据.     

芦山MS7.0地震前水平位移和同震水平位移研究

本文利用“中国地壳运动观测网络（二期）”多个GPS连续观测站观测数据处理结果,将2013年4月20日四川芦山MS7.0地震区域参考框架同震水平位移与全球参考框架同震水平位移进行比较,结果表明两组框架解一致,说明两种参考框架均可当作位错参考框架,也即全球参考框架同震水平位移也可视为区域参考框架同震水平位移．区域参考框架下GPS连续观测站地震前的水平位移和同震水平位移结果表明,震前数年,SCTQ站西侧的GPS站构造运动十分显著,而该站水平位移却很小,即出现反常的闭锁．但该站的同震水平位移使其弹性回跳至正常构造水平位移水平,因此SCTQ站震前的位移闭锁是水平位移空间分布中的异常,是芦山MS7.0地震的前兆．水平位移时空变化表明,该站震前和震时位移完全符合里德的弹性回跳理论．区域参考框架中位移时间系列和同震水平位移的综合研究有助于对芦山地震地壳运动前兆的认识和解释．尽管本文未能直接采用其它GPS连续观测站的资料,但结合本文和其它研究结果可以证实,震中附近其它站地震前后的变化与SCTQ站类似．基于芦山地震前水平位移和同震水平位移及其与前兆关系的研究,本文进一步讨论了GPS监测网的布设、 数据处理和分析等问题．      

On the reliability of long‐period response spectral ordinates from digital accelerograms

Using records from co‐located broadband and digital strong motion (SM) instruments, it is first shown that the displacement waveforms obtained by double integration of the accelerogram need not be free of unrealistic baseline drift to yield reliable spectral ordinates up to at least 10 s. Secondly, to provide objective criteria for selecting reliable digital SM records for ground motion predictions at long periods, a set of synthetic accelerograms contaminated by random long‐period noise has been used, and the difference between the original accelerograms and the spurious ones in terms of response spectra has been quantified, by introducing a noise index that can be easily calculated based on the velocity waveform of the record. The results of this study suggest that high‐pass filtering the digital acceleration record from a cutoff period selected to suppress baseline drifts on the displacement waveform appears to be in most cases too conservative and unduly depletes reliable information on long‐period spectral ordinates. Copyright © 2007 John Wiley & Sons, Ltd.     

基于GPS观测分析日本9.0级地震同震位错与近场形变特征

2011年3月11日日本本州宫城县东海岸近海发生MW9.0级地震,本文在对GPS同震位移场分布及误差特征分析的基础上,反演了同震位错分布.误差分析结果表明震源北西向300 km、北北西向550 km、南西向700 km范围内的同震位移量值明显大于误差,可以为位错反演提供有效的地表位移约束.沿震源北西向GPS剖面结果和位错反演位移剖面结果均表明同震近场位移符合指数衰减特征.位错反演结果表明,日本9.0级地震最大同震位错为25.8 m,位于震中附近;位错量大于10 m的同震破裂集中在震中附近400 km范围内;日本海沟南段同震位错量相对较小,此次地震为日本海沟地区典型逆冲型地震.根据此次9.0级地震和该地区以往强震破裂空间分布特征,此次9.0级地震破裂既体现了强震原地复发的特点,又体现了强震破裂的填空性.     

综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

Near-fault acceleration pulses and non-acceleration pulses: Effects on the inelastic displacement ratio

Near-fault ground motions can impose particularly high seismic demands on the structures due to the pulses that are typically observed in the velocity time-histories. The velocity pulses can be further categorized into either a distinct acceleration pulse (acc-pulse) or a succession of high-frequency, one-sided acceleration spikes (non-acc-pulse). The different characteristics of velocity pulses imply different frequency content of the ground motions, potentially causing different seismic effects on the structures. This study aims to investigate the characteristics of the two types of velocity pulses and their impacts on the inelastic displacement ratio (C_R) of single-degree-of-freedom systems. First, a new method that enables an automated classification of velocity pulses is used to compile a ground motion dataset which consists of 74 acc-pulses and 45 non-acc-pulses. Several intensity measures characterizing different seismological features are then compared using the two groups of records. Finally, the influences of acc-pulses and non-acc-pulses on the C_R spectra are studied; the effects of pulse period and hysteretic behavior are also considered. Results indicate that the characteristics of the two types of velocity pulses differ significantly, resulting in clearly distinct C_R spectral properties between acc-pulses and non-acc-pulses. Interestingly, mixing acc-pulses and non-acc-pulses can lead to local “bumps” that were found in the C_R spectral shape by previous studies. The findings of this study highlight the importance of distinguishing velocity pulses of different types when selecting near-fault ground motions for assessing the nonlinear dynamic response of structures.     

Coseismic and postseismic subsurface displacements and strains for a vertical strike-slip fault in a three-layer elastic medium

A three-layer elastic-gravitational fault displacement model using dislocation theory has been developed and used to examine the effect of layering of earth elastic moduli on surface and subsurface displacement fields for a vertical strike-slip fault. The model has been used to examine the effect of depth variation of elastic properties at coseismic and postseismic time scales. For pure strike-slip motion the effect of gravity on coseismic and postseismic horizontal deformation is negligible. For coseismic deformation the model predicts that (for constant Poisson's ratio) an increase in elastic moduli with depth attenuates the displacements within the upper layers with respect to displacement distribution for a uniform half-space, while an inclusion of a soft layer between the top layer and lower half-space amplifies upper layer displacements. The effect of variation in Poisson's ratio on surface and subsurface displacements has also been examined.The effect of postseismic stress relaxation on surface and subsurface displacements for a three-layer model has been calculated and compared with that of a uniformly relaxed half-space model. Layer 1 is assumed to correspond to the upper crust, layer 2 the lower crust and layer 3 the upper mantle. The effect of postseismic stress relaxation within a uniform half-space and within just the lower crust and upper mantle has been examined. Stress relaxation within the whole half-space decreases the amplitude and shortens the wavelength of displacements, while stress relaxation within the lower two layers increases the amplitude and broadens the wavelength of displacements. The difference between uniform and layered postseismic relaxation is particularly pronounced at the base of the crust.Coseismic and postseismic normal and volumetric strains for a vertical strike-slip fault have also been examined. For a uniformly relaxed half-space model, an increase in normal strains is shown with respect to the coseismic elastic solution, whereas the postseismic volumetric strain is effectively zero. For a three-layer model with stress relaxation in the lower layers only, the normal and volumetric strains within the top elastic layer resemble coseismic strains, while in the lower layers which suffer a rigidity decrease, the postseismic volumetric strain is effectively zero.