2011年新西兰Mw6.1地震震源过程及强地面运动特征初步分析

孟令媛，史保平

1 中国科学院研究生院，北京 100049

2 中国地震局台网中心，北京 100036

2011年新西兰Mw6.1地震震源过程及强地面运动特征初步分析

孟令媛1，2，史保平1

1 中国科学院研究生院，北京 100049

2 中国地震局台网中心，北京 100036

2011年2月21日新西兰南岛发生Mw6.1地震，震中位置43.58°S，172.70°E，震源深度约5.0km，发震断层为新西兰第三大城市克赖斯特彻奇（Christchurch）南约9km一条近东西走向的未知隐伏断层，为逆冲断层机制.地震已经造成163人遇难，Christchurch城内多处建筑物严重损毁，距震中约1km的Heathcote Valley Primary School（HVSC）台站强地面运动峰值加速度（PGA）高达2.0g.针对新西兰Mw6.1地震近场强地面运动偏高这一现象，利用HVSC台站的强震观测记录，计算地震震源谱参数，应用Brune圆盘模型估算其近场强地面运动的理论值，并建立动态复合震源模型进行模拟计算.研究结果表明，新西兰Mw6.1地震近断层强地面运动偏高的主要原因为复杂震源破裂过程中有效应力降（动态应力降，Δσd）过高造成的.未来工作中，需要加强对可能发生的、距离城市较近的中小型地震的重视，防止地震对城市的加强型破坏.

新西兰Mw6.1地震，强震观测记录，有效应力降，强地面运动，动态复合震源模型

1 引 言

2011年2月21日23时51分43秒新西兰南岛Greendale的SE方向约30km处附近发生Mw6.1地震（United State Geology Survey，USGS，2011），震中位置43.58°S，172.70°E，距离新西兰第三大城市克赖斯特彻奇（又名基督城，Christchurch）约6km，震源深度约5.0km，发震断层为Christchurch南约9km处一条近东西走向的隐伏断层，为逆冲断层机制，其后发生4级以上余震23次，5级以上余震4次，最大余震5.7级.截至2011年3月4日，地震已经造成163人遇难，最终死亡人数预计将超过200人，地震发生后，Christchurch市中心多处建筑物受损，城市阶段路面多处出现裂缝，地标性建筑Christchurch大教堂也部分被毁（GNS Science，Te PüAo，New Zealand，2011）［1］.研究结果认定新西兰Mw6.1地震为2010年新西兰Mw7.0地震的余震，且该余震是迄今为止新西兰Mw7.0地震发生后的最大的一次余震（见图1）.2010年9月3日16时35分46秒新西兰南岛Greendale附近发生Mw7.0地震（USGS，2010），震中位置43.53°S，172.12°E，震源深度约10.0km，致使2人受伤，并造成Christchurch城内多处建筑物损毁［2－3］.

图1 2011年新西兰Mw6.1地震震中及发震断层地表出露位置Fig.1 Epicenter location of the 2011 Mw6.1，New Zealand Earthquake，and the distribution of the blinded fault

2010年新西兰Mw7.0地震的零死亡及建筑物低损失一度引起多方对新西兰防震减灾技术的关注.随着2011年新西兰Mw6.1地震的发生，并造成了较之Mw7.0主震巨大的人员伤亡和财产损失，又一次引起了国内外学者对新西兰Mw6.1地震的广泛关注和讨论.新西兰位于太平洋地震带上太平洋板块（Pacific plate）与澳大利亚板块（Australia plate）的交界处，受两板块的相互作用，地震活动十分强烈［2－4］.自1840年以来，新西兰共发生22次7级以上大震，1855年1月23日的惠灵顿大地震（Mw8.2）是新西兰迄今已知的最大地震，2011年新西兰Mw6.1地震共造成死亡人最终将超过200人，是新西兰历史上地震灾害最重的一次地震（中国地震信息网，2011）［5］.

2011年2月21日新西兰Mw6.1地震发生后，新西兰本土多处强震台站记录到了强震数据.截止到2011年5月Center for Engineering Strong Motion Data（CESMD）公布的数据中，记录到的强地面运动加速度的峰值（Peak Ground Acceleration，PGA）达0.001g及以上的台站数据共有106个，其中Heathcote Valley Primary School（HVSC）台站距离震中最近，仅约为1.0km（见图1），该台站记录到的PGA峰值高达2.0g（垂向分量），为新西兰Mw6.1地震发震后该项记录的最高值，此外，新西兰Mw6.1地震记录到的PGA达0.5g及以上的台站数据共5个.尽管目前记录到的高的强地面运动峰值，可以对该地震带来的巨大破坏和引发的人员伤亡给出一定的解释，但对2011年新西兰Mw6.1地震震源过程进一步的探讨，可以加深对近城市的中小型地震潜在破坏力的认识［6］.

我们针对2011年新西兰Mw6.1地震的关注和研究，一定程度上源自于对2010年新西兰Mw7.0地震研究的延续，前文中已提到，2010年新西兰Mw7.0地震即目前认为的2011年新西兰Mw6.1地震的主震.尽管2010年新西兰Mw7.0主震并没有造成过大的人员和财产的损失，但值得关注的是，该地震所造成的强地面运动也很大（CESMD，2010），最开始公布的数据最高值约为0.77g，该记录已略高于新一代的衰减关系（Next Generation Attenuation Relation，NGA，2007），2011年新西兰Mw6.1地震强地面运动的观测值则更加偏高［6－8］，这一点也是我们后面会探讨和讨论的问题之一，即为什么一个矩震级6级的隐伏型地震能够造成强地面运动的如此高值？

