Relative motions between cross-fault sites resultd from the 1999 MW7.6 Taiwan Chi-Chi earthquake
-
摘要:
为了获得近断层地震动准确的残余位移信息，提出了一种变步长网格搜索地震动基线校正方法。首先基于1999年我国台湾MW7.6集集地震强震动数据，探讨了近断层地震动的位移特征，然后通过地震动时程初始时刻的校正获得了集集地震多组上、下盘台站之间的相对运动时程数据，分析了分别靠近断层上、下盘场地之间相对运动的位移时程特征。结果显示：近断层地震动的滑冲效应显著，断层的平均滑冲速度可高达1.4 m/s；断层两盘间的相对运动类似于近断层地震动的运动特征，且幅值更大；近断层的残余位移约是断层两盘间最大相对位移的80%.
-
关键词:
- 1999年集集MW7.6地震 /
- 近断层地震动 /
- 跨断层场地 /
- 相对运动 /
- 断层作用
Abstract:In order to obtain more accurate residual displacement information generated by near-fault ground motion, this paper proposed a variable step-grid search method for ground motion baseline correction. The core of this method lies in preserving and accurately reflecting the residual displacement information of faults, and solving the problem of extremely low computational efficiency in traditional grid search methods when there are many digital seismic sampling points. Based on the strong ground motion data of the 1999 MW7.6 Chi-Chi earthquake in Taiwan region, this paper uses the variable step-grid search method for ground motion baseline correction to conduct in-depth research on the ground motion displacement characteristics near faults. By accurately correcting the initial moment of ground motion time history, this study successfully obtained relative ground motion time history data between multiple sets of the hanging and foot wall stations in the Chi-Chi earthquake. And based on this data, the displacement time history characteristics of the relative motion between the hanging and foot wall sites near the fault were analyzed. The research results indicate that the fling-step effect of ground motion is particularly significant in the Chi-Chi earthquake. The maximum displacement and residual displacement of near fault ground motion exceed 8 m and 6 m, respectively, and the average sliding velocity of the fault is as high as 1.4 m/s, fully demonstrating the strong destructive force of the earthquake. This poses challenges to the numerical simulation of fault rupture and related experimental researches. In addition, the relative motion between fault walls also exhibits similar motion characteristics to near fault ground motion, and its amplitude is larger. Therefore, the relative motion data between earthquake fault walls can refer to the time history of near fault ground motion as the input basis for engineering design and analysis. In addition, the residual displacement of the fault is about 80% of the maximum relative displacement between the fault walls. This study provides important information on residual displacement of fault ground motion, which is of great significance to the researches on fields such as earthquake engineering, earthquake warning, and disaster assessment.
-
引言
在日本Kobe地震(Kawase，1996)、美国Northridge地震(Graves et al，1998)和墨西哥Michoacan地震(Reinoso，Ordaz，1999)中均出现了沉积盆地内的震害比盆地周边地区更严重的现象，这是因为盆地内的沉积层对地震动的振幅和持时有明显的放大作用(Olsen et al，1995)，加重了对地表结构物的破坏(Olsen，2000；Komatitsch et al，2004；Shani-Kadmiel et al，2011)，因此要预测沉积盆地内的地震动，需要建立精确的盆地内沉积层三维速度结构模型(高孟潭等，2002；Kohler et al，2003；Aoi et al，2008).同时由于多数城市均位于沉积盆地内，这也使得沉积盆地三维速度结构建模的研究工作变得尤为重要(Koketsu et al，2009).
沉积盆地三维速度结构模型的准确度主要由建模数据的精度和数量所控制.近几年，一些研究人员基于丰富的钻孔、地震测线、重力勘探、地脉动、强震记录及地震波形等数据，在一些沉积盆地内建立了准确的三维速度结构模型，例如Koketsu等(2009)建立的Kanto盆地速度结构模型和Dhakal等(2011)建立的Tokachi盆地速度结构模型.在各类建模数据中，由于获取地脉动数据比获取其它建模数据更加经济和方便，使得地脉动数据被广泛地应用于地下速度结构的确定(Uebayashi et al，2012).目前采用地脉动数据建立地下速度结构模型的方法可分为两类：利用地脉动台阵法反演浅层剪切波速度(Arai，Tokimatsu，2005；Huang，Wu，2006；何正勤等，2007；王伟君等，2009；Claprood et al，2011；徐佩芬等，2012)和利用单点地脉动的水平向与垂直向H/V谱比反演地下速度结构(Delgado et al，2000；Dolenc，2005；Özalaybey et al，2011)，但通过这两种方法得到的速度结构均存在一定的误差(Arai，Tokimatsu，2004；Claprood, Asten，2010)，尤其是利用单点地脉动H/V谱反演得到的速度结构误差更大(王伟君等，2011).此外，单点地脉动H/V谱也常被用于评估场地的放大效应(陈棋福等，2008).
随着对地脉动研究的逐渐深入，Rhie和Dreger(2009)提出可通过对比理论地脉动谱与实测地脉动H/V谱的波形来修正沉积盆地的速度结构模型.但无论是基于一维速度结构模型推导出的地脉动H/V谱(Lunedei，Albarello，2010；Sánchez-Sesma et al，2011)，还是基于三维速度结构模型模拟出的地脉动H/V谱(Guéguen et al，2007；Barnaba et al，2010)，均与实测地脉动H/V谱的振幅存在偏差，所以基于对比地脉动H/V谱的波形变化对模型进行修正的方法并不成熟.此外，在模型修正环节中为避免模拟地脉动H/V谱，Suzuki等(2005)根据地震数据的H/V谱与基阶瑞雷波的H/V谱形状变化相似的原则对模型进行修正. Aoi等(2008)建立的Tokachi盆地三维速度结构模型以及Koketsu等(2009)建立的Kanto盆地三维速度结构模型，均利用了Suzuki等(2005)提出的模型修正方法，并得到了较精确的沉积盆地三维速度结构模型，证明了依据地震数据的H/V谱与基阶瑞雷波H/V谱形状变化相似的原则对模型进行修正的方法是可行的，但该修正方法依赖于丰富的地震数据，对于缺少地震记录的地区并不适用.在模型检验环节，Iwata等(2008)以及Iwaki和Iwata(2010)通过对比模拟的地震动波形与实测的地震动波形来检验沉积盆地模型建立的准确性，并取得了很好的应用效果，但该方法仅适用于强震记录丰富的地区.
鉴于此，为了避免以往建模过程中对地震数据，特别是强震数据的依赖, 本文拟提出一种包含数据预处理、模型建立、模型修正和模型检验的完整的建模方法.基于各类数据间不同的可信度，提出消除不同数据之间速度偏差的折减函数.为了避免以往在模型修正环节中对地震波形数据的依赖以及地脉动H/V谱的模拟等问题，本文提出一种改进的模型修正方法，即根据基阶瑞雷波H/V谱(H/V)R0与实测地脉动H/V谱(H/V)m形状变化相似(郭明珠等，2004；Arai，Tokimatsu，2004)的原则，对模型进行修正.由于地脉动数据的获取比较经济和方便，使得该模型修正方法具有广泛的适用性.本文以玉溪盆地为例进行建模，玉溪盆地位于云贵高原西南部的盆岭分布区，四周群山环绕，盆地内地形平坦，是一个典型的沉积盆地.
1. 建模过程
地震动模拟的关键是基岩面以上沉积层速度结构(Dhakal，Yamanaka，2013)，且基岩面以下缺少详细的建模数据，故本文仅建立该盆地内基岩面以上沉积层的三维速度结构模型.
在建立玉溪盆地模型的过程中，可用于盆地速度结构建模的数据有钻孔资料、地震勘探数据和地脉动数据.因为在地质分层较多的情况下，通过地脉动数据反演得到的地下速度结构存在较大误差，所以本文并未将地脉动数据用于初始建模，而仅用于模型修正.玉溪盆地的建模过程为：①数据预处理.依据经验公式(朱广生等，1995)将地震勘探得到的P波速度转换成S波速度.当钻井资料和地震勘探数据之间存在速度偏差时，考虑到各类数据间不同的可信度，对数据进行修正，以消除各类型数据间的速度偏差. ②模型建立.综合考虑钻孔资料和地震勘探数据，以剪切波速度为主要分层依据，同时参考地质序列.将钻孔资料和地震勘探数据作为模型层面的“硬性约束”，采用克里格插值方法，获得各沉积层面(各沉积层之间的分界面)的几何形状和展布形态，得到初步构建的盆地分层模型. ③模型修正.通过对比由盆地初始分层模型计算得到的基阶瑞雷波H/V谱(H/V)R0与实测地脉动H/V谱(H/V)m的波峰周期和波谷周期，对模型进行修正，得到最终模型. ④模型检验.用未参与模型修正的地脉动测点数据对修正后的模型进行检验，检查(H/V)R0与(H/V)m的波峰周期和波谷周期是否一致，进而确定本文所建立的速度结构模型的合理性.若检验结果显示不一致，则进一步修正模型.
2. 玉溪盆地三维速度结构模型
2.1 数据预处理
用于玉溪盆地速度结构建模的数据包括钻孔(深度为80 m左右)、地震勘探测线和地脉动，如图 1所示.由于地表附近(深度小于80 m)的钻孔资料与地震勘探数据的精度不同，造成即使在同一位置由地震勘探得到的剪切波速度远大于由钻孔得到的剪切波速度.
选取图 1中7个青色钻孔及其邻近地震勘探测线上距钻孔最近的7个点，计算这7组地震勘探数据与钻孔资料之间的剪切波速度偏差的平均值dv.使用R语言对剪切波速度偏差dv与深度x作相关性分析，得到两者的相关系数为-0.935，具有很好的负相关性，即地震勘探数据与钻孔数据之间的剪切波速度偏差dv随深度x的增加而递减.产生这种现象的原因是在通过钻孔实测数据获取剪切波速度数据的过程中不确定性因素较少，所以由钻孔实测数据获取剪切波速度的准确度较高.另一方面通过地震勘探数据获取剪切波速度的过程中不确定性因素较多，例如数据采集时外界环境的干扰、反演结果的不唯一性和用于反演的震相射线路径的平均效应等问题，均会造成由地震勘探数据获取剪切波速度的准确度较低；然而随着深度的增加，外界干扰相对降低，由地震勘探数据得到的剪切波速度的精度相对于地表处有所增加.
由于钻孔资料的精度较高，本文将由钻孔资料得到的剪切波速度视为准确的剪切波速度，使用剪切波速度偏差dv修正由地震勘探得到的剪切波速度.用R语言作剪切波速度偏差dv与深度x的多项式拟合(图 2)，得到一个速度偏差dv与深度x的三次多项式，即
![]()
图 2 地震勘探数据与钻孔数据之间的速度偏差dv (黑色圆点)随深度x的变化曲线Figure 2. Curve of velocity deviation dv of seismic exploration data from borehole data versus depth x
(1) 式(1)中各项系数(从3次项依次到0次项)t检验的概率值p依次为0.02349，0.00173，1.81×10－6和2×10－16，均小于0.025，因此拒绝各项系数为0的假设，即各项系数显著不为0.多项式的可决系数为0.9828，调整可决系数为0.9786.对多项式的整体显著性进行检验，即对残差进行F检验，得到p为7.408×10－11，小于0.01，拒绝速度偏差dv与深度x的精度存在显著性差异的假设，说明该拟合结果是可靠的.同时得到各项系数95%的置信区间分别为(-1.24×10-3, -1.08×10-4)，(5.44×10－2, 0.18)，(-10.13，-6.03)，(510.08，544.44).由式(1)可知，当深度约为131 m时dv减小至0.本文假定当深度大于131 m时由地震勘探得到的剪切波速度是准确的，可直接用于速度建模；当深度在131 m以内时，由地震勘探得到的剪切波速度需使用式(1)进行速度修正后，才可用于速度建模.
2.2 模型建立
何正勤等(2013)综合考虑地震勘探数据、钻孔资料和沉积层的地质序列，提出了玉溪盆地普渡河断裂带及其附近的速度结构模型.由于普渡河断裂带南段穿过玉溪盆地，所以该断裂带及其附近的速度结构能很好地反映盆地内的速度结构.本文在该速度结构模型的基础上，综合考虑已获得的各类数据，以剪切波速度为主要分层依据，初步假定玉溪盆地内存在6个沉积层，各沉积层的剪切波速度分别为400, 600, 800, 1200, 1400和1600 m/s.玉溪盆地的基岩主要为震旦系澄江组及中元古界昆阳群灰岩、砂岩和板岩(何正勤等，2013).本文将基岩面以下地层的剪切波速度定为3000 m/s(Koketsu et al，2009).
各沉积层面的形状变化满足空间自相关关系，因此可使用克里格插值方法对由钻孔资料和地震勘探测线得到的剪切波速度进行插值，得到玉溪盆地模型内各沉积层面的几何形状和展布形态，其中在克里格插值过程中所需的变差函数可通过拟合已知的剪切波速度数据来确定.
2.3 模型修正
目前计算理论地脉动H/V谱的方法仍不成熟，如果通过对比理论地脉动H/V谱与实测地脉动H/V谱的振幅变化的方法来修正模型，会使修正后的模型存在较大误差.本文根据基阶瑞雷波H/V谱(H/V)R0与实测地脉动H/V谱(H/V)m的变化形状相似(波峰周期和波谷周期相同，振幅不同)的原则，对模型进行修正.基阶瑞雷波H/V谱(H/V)R0为(Arai，Tokimatsu，2004)
(2) 
(3) 
(4) 式中：D为介质响应因子(Harkrider, 1964, 1970)；κ为波数，κ=(2/h)e－4πh，h为地层散射阻尼系数；PHR0(ω)和PVR0(ω)分别为瑞雷波水平向和垂直向功率谱；LV为当地脉动测点到地脉动加载源的距离大于瑞雷波的波长时地表处垂直分量的随机分布；θ为地脉动加载源水平向与垂直向之比，
(5) 
(6) 
(7) 式中: γn=2(βn/c)2，α和β分别为P波和S波的速度，c为瑞雷波的相速度; ρn为第n层土层的密度；AR=an－1an－2…a1，ai为第i层瑞雷波传递矩阵; 当c>αi时，rαi=
基阶瑞雷波H/V谱的变化形状主要由盆地内各沉积层的厚度和剪切波速度控制，因此，可以通过改变盆地速度模型内各沉积层的厚度或剪切波速度的方式来调整基阶瑞雷波H/V谱(H/V)R0的变化形状，使其与实测地脉动H/V谱(H/V)m的变化形状相同(振幅不同)，以达到修正模型的目的.若通过改变剪切波速度的方法修正模型，只是改变了单个地脉动测点所在位置处地下速度结构中某一个或若干个沉积层的剪切波速度，修改后的速度数据无法利用克里格插值得到新的各沉积层面的展布形态和几何形状；而若通过改变沉积层厚度的方法修正模型，则可改变单个地脉动测点所在位置处地下速度结构中某一个或若干个沉积层分界面的埋深，修改后的分界面埋深数据可直接用于克里格插值，从而得到各沉积层面的展布形态和几何形状，因此本文通过改变各沉积层厚度(即各沉积层分界面埋深)的方法进行模型修正.
玉溪盆地内地脉动测点的分布情况如图 1所示.以图 1中蓝色圆点代表的地脉动测点为例，提出该模型修正方法的依据及具体步骤.在初始建立的速度结构模型中，地脉动测点所在位置的地下速度结构显示，剪切波速度分别为400，600，800，1200，1400，1600 m/s的沉积层分界面的埋深分别为63，116，162，502，663，862 m，基岩面埋深为496 m(何正勤等，2013).剪切波速度为1200 m/s的沉积层分界面的埋深(502 m)已超过基岩面的埋深(496 m)，所以该点地下沉积层速度结构中剪切波速度的最大值为1200 m/s.面波H/V谱的形状对浅层速度结构的变化比较敏感，即浅部速度结构对面波H/V谱的影响更显著(Arai，Tokimatsu，2004)，并且本文只对盆地内沉积层的速度结构进行建模，未考虑基岩面以下地层的速度结构，所以在计算(H/V)R0时，将最深层的沉积层视为无限半空间.在计算地脉动测点所在位置处地下速度结构的(H/V)R0时，将剪切波速度为1200 m/s的沉积层视为无限半空间.
模型修正时只需(H/V)R0与(H/V)m的波峰周期和波谷周期的变化形状相同即可，并不要求振幅相同.基阶瑞雷波H/V谱的公式是依据瑞雷波在水平层状半空间介质中的传播理论推导出来的，使用单点地脉动反演浅层速度结构时，将单个地脉动测点所在位置处的地下速度结构视为水平层状介质，从而反演出单个地脉动测点处的地下速度结构模型.经多次模型修正，提出修正的原则：基阶瑞雷波H/V谱的波峰周期主要由剪切波速度为800，1200和1400 m/s的沉积层分界面的深度所控制，而基阶瑞雷波H/V谱的波谷周期主要由剪切波速度为400和600 m/s的沉积层分界面的深度所控制.若(H/V)R0的波峰周期或波谷周期大于(H/V)m的波峰周期或波谷周期时，各沉积层面的深度要相应减小；反之，各沉积层面的深度则相应增大.一般情况下，玉溪盆地中单个地脉动测点所在位置处的地下速度结构中各沉积层分界面的深度每增加约15 m时，由该点地下速度结构得到的(H/V)R0波峰周期和波谷周期则会增加约0.1 s.
图 3a给出了计算得到的模型修正前地脉动测点(蓝色圆点)的(H/V)R0和(H/V)m。从图中可以看出，初始模型中(H/V)R0的波谷周期比(H/V)m的波谷周期小0.1 s，根据模型修正原则，将剪切波速度为400和600 m/s的沉积层分界面的深度均增加15 m左右，依次修正为75和126 m；(H/V)R0与(H/V)m的波峰周期相差0.2 s，根据模型修正原则，将剪切波速度为800 m/s的沉积层分界面的深度增加30 m左右，修正为192 m.模型修正后的(H/V)R0与(H/V)m如图 3b所示，(H/V)R0与(H/V)m的波峰周期和波谷周期均很一致，证明该地脉动测点所在场地的速度结构模型是合理的.其余各地脉动测点所在场地的速度结构模型也通过该方法进行多次修正，将修正后得到的各地脉动测点的速度结构模型作为新的已知数据，结合钻井资料和地震勘探数据重新进行克里格插值，得到修正后的玉溪盆地速度结构模型.
![]()
图 3 模型修正前(a)、后(b)本文示例中地脉动测点所在场地的地脉动H/V谱(H/V)m(实线)和基阶瑞雷波H/V谱(H/V)R0(虚线)Figure 3. The H/V spectra of microtremors (H/V)m (solid lines) and those of fundamental-mode Rayleigh (H/V)R0 (dashed lines) at the site used for exmaple before (a) and after (b) model updating2.4 模型检验
由于玉溪盆地内缺少强震记录，不能通过对比模拟地震动与实测地震动的波形变化来检验模型的准确度; 而本文通过对比(H/V)R0与(H/V)m的形状变化来检验修正后的玉溪盆地的速度结构模型，则可避免对强震数据的依赖.使用图 1中预留的6个未参与模型修正的黄色圆点(A-F)所代表的地脉动测点，检验修正后的盆地速度模型是否准确。6个地脉动测点所在场地的(H/V)R0和(H/V)m如图 4所示。
![]()
图 4 场地A-F的地脉动H/V谱(H/V)m(实线)和基阶瑞雷波H/V谱(H/V)R0(虚线)Figure 4. The H/V spectra of microtremors (H/V)m (solid lines) and those of fundamental-mode Rayleigh (H/V)R0 (dashed lines) at sites A-F从图 4可以看出，6个场地的(H/V)R0与(H/V)m的形状变化均很相似.场地C和D的(H/V)R0波峰周期比(H/V)m波峰周期小，主要是因为场地C和D附近用于建模的数据均较少，场地C附近的钻孔和地震勘探数据均比较缺乏, 场地D附近仅有的建模数据均为地震勘探得到的剪切波速度，且地震勘探数据本身存在一定误差，所以在初始建立的模型中，场地C和D处存在一些误差，并且玉溪盆地的中东部地区缺少地脉动数据，使得场地D附近缺少修正模型的数据，导致其速度结构并未被修正.由于场地A，B，E和F附近用于初始建模的钻孔和地震勘探数据比较丰富，所以这4个场地的速度结构模型在初始模型中比较准确.此外，场地A，B，E和F附近用于模型修正的地脉动数据也较充足，因此这4个场地中(H/V)R0与(H/V)m的变化形状很相似.除了场地C和D外，场地A，B，E和F的速度模型均能准确地反映真实的地下速度结构，所以使用本文提出的新的建模方法得到的玉溪盆地速度结构模型的准确度较高，证明该建模方法是可靠的.
图 5给出了修正后的玉溪盆地内沉积层的三维速度结构模型.图中沉积层分界面的剪切波速度由上到下依次为400，600，800，1200，1400 m/s，基岩面与剪切波速度为1400 m/s的沉积层分界面之间的沉积层的剪切波速度为1600 m/s.
![]()
图 5 玉溪盆地三维速度结构. x轴为经度方向，y轴为纬度方向，z轴为深度方向Figure 5. The 3-D velocity structures of the Yuxi basinx-axis is the longitude direction, y-axis is the latitude direction, and z-axis is the depth direction玉溪盆地是发育于普渡河断裂带南端的沉积盆地，沉积层厚度变化较大，盆地边缘沉积层的厚度仅约100 m，而盆地西北部厚度可达800 m，根据岩性差异可分为5—6层(何正勤等，2013).这与本文建立的玉溪盆地速度模型(图 5)相符：该模型中基岩面埋深的变化较大，盆地边缘地区基岩面埋深较浅，由于盆地中部地区太极山的出露使得该地区的基岩埋深较浅；该模型在沉积层厚度较小的区域有4—5个沉积层，在沉积层厚度较大的中部区域有6个沉积层，与玉溪盆地有5—6层不同岩性沉积层的地质情况相符合.
3. 讨论与结论
沉积层对地震动的振幅和持时有显著的放大作用，增大了位于盆地内的城市震害.盆地内的地震动模拟需要建立精确的沉积盆地三维速度结构模型.本文基于已有的数据，提出了一种包括数据预处理、模型建立、模型修正和模型检验的沉积盆地三维速度结构的建模方法.在数据预处理时，依据各类数据间不同的可信度，拟合得到消除速度偏差的折减函数, 然后根据(H/V)R0与(H/V)m的变化形状相似的原则对模型进行修正和检验.该模型修正方法避免了对地震波形数据的依赖, 使得该模型修正方法具有广泛的适用性.使用该模型修正方法得到修正后的玉溪盆地三维速度结构模型，经检验准确度较高，说明了该模型修正方法的合理性.
本文通过地脉动H/V谱的形状变化对比的方法对模型进行修正，未考虑地脉动H/V谱振幅变化的影响，造成修正后的模型还存在一些误差.同时在模型修正时，存在主观因素的干扰，如何避免此类干扰是一个亟需解决的问题.另外，地脉动数据的缺少也极大地限制了模型修正的精度，例如玉溪盆地的中东部地区缺少地脉动数据，造成该地区模型无法进行修正，因此该地区模型的准确度较低.
-
![]()
图 1 1999年9月21日台湾集集MW7.6地震断层和台站分布示意图
Figure 1. Schematic diagram of the fault and stations distribution of the MW7.6 Chi-Chi earthquake in Taiwan on 21 September 1999
![]()
图 2 速度时程基线校正示意图
Figure 2. Schematic diagram of baseline correction for velocity time history
![]()
图 3 本文方法与太平洋地震工程研究中心方法处理的TCU068台站近断层地震动位移时程比较
（a） NS分量；（b） EW 分量；（c） UD 分量
Figure 3. Comparison of near fault ground motion displacement time histories of the station TCU068 processed by the method proposed in this study with those from Pacific Earthquake Engineering Research Center
（a） NS component；（b） EW component；（c） UD component
![]()
图 4 上盘台站强震动记录三分量EW （左），NS （中），UD （右）分量的位移时程
蓝线表示位移的最值，红线表示位移初始值或最终值，下图同
Figure 4. Displacement time histories of EW （left），NS （middle） and UD （right） components of strong ground motion records at the hanging wall stations
The blue line indicates the extreme displacement，and the red line indicates the initial or final displacement，the same below （c） TCU072；（d） TCU071；（e） TCU089；（f） TCU078
![]()
图 5 强震动位移时程形态典型特征
（a） TCU052台站EW分量；（b） TCU054台站NS分量；（c） TCU075台站UD分量
Figure 5. Typical characteristics of strong ground motion displacement time histories
（a） EW component at the station TCU052；（b） NS component at the station TCU054；（c） UD component at the station TCU075
![]()
图 6 六对上、下盘台站的EW （左），NS （中）和UD （右）分量的相对位移时程
Figure 6. The relative displacement time history of EW （left），NS （middle），and UD （right） components of the six pairs of the hanging wall and footwall stations
（e） TCU089-TCU076；（f） TCU078-TCU122
![]()
图 7 上下盘台站最大相对位移与相对残余位移关系
Figure 7. Relationship between maximum relativedisplacement and relative residual displacement of hanging wall and footwall stations
![]()
图 8 上、下盘台站强震动相对位移时程典型特征
（a） TCU052-TCU054台站的EW分量；（b） TCU072-TCU067台站的EW分量；（c） TCU072-TCU067台站的UD分量
Figure 8. Typical characteristics of relative displacement time history of strong ground motion records at hanging wall and footwall stations
（a） EW component of the station TCU052-TCU054；（b） EW component of the station TCU072-TCU067；（c） UD component of the station TCU072-TCU067
![]()
图 9 集集地震土木工程断层破坏7例
（a） 石岗大坝（TCU068-TCU103）；（b） 北丰桥（TCU068-TCU103）；（c） 欢乐谷公寓群（TCU068-TCU103）；（d） 第一公园大桥（TCU052-TCU054）；（e） 光复中学（TUC071-TCU075）；（f） 乌溪桥（TUC071-TCU075）；（g） 刘眉桥（TCU078-TCU122）
Figure 9. Seven cases of civil engineering fault failure in Chi-Chi earthquake
（a） Shigang Dam （TCU068-TCU103）；（b） Beifeng Bridge （TCU068-TCU103）；（c） Happy Valley Apartments （TCU068-TCU103）；（d） The First Park Bridge （TCU052-TCU054）；（e） Guangfu Middle School （TUC071-TCU075）；（f） Wuxi Bridge （TUC071-TCU075）；（g） Liumei Bridge （TCU078-TCU122）
![]()
图 10 四对上、下盘台站的空间相对位移时程
Figure 10. Spatial relative displacement time histories of the four pairs of hanging wall and footwall stations
（a） TCU068-TCU103；（b） TCU052-TCU054；（c） TCU071-TCU075；（d） TCU078-TCU122
![]()
图 11 跨断层场点间相对位移与场点间距的关系
Figure 11. Relationship between relative displacement and field spacing across fault sites
![]()
图 12 逆断层盘间相对位移趋势示意图
Figure 12. Schematic diagram of relative displacement trend between reverse fault walls
表 1 1999年台湾集集MW7.6地震强震动台站位置
Table 1 Location of ground motion stations of 1999 MW7.6 Chi-Chi earthquake in Taiwan region
断层盘 台站 东经/° 北纬/° 海拔/km 断层距/km 上盘 TCU068 120.8 24.3 0.276 0.071 TCU052 120.7 24.2 0.170 0.092 TCU072 120.8 24.0 0.363 13.021 TCU071 120.8 24.0 0.187 10.003 TCU089 120.9 23.9 0.020 15.277 TCU078 120.8 23.8 0.272 14.402 下盘 TCU103 120.7 24.3 0.222 5.982 TCU054 120.7 24.2 0.097 7.029 TCU067 120.7 24.1 0.073 0.716 TCU075 120.7 24.0 0.096 1.209 TCU120 120.6 24.0 0.228 7.793 TCU118 120.4 24.0 0.008 30.496 TCU076 120.7 24.0 0.103 3.653 TCU129 120.7 23.9 0.110 3.059 TCU122 120.6 23.8 0.075 10.996 CHY024 120.6 23.8 0.085 12.255 注：数据来源于台湾气象厅公开发行的光盘（ Lee et al，2001 ）。
下载: 导出CSV
表 2 1999年集集MW7.6地震上、下盘台站本文所选组对情况
Table 2 Group matching of hanging wall and footwall stations of the 1999 MW7.6 Chi-Chi earthquake in this paper
台站 组队 上盘 TCU068 TCU052 TCU072 TCU071 TCU089 TCU078 下盘 TCU103 TCU054 TCU067 TCU075 TCU075 TCU122
下载: 导出CSV
表 3 上盘台站强震动记录位移最大值及残余位移
Table 3 Maximum displacement and residual displacement of strong ground motion records at the hanging wall stations
上盘台站 EW分量 NS分量 UD分量 最大位移/m 残余位移/m 最大位移/m 残余位移/m 最大位移/m 残余位移/m TCU068 7.06 5.81 8.61 5.94 4.52 3.45 TCU052 4.99 3.94 7.18 6.78 3.92 3.16 TCU072 2.02 1.80 2.35 2.00 1.33 1.16 TCU071 1.76 1.19 2.61 0.15 2.26 1.74 TCU089 1.91 1.81 1.37 1.18 0.24 0.03 TCU078 1.21 1.19 0.85 0.66 0.35 0.25
下载: 导出CSV
表 4 六对上、下盘台站强震动相对位移最大值及相对残余位移
Table 4 The maximum relative displacement and relative residual displacement of strong ground motion records at six pairs of hanging wall and footwall stations
上、下盘台站 EW分量 NS分量 UD分量 最大相对位移/m 相对残余位移/m 最大相对位移/m 相对残余位移/m 最大相对位移/m 相对残余位移/m TCU068-TCU103 7.54 6.21 9.11 6.47 4.09 3.33 TCU052-TCU054 5.57 4.50 7.81 7.08 3.81 3.29 TCU072-TCU067 3.46 3.02 3.27 2.97 1.59 1.41 TCU071-TCU075 2.89 2.28 2.99 2.52 2.55 1.74 TCU089-TCU076 2.82 2.50 1.72 1.42 0.39 0.26 TCU078-TCU122 2.06 1.82 1.26 1.16 0.35 0.18
下载: 导出CSV
表 5 工程破坏现场地表位错与附近断层上、下盘台站位错计算值的比较
Table 5 Comparison of calculated values of surface dislocations at the site of engineering failure with those at the hanging wall and footwall stations nearby faults
工程名称 上、下盘台站 台站间距/km 位错观察值/m 位错计算值/m 数据记录方向 石岗大坝 TCU068-TCU103 6.05 8.00 3.33 上下 北丰桥 5.50 上下 欢乐谷公寓群 3.50 上下 第一公园大桥 TCU052-TCU054 7.12 4.00 3.29 上下 光复中学 TCU071-TCU075 11.21 2.50 1.74 上下 乌溪桥 2.00 2.61 水平 刘眉桥 TCU078-TCU122 25.40 0.20 0.18 上下
下载: 导出CSV
-
曹志磊,周琼,葛计划,孙军,隆爱军,龙剑锋,赵希磊. 2019. 郯庐断裂带安徽段及“霍山窗”断层活动特征与地震关联性研究[J]. 大地测量与地球动力学,39(7):681–685. Cao Z L,Zhou Q,Ge J H,Sun J,Long A J,Long J F,Zhao X L. 2019. Fault activity and correlation study of Tan-Lu fault zone and “Huoshan Seismic Window"[J]. Journal of Geodesy and Geodynamics,39(7):681–685 (in Chinese).
陈勇,陈鲲,俞言祥. 2007. 用集集主震记录研究近断层强震记录的基线校正方法[J]. 地震工程与工程振动,27(4):1–7. doi: 10.3969/j.issn.1000-1301.2007.04.001 Chen Y,Chen K,Yu Y X. 2007. Base line correction method for near-fault accelerograms using Chi-Chi main shock record[J]. Journal of Earthquake Engineering and Engineering Vibration,27(4):1–7 (in Chinese).
高波,王峥峥,袁松,申玉生. 2009. 汶川地震公路隧道震害启示[J]. 西南交通大学学报,44(3):336–341. Gao B,Wang Z Z,Yuan S,Shen Y S. 2009. Lessons learnt from damage of highway tunnels in Wenchuan earthquake[J]. Journal of Southwest Jiaotong University,44(3):336–341 (in Chinese).
胡进军. 2009. 近断层地震动方向性效应及超剪切破裂研究[D]. 哈尔滨: 中国地震局工程力学研究所: 126–133. Hu J J. 2009. Directivity Effect of Near-Fault Ground Motion and Super-Shear Rupture[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 126–133 （in Chinese）.
胡聿贤. 2006. 地震工程学[M]. 第二版. 北京: 地震出版社: 13–14. Hu Y X. 2006. Earthquake engineering[M]. Second edition. Beijing: Seismological Press: 13–14 （in Chinese）.
黄润秋,李为乐. 2009. 汶川大地震触发地质灾害的断层效应分析[J]. 工程地质学报,17(1):19–28. Huang R Q,Li W L. 2009. Fault effect analysis of geo-hazard triggered by Wenchuan earthquake[J]. Journal of Engineering Geology,17(1):19–28 (in Chinese).
李爽,周洪圆,刘向阳,贾俊峰. 2020. 基于中国规范的近断层区抗震设计谱研究[J]. 建筑结构学报,41(2):7–12. Li S,Zhou H Y,Liu X Y,Jia J F. 2020. Study on near-fault seismic design spectra based on Chinese Code for Seismic Design of Buildings[J]. Journal of Building Structures,41(2):7–12 (in Chinese).
王栋,谢礼立,胡进军. 2008. 倾斜断层不对称分布引起的几何效应:上下盘效应[J]. 地震学报,30(3):271–278. doi: 10.3321/j.issn:0253-3782.2008.03.006 Wang D,Xie L L,Hu J J. 2008. Geometric effects resulting from the asymmetry of dipping fault:Hanging wall/footwall effects[J]. Acta Seismologica Sinica,30(3):271–278 (in Chinese).
王永安,李琼,刘强. 2011. 跨断层形变累积率的变化特征与云南地区强震[J]. 地震研究,34(2):136–142. Wang Y A,Li Q,Liu Q. 2011. Variation of the cumulative rate of cross-fault deformation and strong earthquakes in Yunnan[J]. Journal of Seismological Research,34(2):136–142 (in Chinese).
谢礼立,徐龙军,陶晓燕,杨绪剑. 2021. 跨断层土木工程研究与实验装置研发现状[J]. 工程力学,38(4):20–29. Xie L L,Xu L J,Tao X Y,Yang X J. 2021. Research status of civil engineering structures across faults and the development of experimental devices for fault simulation[J]. Engineering Mechanics,38(4):20–29 (in Chinese).
喻畑,李小军. 2012. 基于NGA模型的汶川地震区地震动衰减关系[J]. 岩土工程学报,34(3):552–558. Yu T,Li X J. 2012. Attenuation relationship of ground motion for Wenchuan earthquake region based on NGA model[J]. Chinese Journal of Geotechnical Engineering,34(3):552–558 (in Chinese).
俞言祥,高孟潭. 2001. 台湾集集地震近场地震动的上盘效应[J]. 地震学报,23(6):615–621. Yu Y X,Gao M T. 2001. Effects of the hanging wall and footwall on peak acceleration during the Jiji (Chi-Chi),Taiwan Province,earthquake[J]. Acta Seismologica Sinica,14(6):654–659.
张红艳,谢富仁. 2013. 天山地区跨断层形变观测与地壳应力特征[J]. 西北大学学报(自然科学版),43(4):617–622. Zhang H Y,Xie F R. 2013. Cross-fault deformation observation and characteristics of crustal stress in Tianshan region[J]. Journal of Northwest University (Natural Science Edition),43(4):617–622 (in Chinese).
周云好,陈章立,缪发军. 2004. 2001年11月14日昆仑山口西MS8.1地震震源破裂过程研究[J]. 地震学报,26(增刊):9–20. Zhou Y H,Chen Z L,Miao F J. 2004. Source process of the 14 November 2001westhern Kunlun mountain MS8.1 earthquake[J]. Acta Seismological Sinica,26(S1):9–20 (in Chinese).
Abrahamson N A. 2000. Near-fault ground motions from the 1999 Chi-Chi earthquake[C]//Proc. of US-Japan Workshop on the Effects of Near-Field Earthquake Shaking. San Francisco, California: Pacific Earthquake Engineering Research Center: 11−13.
Bolt B A. 1999. Earthquakes[M]. Fourth edition. New York: W.H. Freeman and Company: 52–53.
Boore D M. 2001. Effect of baseline corrections on displacements and response spectra for several recordings of the 1999 Chi-Chi,Taiwan,earthquake[J]. Bull Seismol Soc Am,91(5):1199–1211.
Bray J D,Rodriguez-Marek A. 2004. Characterization of forward-directivity ground motions in the near-fault region[J]. Soil Dyn Earthq Eng,24(11):815–828. doi: 10.1016/j.soildyn.2004.05.001
Chao W A,Wu Y M,Zhao L. 2009. An automatic scheme for baseline correction of strong-motion records in coseismic deformation determination[J]. J Seismol,14(3):495–504.
Chu D B, Brandenberg S J, Lin P S. 2008. Performance of bridges in liquefied ground during 1999 Chi-Chi earthquake[C]//The 14th WCEE. Beijing: CAEE and IAEE: 13–17.
Dong J J, Wang C D, Lee C T, Liao J J, Pan Y W. 2004. The influence of surface ruptures on building damage in the 1999 Chi-Chi earthquake: A case study in Fengyuan City[J]. Engineering Geology, 71（1/2）: 157–179.
Iwan W D,Moser M A,Peng C Y. 1985. Some observations on strong motion earthquake measurement using a digital accelerograph[J]. Bull Seismol Soc Am,75(5):1225–1246. doi: 10.1785/BSSA0750051225
Kawashima K. 2002. Damage of bridges resulting from fault rupture in the 1999 Kocaeli and Duzce, Turkey earthquakes and the 1999 Chi-Chi, Taiwan earthquake[J]. Structural Eng/Earthq Eng, 19（2 Special Issue）: 179−197.
Lee W H K, Shin T C, Kuo K W, Chen K C, Wu C F. 2001. CWB free-field strong-motion data from the 21 September Chi-Chi, Taiwan, earthquake[J]. Bull Seismol Soc Am, 91（5）: 1370−1376.
Lin M L,Lin C H,Li C H,Liu C Y,Hung C H. 2021. 3D modeling of the ground deformation along the fault rupture and its impact on engineering structures:Insights from the 1999 Chi-Chi earthquake,Shigang District,Taiwan[J]. Eng Geol,281:105993. doi: 10.1016/j.enggeo.2021.105993
Ma K F,Mori J,Lee S J,Yu S B. 2001. Spatial and temporal distribution of slip for the 1999 Chi-Chi,Taiwan earthquake[J]. Bull Seismol Soc Am,91(5):1069–1087.
Mccomb H,Ruge A,Neumann F. 1943. The determination of true ground motion by integration of strong-motion records:A symposium[J]. Bull Seismol Soc Am,33(1):1. doi: 10.1785/BSSA0330010001
Ota Y,Watanabe M,Suzuki Y,Yanagida M,Miyawaki A,Sawa H. 2007. Style of the surface deformation by the 1999 Chichi earthquake at the central segment of Chelungpu fault,Taiwan,with special reference to the presence of the main and subsidiary faults and their progressive deformation in the Tsauton area[J]. J Asia Earth Sci,31(3):214–225. doi: 10.1016/j.jseaes.2006.07.030
Somerville P G,Smith N F,Graves R W,Abrahamson N A. 1997. Modification of empirical strong ground motion attenuation relations to include the amplitude and duration effects of rupture directivity[J]. Seismol Res Lett,68(1):199–222. doi: 10.1785/gssrl.68.1.199
Wang G Q,Zhou X Y,Zhang P Z,Igel H. 2002. Characteristics of amplitude and duration for near fault strong ground motion from the 1999 Chi-Chi,Taiwan earthquake[J]. Soil Dyn Earthq Eng,22(1):73–96. doi: 10.1016/S0267-7261(01)00047-1
Wang R J,Schurr B,Milkereit C,Shao Z G,Jin M P. 2011. An improved automatic scheme for empirical baseline correction of digital strong-motion records[J]. Bull Seismol Soc Am,101(5):2029–2044. doi: 10.1785/0120110039
Wang W L,Wang T T,Su J J,Lin C H,Seng C R,Huang T H. 2001. Assessment of damage in mountain tunnels due to the Taiwan Chi-Chi earthquake[J]. Tunn Undergr Space Technol,16(3):133–150. doi: 10.1016/S0886-7798(01)00047-5
Wu Y M, Wu C F. 2007. Approximate recovery of coseismic deformation from Taiwan strong-motion records[J]. J Seismol, 11（2）: 159–170.
Xu L J,Zhao G C,Gardoni P,Xie L L. 2018. Quantitatively determining the high-pass filter cutoff period of ground motions[J]. Bull Seismol Soc Am,108(2):857–865. doi: 10.1785/0120170074
Zhao G C,Xu L J,Gardoni P,Xie L L. 2019. A new method of deriving the acceleration and displacement design spectra of pulse-like ground motions based on the wavelet multi-resolution analysis[J]. Soil Dyn Earthq Eng,119:1–10. doi: 10.1016/j.soildyn.2019.01.008
-
期刊类型引用(2)
1. 师黎静，宋健，党鹏飞，刘佳轩. 区域场地近地表速度结构建模研究. 岩土工程学报. 2022(02): 360-367 . 
百度学术
2. 万远春，于彦彦，丁海平，胡颖平. 考虑不均匀地壳构造的四川盆地地震动模拟研究. 自然灾害学报. 2022(02): 204-214 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 396
- HTML全文浏览量: 98
- PDF下载量: 103
- 被引次数: 5