本文中我们着重讨论新西兰Mw6.1地震的震源过程，尝试性地解释该地震造成强地面运动峰值偏高、进而引发高破坏力可能的原因.工作主要分为以下几个部分进行：首先，利用HVSC台站的强震观测数据，计算该地震的震源谱等参数，如速度平方谱、位移平方谱等，估算该地震的震源参数，并结合已有的反演结果，进行对比和比较；第二步，由于该地震的发震断层为隐伏型断层，利用经验关系回归该地震的断层面积，应用已有震源参数结合理论模型，计算该地震强地面运动的理论值；第三步，引入新一代衰减关系（NGA）对观测数据进行整体分析，建立有限断层运动学模型模拟该台站的强地面运动，并与观测记录进行对比，探讨该地震发震后近场强地面运动高值的可能成因；最后进行讨论并给出结论.

2 震源参数及强地面运动特征

2.1 利用观测数据计算震源谱参数

前文中已经提到，该地震发生在新西兰南岛Greendale的SE方向约30km处附近，距离Christchurch约6km，发震断层为Christchurch南约9km处一条近东西走向的逆冲、隐伏型断层（见图1），新西兰Mw6.1地震的震源参数详见表1（GNS，2011；USGS，2011）［1－2］.

表1 2011年新西兰Mw6.1地震震源参数Table 1 Source parameters of the 2011 Mw6.1 New Zealand Earthquake

本文在研究过程中首先从强震观测数据出发，选取HVSC台站的强震观测记录作为参照（见图2）.图2中分别给了HVSC台站加速度三分量的观测记录，由图2中可以看出，由于该台站距离震中很近，约为1km，且发震断层为隐伏型断层，因此该台站距发震断层的距离尚不明确.其中垂向（Vertical）、S26W向和S64E向三个分量加速度的峰值分别为2160.5cm·s－2、1646.8cm·s－2及1244.6cm·s－2，速度的峰值分别为97.3cm·s－1、62.3cm·s－1及39.2cm·s－1（CESMD，2011）［6］.由图2（a、b、c）可以看出，HVSC台站的强震观测记录中较为明显地呈现了高频成分相对集中的特征，主要集中在10～20s这一时间段.也就是说，参照HVSC台站计算新西兰Mw6.1地震的地震波辐射能，可通过计算10～20s这一时间区域内的地震波辐射能来近似获得.

简单地说，地震矩率函数在时间域中表示为标量地震矩乘以震源时间函数，在频率域中则表示为震源谱［9－11］：

式（1）中M0为地震矩，f为频率，f0为拐角频率（corner frequency），γ为高频衰减常数.通常把振幅谱高频渐近趋势（包络线）和低频趋势（零频水平）的交点称作拐角（corner），与拐角相应的频率称作拐角频率（corner frequency）.一般情况下，f0很难在频谱图上直接准确测得，这里我们采用Brune模型［9－10］（1970，1971），即令γ＝2.当震源距（r）很小，非弹性衰减的影响很小且可以忽略时，Andrews提出［12］可以通过积分的方式求解震源谱的参数，如速度平方谱及位移平方谱等，即对于近震观测的情况，强地面运动位移的频谱可以表示为

其中D（f）为位移的频谱，D2（f）为位移的平方谱（D2（f）＝V2（f）·（2πf）－2，V2（f）为速度的平方谱），f为频率，则SD2为位移平方谱的积分.同理有

且有拐角频率f0可通过公式（2）和（3）的积分求得，即

Andrews指出［12］上式中f0的计算方式可以应用于任何频谱，因此可将其作为拐角频率的定义式，Andrews同时提出，对于点源或圆盘断层模型，其地震波辐射能可以有下式近似求得：

上式中ρ为密度，β为S波的波速.

图2 2011年新西兰Mw6.1地震HVSC台站强震观测记录，（a，b，c）中上为0～50s的观测记录，下为10～20s的观测记录.（a）垂向分量，观测记录PGA为2160.5cm·s－2；（b）S26W分量，观测记录PGA为1646.8cm·s－2；（c）S64E分量，观测记录PGA为1244.6cm·s－2Fig.2 Observed accelerations of the three－component time－histories at the station HVSC，New Zealand，Mw6.1 Earthquake.（a）Vertical component，the PGA is 2160.5cm·s－2；（b）S26Wcomponent，the PGA is 1646.8cm·s－2；（c）S64Ecomponent，the PGA is 1244.6cm·s－2

在研究过程中，我们首先从HVSC台站的观测数据出发，截取加速度时程曲线中高频成分集中的10～20s一段（详见图2），将该段持续时间为10s的区域作为我们计算区域；第二步，通过积分分别求得10～20s时域内的速度平方谱和位移平方谱，即应用式（2）和式（3）分别计算位移平方谱的积分和速度平方谱的积分；最后应用式（4）和式（5）分别计算拐角频率f0和地震波辐射能Erad，得到f0约为0.22Hz，Erad约为1.2×1014J.结合表1中给出的远场地震波辐射能（Es）的反演结果，Es＝（2.4～3.6）×1014J.考虑到我们利用HVSC台站强震观测记录进行计算过程中，仅截取了10s持续时间范围内的地震波辐射能，计算过程中仅考虑S波的地震波辐射能.HVSC台站的震源距约为1km，该台站受到地震波辐射影响的持续时间也就相应较短，因此通过该台站实际的观测记录所计算得到的Erad，约为远场反演结果Es的35%～46%，仍可视为合理.对于f0的计算值，后面在估算强地面运动的理论值中将会得到应用.

我们尝试从远场模型及相关方法入手，对新西兰Mw6.1地震HVSC台站记录到的强震观测数据，尝试性地用频谱的相关理论和方法进行分析，即Brune模型的远场RMS平均谱公式（式（1））及相应Andrews的方法（式（5））［9，12］.对于选择台站及观测数据，理想情况是选择尽量相对远场的台站及数据，根据表1中给出的断层尺度（15km×9km），最好选择也要在30km以外的台站，仔细分析2011年新西兰Mw6.1地震PGA观测值较大的前20个台站，可知均位处震源≤30km的范围内.因此，我们又仔细对现有的震源反演数据进行了深入的分析，详见图3a（下图），断层面上主要的Asperity（滑动集中区）仅有约4km×4km在≥1.5m的范围内，≥2.0m范围内的Asperity则仅有约2km×2km（见图3a中红色圆圈标示）.

图3a为GNS给出的2011年新西兰Mw6.1地震断层面上的反演结果［1］，红色圆圈所标示的位置显示了滑动位移的集中区，从而表明地表运动的主要影响来自于这一高滑移区.实际上我们做了一定的近似，近似地认为新西兰Mw6.1地震断层面上起到主要影响的Asperity尺度半径约为2～4km.对于Brune近场圆盘模型，大致可近似认为是远场的范围就相应扩大，从另一个角度讲，近场和远场均是相对而言的，尽管HVSC台站距离震中约1.0km，但毕竟2011年新西兰Mw6.1地震的发震断层是未出露地表的隐伏型断层［1，2，5］，图3a显示断层顶端大致距离地表1.0km（这个深度可能更大），震中位置在发震断层面上沿倾向（dip）约7.0km.因此，我们大胆假设HVSC台站可能恰好在远场和近场的分界处，加之HVSC台站的观测（未滤波的垂向分量）PGA值最高，这样一来，选择HVSC台站的观测数据作为分析的重点，则更加有意义.

图3 （a）2011年新西兰Mw6.1地震发震断层面上反演结果（GNS，2011）；（b）应用本文建立的动态复合震源模型给出的断层面上滑动位移的不均匀分布特征；（c）动态复合震源模型上走向、倾向及滑动角等参数平面示意图；（d）动态复合震源模型参数空间分布图Fig.3 （a）The slip distributions on the main fault generated for the 2011 Mw6.1New Zealand Earthquake（GNS，2011）；（b）Inhomogeneous distribution of slip on the main fault in our dynamical composite source model；（c）Simplified spatial distribution of slip－displacement on the fault plan used（Dip＝68°，Strike＝58°，Rake＝45°，Rake＝45°，Rake＝135°，Rake＝0°）in our model；（d）Spatial distribution of slip－displacement on the fault plan same as（c）

具体意义在于以下几点：其一，前文中我们在强调成文意义时已提出为何一个矩震级6级的地震能够造成如此高的PGA；其二，对于Brune圆盘模型的假设，其中最重要即是远场角频率（f0）概念的应用，在应用DCSM（Dynamical Composite Source）模型的过程中，对于近场强地面运动的合成，也是基于远场f0的假设，即当我们将断层面上的大个子源认为是小个子源的合成结果，则可以进一步将小子源认为是符合远场理论的点源，只有承认了这个假设，其后的模拟工作才有意义；最后，需要说明的是，在DCSM模型中，图3a中标示的红色圆圈这一高的滑移区可有几个不同大小的Brune Model合成，而地表运动则是不同大小的Brune Pulse线性叠加的结果，如果采用了相同的Green′s Function，在理论分析中为了简化问题，我们往往近似用一个Circular Model进行理论值的估算.总体而言，文章中选择用HVSC台站的分析和应用，在一定的近似简化下可认为不存在绝对的近场和远场的矛盾.

2.2 估算强地面运动的理论值

进一步估算新西兰Mw6.1地震强地面运动的理论值，由于该地震的发震断层为隐伏型断层，且未破裂至地表，需要对其断层的破裂面积进行一定的验算.应用Hanks和Bakun［13］2002年给出的地震矩震级与断层破裂面积之间的双线性经验关系，对该地震的断层破裂面积进行初步的计算：

Mw＝lgA＋3.98±0.03，A≤537km2，（6）上式中Mw为地震矩震级，A为断层的破裂面积，式（6）在Hanks和Bakun［13］文中原为分段函数，即以A取537km2为分段点，当A＝537km2时，Mw≈6.71，由于本文涉及的2011年新西兰地震矩震级小于6.71，因此我们在验算该地震的破裂面积时，仅采用式（6）中给出的这一部分.

参照表1和式（6）可以推算出，新西兰Mw6.1地震发震断层的断层面积约为123～141km2，这一结果与GNS最新的反演结果相一致（详见表1，图3a）.新西兰Mw6.1地震发震断层面积相对较小，研究过程中将其发震断层近似为圆盘模型进行理论计算，首先对圆盘模型的半径（R）进行初步的估算［9，10］：

上式中β为S波的波速，且β＝3×105cm·s－1，f0为拐角频率，前文中给出的计算结果为0.22Hz，则有R≈5.0km，令圆盘模型的面积为Ac，则有Ac＝πR2≈80km2，Ac的估算结果与图3a中滑动集中区的有效破裂面积相一致，因此可视为理论估算模型合理.

Andrews［12］对地震强地面运动的频谱响应进行分析时指出：

其中σB为Brune圆盘模型下的有效应力降，即动态应力降Δσd，AFS为S波辐射图型因子，其在震源球球面上的均方根为AFS＝，σa为视应力［14］（σa＝μEs／M0，其中μ为剪切模量，Es为地震波辐射能，M0为地震矩）.参照表1，新西兰Mw6.1地震的σa约为4.5MPa，则Δσd约为19.35MPa，进一步采取Brune［15］圆盘模型下给出的近场强震质点运动速度的估算方法，即

Brune圆盘模型下的理论估算结果有些偏高（参照图2），且由式（9）可以发现，采用Brune圆盘模型对质点加速度和质点速度进行估算，主要决定于Δσd的大小.因此，通过理论模型的估算，我们可以初步给出这样的假设，对于2011年新西兰Mw6.1地震，尽管其矩震级尺度为6级，但其造成的近场较高的强地面运动（GNS，2011；CESMD，2011）［1，6］有可能来自于震源破裂过程中动态应力降（Δσd）偏高的影响.关于Δσd取值的高低是否影响到近场强地面运动的强弱，我们在下面将通过模型计算并对比观测结果来进行详细的讨论.

2.3 浅层速度结构（VS30）的影响

2011年新西兰Mw6.1地震发生后，CESMD目前公布的数据中，记录到的PGA达0.001g及以上的台站数据共有106个［6］.由于2011年新西兰Mw6.1地震的发震断层为隐伏型断层，目前已公开发表的数据均未给出距发震断层的距离（CESMD，2011），且观测数据相对集中，应用Boore等2007年提出的新一代衰减关系（NCA）［8］对强震观测数据进行分析过程中，利用台站到震中的距离近似为距发震断层的距离.因此，我们需要同时应用2010年新西兰Mw7.0地震的强震观测数据与NGA进行对比.2010年新西兰Mw7.0发生后，新西兰本土多处强震台站记录到了强震数据.截止到2011年5月CESMD公布的数据中，记录到的水平峰值加速度（PGA）达0.001g及以上的台站数据共有120个［7］.

分别利用现有的2010年新西兰Mw7.0地震的120个强震观测数据及2011年新西兰Mw6.1地震的106个强震观测数据，地震动衰减关系分别采用NGA.在Boore等2007年给出的NGA中，断层类型（走滑，正断层或逆冲断层）的影响及浅层速度（VS30）给予了充分考虑.图4分别给出了Mw7.0及Mw6.1地震走滑型断层的NGA及相应强震观测记录的对比分析.

图4a、4b和4c分别为不同浅层速度结构（VS30）条件下NGA与2010年新西兰Mw7.0地震强震观测PGA的对比关系图.图4b中99%以上的PGA观测记录（～200km）落在衰减关系曲线±95%的置信区间内，且距断层最近的台站PGA观测数据与其NGA对应程度最高.相比之下图4a和图4c中的NGA与2010年新西兰Mw7.0地震强震观测PGA则偏低或偏高，综合以上比较不难发现，对于2010年新西兰Mw7.0地震发震区域，VS30＝360m／s的NGA与强震观测记录更加符合，这一点与USGS给出的该区域VS30取值范围在300～360m／s相一致［16］.

图4 2010年新西兰Mw7.0地震及2011年新西兰Mw6.1地震强震观测数据与NGA对比图.（a），（b），（c）分别为2010年新西兰Mw7.0地震VS30＝760m／s、VS30＝360m／s和VS30＝180m／s PGA的NGA曲线；（d），（e），（f）分别为2011年新西兰Mw6.1地震VS30＝760m／s、VS30＝360m／s和VS30＝180m／s PGA的NGA曲线；（g），（h），（i）分别为2010年新西兰Mw7.0地震VS30＝760m／s、VS30＝360m／s和VS30＝180m／s PGV的NGA曲线Fig.4 Comparison of the observed PGAs and PGVs，New Zealand Earthquake of 2010，Mw7.0and 2011，Mw6.1with NGA.（a），（b），（c）are NGA of PGAs from New Zealand Earthquake of 2010，Mw7.0with VS30＝760m／s，VS30＝360m／s，and VS30＝180m／s，respectively；（d），（e），（f）are NGA of PGAs from New Zealand Earthquake of 2011，Mw6.1with VS30＝760m／s，VS30＝360m／s，and VS30＝180m／s，respectively；（g），（h），（i）are NGA of PGVs from New Zealand Earthquake of 2011，Mw6.1with VS30＝760m／s，VS30＝360m／s，and VS30＝180m／s，respectively

图4d、4e和4f分别为不同浅层速度结构（VS30）条件下NGA与2011年新西兰Mw6.1地震强震观测PGA的对比关系图，4g、4h和4i则分别为不同VS30取值下NGA与强震观测PGV的对比关系图.参照前面对图4a、4b和4c的分析，这里着重分析图4e和4h.图4e中PGA观测记录在≥30km的远场区域与VS30＝360m／s条件下NGA符合程度较好，约90%以上的观测数据在NGA±95%的置信区间内；尽管PGA观测记录在≤30km的近场区域与NGA的符合程度略有下降，但是值得注意的是，这些PGA的观测数据约有90%落在了NGA与＋95%的置信区间内.即2011年新西兰Mw6.1地震强震观测PGA在相对近场区域明显偏高，这一现象在图4h中更加明显.

2.4 建立有限断层模型模拟计算

最后，针对新西兰Mw6.1地震建立动态复合震源模型［17－18］，模型参数的设定综合考虑前文中提到的反演结果（图3a）、强震观测记录及理论模型，详见表2.作者在2011年针对2008年汶川地震的研究中［10］，修正了复合震源模型，提出可以根据已有反演数据，建立动态化的复合震源模型（DCSM），即断层面上的走向、倾向和滑动角参数不再单一取值，而是可以动态化的进行赋值.由于新西兰Mw6.1地震的主发震断层未破裂至地表，为隐伏型断层，因此，本文针对该地震建立的DCSM的动态化主要体现在滑动角的变化上（详见图2b、2c、2d）.

前文中提到理论的Brune圆盘模型下［9－10］，半径R≈5.0km，因此，可估算其静态应力降（Δσs）作为应用动态复合震源模型Δσd取值的参考：

表2 新西兰Mw6.1地震模型输入参数Table 2 Parameters for main fault of the dynamical composite source model，New Zealand，Mw6.1Earthquake

上式中M0为地震矩，参照表1取值，R取值5.0km，则有Δσs≈7.0MPa，而对应于圆盘模型的平均位移（¯D≈0.73RΔσs／μ）［19］，其中μ＝3·1011dyne·cm－2，则有¯D≈85cm，参照图2b，可视为Δσs的取值合理.具体的断层面上的几何和物理参数在表2中给出.

台站的强地面运动受到很多因素的影响，如场地条件、浅层速度结构及地形等因素等.HVSC台站位于新西兰Christchurch城内的一所学校内，尽管CESMD尚未给出该台站具体的场地条件，但其地形条件相对平坦，前文中在分析图4时已指出USGS给出的该区域VS30取值范围在300～360m／s，因此，HVSC台站很可能受到浅层速度结构（VS30）的影响［20－21］.本文在应用动态复合震源模型对该台站进行模拟计算的过程中，在考虑动态应力降Δσd取值的同时，用样参照图4d、e、f、g、h和i考虑了VS30不同取值的影响.其中Δσd的取值为三种情况，即Δσd＝σa≈4.5MPa、Δσd＝Δσs≈7.0MPa及Δσd＝σB≈19.35MPa，VS30的取值仍然为前文中提到的三种取值情况，即VS30＝760m／s、VS30＝360m／s、及VS30＝180m／s.模拟计算结果详见表3，表3中给出了Δσd和VS30分别取值下，应用文中所建立的动态复合震源模型计算HVSC台站峰值加速度（PGA）和峰值速度（PGV）的结果.前文中分析HVSC观测数据时提到，该台站垂向（Vertical）、S26W向和S64E向三个分量加速度的峰值分别为2160.5cm·s－2、1646.8cm·s－2及1244.6cm·s－2，速度的峰值分别为97.3cm·s－1、62.3cm·s－1及39.2cm·s－1（CESMD，2011）［6］.

分析表3给出的9种模拟结果可以发现，当VS30取值相同时，Δσd取值的变化对强地面运动质点加速度和速度的影响很大，以No.2、No.5和No.8三组模拟结果为例，其PGA的S26W、S64E和Vertical三分量之比约为1.5∶2.3∶6.2，同理PGV三分量的比值约为1.4∶2.3∶6.3.如果我们令Δσd取值相同，比较VS30取值变化对强地面运动的影响，以No.8模拟结果为例，其PGA和PGV三分量之比约分别为2.6∶3.2∶4.2和2.8∶3.0∶4.2.因此，HVSC台站强地面运动的PGA和PGV随着Δσd取值的变化程度要明显高于随着VS30取值的变化程度，同样以No.8模拟结果的PGA为例，其随着Δσd取值的变化最高幅度为47%，而随着VS30取值变化的最高幅度仅为16%.前文中针对强地面运动进行理论估算时已经指出，对于2011年新西兰Mw6.1地震，较高的Δσd很可能是影响近场强地面运动的主要原因之一，这一点与对表3的分析相一致.

图5中给出了HVSC台站强地面运动质点加速度三分量观测数据和模拟计算结果，及各自相应的Fourier频谱对比分析曲线，其中模拟结果采用的表3中No.8的结果.分别对比图5a和5d、5b和5e及5c和5f，可以发现模拟结果的时程曲线与实测记录具备较好的相似性，波形特征也较为相似，主要为高频成分且持续时间较短，大约为～10s；二者持续时间的一致性是比较高的，观测记录集中在10～20s之间，模拟结果则是在7～17s之间；模拟结果水平分量的振幅偏小，尤其是S64E分量，而垂向分量的振幅则略有偏大；进一步对比图5g、5h及5i中三分量模拟结果与观测记录Fourier频谱分别的对比图，整体一致性较高，Vertical方向的模拟结果在高频部分一致性要好于低频部分，而两个水平分量S26W和S64E，则低频部分的一致程度好于高频部分.

尽管图5中应用DCSM模拟HVSC台站强地面运动的结果与观测记录具备较好的一致性，但也存在一定的不足之处.本文所建立的DCSM作为运动学模型的一种采用的是1－D的速度结构，复杂地壳介质中（包括盆地效应，地形影响，高频地震波散射等）地震波传播过程不能完全地反映，远场强地面运动的模拟效果可能会偏离真实地表运动特征，尽管能够得到包含高频成分的波形模拟结果，但局限性是存在的.虽然实际发震断层的发震模式存在着不确定性和复杂性，应用动态复合震源模型模拟强地面运动也会存在一定的不可避免的问题，然而相比于点源模型仅能给出一条水平分量的时程曲线模拟结果，动力学模型仅能对低频段进行模拟，运动学模型存在一定的优势［17－18］.

3 讨论和结论

本文在引言中简单介绍了2011年新西兰Mw6.1地震具备典型的“震级小灾情重”的特征，造成这种现象主要有三个因素：（1）震源深度浅，通常发生在城市附近的地震，其震源深度越浅，则破坏力越大；（2）震中位置近，Mw6.1地震震中位于Christchurch西南方向不足10km处，且地震发生时正值人口聚集在市区的时间（新西兰当地时间中午12时51分）；（3）建筑物已脆弱，经历2010年Mw7.0主震及之后数百次的余震后，Christchurch城里的很多建筑物已经不堪一击，很多建筑物仍在重建，尚未完全修复，经不起二次打击（中国地震信息网，2011）［5］.

上述原因在一定层面上可以解释新西兰Mw6.1地震发震后造成破坏的原因，但没有从震源物理的角度说明一个6级的地震为什么能产生高达2.0g的近场质点加速度.本文在针对2011年新西兰Mw6.1地震的研究过程中，重点在于研究其震源的物理过程，及其地震发生后所造成的近场强地面运动的特征，尝试性地分析该地震发生后之所以产生近场强地面运动高值的原因.分别通过强震观测数据进行震源谱分析，并通过理论模型对其强地面运动的估算，以及最后通过建立动态复合震源模型进行模拟计算.

新西兰Mw6.1地震具备一个明显的特征，即其震源过程中有效应力降（动态应力降Δσd）偏高，将近20MPa.因此本文在研究过程中也从不同的角度对其进行了验证，结果表明，研究过程中Δσd的取值具备观测、理论及模型验证相一致的合理性.前文中通过对图5和表3的分析中可以发现，我们考虑了对近场强地面运动可能产生影响的两个方面Δσd及VS30，显然Δσd取值的变化对强地面运动的影响更明显，约为47%.针对Δσd取值的讨论，我们考虑了两个下限值，即令Δσd＝σa和Δσd＝Δσs，其中σa＝4.5MPa、Δσs＝7.0MPa.Savage和Wood（1971）［21］给定2σa／Δσs＜1或2σa／Δσs＞1对应于断层动态摩擦过程中的应力上调（overshoot）和下调（undershoot），显然新西兰Mw6.1地震属于应力下调，即分数应力降的情况［2］，这与前文中提到发震断层为隐伏型且能量震级（Me＝6.7）大于地震矩震级（Mw＝6.1）的反演结果相一致.

图5 HVSC台站地面质点加速度三分量模拟（Δσd＝19.35MPa，VS30＝360m／s）与实测时程曲线对比及Fourier频谱分析对比，灰色虚线为实测结果，黑色实线为模拟结果.（a）实测记录，Vertical分量，PGA为2160.5cm·s－2；（b）实测记录，S26W分量，PGA为1646.8cm·s－2；（c）实测记录，S64E分量，PGA为1244.6cm·s－2；（d）模拟结果，Vertical分量，PGA为2198.0cm·s－2；（e）模拟结果，S26W分量，PGA为1349.1cm·s－2；（f）模拟结果，S64E分量，PGA为607.4cm·s－2；（g）实测记录与模拟结果Fourier频谱对比图，Vertical分量；（h）实测记录与模拟结果Fourier频谱对比图，S26W分量；（i）实测记录与模拟结果Fourier频谱对比图，S64E分量Fig.5 Simulated（black real line），Δσd＝19.35MPa，VS30＝360m／s，and recorded（gray dash line）acceleration of threecomponent time－histories，Fourier spectra comparison at the station HVSC.（a）Vertical component，recording，the PGA is 2160.5cm·s－2；（b）S26Wcomponent，recording，the PGA is 1646.8cm·s－2；（c）S63Ecomponent，recording，the PGA is 1244.6cm·s－2；（d）Vertical component，simulation，the PGA is 2198.0cm·s－2；（e）S26Wcomponent，simulation，the PGA is 1349.1cm·s－2；（f）S63Ecomponent，simulation，the PGA is 607.4cm·s－2；（g）Fourier spectra comparison，Vertical component；（h）Fourier spectra comparison，S26Wcomponent；（i）Fourier spectra comparison，S64Ecomponent

如果选择M模型（M－model）［22－23］作为理论计算的模型，由于Madariaga在1976年的文章中指出，M－model的情况下Δσd总是小于Δσs［22］，则需进一步假设断层动态摩擦过程中对应于overshoot，即破裂过头的情况，则与前文中公式（7）相对应公式为［24］

同理上式中β为S波的波速，f0为拐角频率（f0＝0.22Hz），根据前文中公式（10）可推导出Δσs＝0.4375 M0／R3≈40MPa.

断层面上滑动位移的形式及大小等反演结果均存在不唯一性及不确定性，即随着Δσd的增大，震源有效滑动集中区的尺寸（size）则会变小，如果我们假设模型中Δσd设定偏小，则需要调高该值，这样一来，size则相应变小，这一点与前文提到的远近场不存在绝对的矛盾的观点相契合；另一方面，从理论上也对应到Brune理论模型［9－10］与Madariaga理论模型［22－23］中Δσd取值的不同之处，即同样size的震源有效破裂尺度，如果采用Madariaga的理论模型，Δσd的取值则要高于Brune理论的模型.

Shi等［25］2010年文章中关于Δσd与Δσs关系的讨论中指出，在M－model的情况下，当破裂速度（vR）与S波的波速（β）之比约为0.9时（即vR／β≈0.9，vR取值参见表2），则有Δσs／Δσd＝1.2，进一步可得到Δσs≈33.3MPa.根据前文中公式（9），可估算出M－model情况下质点加速度（ü）和质点速度）分别约为6666.7cm·s－2、333.3cm·s－1.也就是说，如果选择M－model作为理论计算的模型，则显然对新西兰Mw6.1地震近场强地面运动会造成过高的估计，相比之下，前文中所选择的Brune的圆盘模型则更为合理.

Pitarka等［26］针对隐伏型断层和破裂至地表断层的两类地震进行了对比研究，其文章中指出，隐伏型断层地震造成的强地面运动观测结果相对较高.Pitarka等［26］还指出，由于隐伏型断层未破裂至地表等原因，往往该类地震伴随着较高的应力降产生，并且较高的应力降被认为是该类地震近断层台站强地面运动变大的关键因素之一.该文章在结论中明确说明，对于震级相对较小的地震，沿发震断层方向距离震中30km内及沿垂直于发震断层方向距离断层10km范围内，隐伏型断层地震所造成的强地面运动明显大于破裂至地表断层地震的强地面运动.

2011年新西兰Mw6.1地震发震过程中Δσd相对较高的现象，与其隐伏型断层发震的地震机制具有一定的一致性，且观测到强地面运动最高值的HVSC台站距离震中仅1km.关于Δσd高值可能产生较强的强地面运动这一现象，本文在研究过程中也在观测数据分析、理论估计及模型计算三个方面给出了相对一致的论证.因此，我们认为新西兰Mw6.1地震近断层强地面运动偏高的主要原因在于，复杂震源过程中Δσd过高造成.

最后需要指出的是，未来的工作中应对可能发生的、距离城市较近的中小型地震给予足够的重视，相比于以往对城市活断层的重视，需加强对城市区域内隐伏断层的勘探，避免城市直下型地震的加强型破坏.并且需要重视城市周边震后的修复，防止二次破坏，重视可能遭受二次破坏的城市建筑，重视震后及时、快速、高标准修复，以便应对未来中小型地震（6级左右）可能对人口密集城市造成的破坏.

（References）

［1］ http：∥www.gns.cri.nz／Home／News－and－Events／Media－Releases／Most－damaging－quake－since－1931／Canterburyquake／Hidden－fault［2011－04－18］.

［2］ http：∥earthquake.usgs.gov／earthquakes／eqinthenews／2011／usc000466f／［2010－02－24］.

［3］ http：∥earthquake.usgs.gov／earthquakes／recenteqsww／Quakes／us2010atbj.php［2010－09－20］

［4］ 石耀霖，范桃园.大洋岩石层拖曳窄条陆壳俯冲的极限尺度分析：以新西兰南岛和大别山超高压变质带为例.地球物理学报，2001，44（6）：754－760.Shi Y L，Fan T Y.Analysis on maximum scale of continental crust sliver driven by subducting oceanic slab：a case study of South Island of New Zealand and UHPM at Dabie，China.Chinese J.Geophys.（in Chinese），2001，44（5）：754－760.

［5］ http：∥www.csi.ac.cn／manage／html／4028861611c5c2ba 0111c5c558b00001／＿content／11＿03／01／1298949448216.html［2011－03－01］

［6］ http：∥www.strongmotioncenter.org／cgi－bin／CESMD／iqr＿dist＿DM2.pl？IQRID＝NewZealand＿21Feb2011＆SFlag＝0＆Flag＝2［2011－02－24］

［7］ http：∥www.strongmotioncenter.org／cgi－bin／CESMD／iqr＿dist＿DM2.pl？iqrID＝NewZealand＿03Sep2010＆SFlag＝0＆Flag＝2［2010－09－09］

［8］ Boore D M，Atkinson G M.Boore－Atkinson NGA ground motion relations for the geometric mean horizontal component of peak and spectral ground motion parameters，PEER Report 2007／01.Pacific Earthquake Engineering Research Center，2007.

［9］ Brune J N.Tectonic stress and the spectra of seismic shear waves from earthquakes.J.Geophys.Res.，1970，75（26）：4997－5009.

［10］ Brune J N.Correction of“tectonic stress and spectra of seismic shear waves from earthquakes”.J.Geophys.Res.，1971，76（20）：5002.

［11］ 陈运泰.地震参数——数字地震在地震预测中的应用.北京：地震出版社，2003：73－81.Chen Y T.Earthquake Parameter：the Digital Seismic Application in Prediction the Earthquake（in Chinese）.Beijing：Earthquake Press，2003：73－81.

［12］ Andrews D J.Objective determination of source parameters and similarity of earthquakes of different size.∥Das S，Boatwright J，Scholz C H eds.Earthquake Source Mechanics.Washington：AGU，1986：259－267.

［13］ Thomas C H，Bakun W H.A bilinear source－scaling model for M－log A Observations of continental earthquakes.Bull.Seism.Soc.Amer.，2002，92（5）：1841－1846.

［14］ Wyss M，Brune J N.Seismic moment，stress，and source dimensions for earthquakes in the California－Nevada region.J.Geophys.Res.，1968，73（14）：4681－4694.

［15］ Brune J N.The Physics of Earthquake Strong Motion.∥Lomnitz C，Rosenblueth E eds.Seismic Risk and Engineering Decisions.New York：Elsevier Sci.，Publ.Co.Press，1976：141－177.

［16］ http：∥ehp1－earthquake.cr.usgs.gov／hazards／apps／vs30／output／HVSC／1313993824.jpg［2010－09－20］

［17］ Zeng Y，Anderson J G，Yu G.A composite source model for computing realistic synthetic strong ground motions.Geophys.Res.Lett.，1994，21（8）：725－728.

［18］ 孟令媛，史保平.应用动态复合震源模型模拟汶川Mw7.9地震强地面运动.地球物理学报，2011，54（4）：1010－1027.Meng L Y，Shi B P.Near－fault strong ground motion simulation of the May 12，2008，Mw7.9Wenchuan earthquake by dynamical composite source model.Chinese J.Geophys.（in Chinese），2011，54（4）：1010－1027.

［19］ Anderson J G.Seismic energy and stress－drop parameters for a composite source model.Bull.Seism.Soc.Amer.，1997，87（1）：85－96.

［20］ Wald D J，Allen T I.Topographic slope as a proxy for seismic site conditions and amplification.Bull.Seism.Soc.Amer.，2007，97（5）：1379－1395.

［21］ Savage J C，Wood M D.The relation between apparent stress and stress drop.Bull.Seism.Soc.Amer.，1971，61（5）：1381－1388.

［22］ Madariaga R.Dynamics of an expanding circular fault.Bull.Seism.Soc.Am.，1976，66（3）：639－666.

［23］ Boatwright J L.Spectral theory for circular seismic sources；Simple estimates of source dimension，dynamic stress drop，and radiated seismic energy.Bull.Seism.Soc.Am.，1980，70（1）：1－27.

［24］ Aki K，Richards F G.Quantitative Seismology Theory and Methods.New York：W H Freeman and Company Press，1980：883.

［25］ Shi B P，Liu B Y，Meng L Y.Bounding of near－fault ground motion based on radiated seismic energy with a consideration of fault frictional mechanisms.Earthquake Science，2010，23（4）：357－368.

［26］ Pitarka A，Dalguer L A，Day S M，et al.Numerical study of ground－motion differences between buried－rupturing and surface－rupturing earthquake.Bull.Seism.Soc.Amer.，2009，99（3）：1521－1537.

Source rupture process of the February 21，2011，Mw6.1，New Zealand earthquake and the characteristics of the near－field strong ground motion

MENG Ling－Yuan1，2，SHI Bao－Ping1
1 Graduate University of Chinese Academy of Sciences，Beijing100049，China
2 China Earthquake Networks Center，Beijing100036，China

The New Zealand earthquake of February 21，2011，Mw6.1occurred in the South Island，New Zealand with the epicenter at longitude 172.70°E and latitude 43.58°S，and with depth of 5km.The Mw6.1earthquake occurred on an unknown blind fault involving obliquethrust faulting，and this fault is 9km south of Christchurch，which is the third largest city of New Zealand，with a striking direction from east toward west.At least 163were killed and lots of structures were destroyed in Christchurch city.The Peak Value of Acceleration（PGA）observed at station Heathcote Valley Primary School（HVSC），which is 1km away from the epicenter，is up to almost 2.0g.The ground－motion observation suggests that the New Zealand earthquake generated higher than normal near－fault ground motion.In this study，we analyzed the source spectral parameters based on the observations，and estimated the near－field groundmotion based on the Brune′s circle model.The results indicate that the larger ground motion may be produced by higher effective stress drop（dynamic stress drop，Δσd）in this buried earthquake.A dynamical composite source model（DCSM）has been developed to simulate the near－fault strong ground motion with associated fault rupture properties from a kinematic point of view.For comparing purpose，we conduct broadband ground motion predictions for the station of HVSC；the synthetic seismogram produced for this station has good agreement with the observations in time histories，waveforms，peak values and frequency contents，which indicate that one of the key factors in ground－motion amplification was the higher effective stress drop.It should be pointed out that，in addition to the study of large earthquakes，the small and moderate earthquakes near the cities should be monitored in order to prevent the devastating destroy of the cities once earthquake occurs.

New Zealand Mw6.1earthquake，Observation strong motions，Effective stress drop，Strong ground motion，Dynamical composite source model

10.6038／j.issn.0001－5733.2012.05.018

P315

2011－08－23，2012－03－28收修定稿

中国科学院知识创新工程（KZCX2－YW－Q08－2）资助.

孟令媛，女，1983年生，博士，主要从事强地面运动和震源动力学研究.E－mail：meng.lingyuan＠hotmail.com

孟令媛，史保平.2011年新西兰Mw6.1地震震源过程及强地面运动特征初步分析.地球物理学报，2012，55（5）：1601－1612，

10.6038／j.issn.0001－5733.2012.05.018.

Meng L Y，Shi B P.Source rupture process of the February 21，2011，Mw6.1，New Zealand earthquake and the characteristics of the near－field strong ground motion.ChineseJ.Geophys.（in Chinese），2012，55（5）：1601－1612，doi：10.6038／j.issn.0001－5733.2012.05.018.

（本文编辑 胡素芳）