2013年蘆山MS7．0地震震源參數(shù)特征及近斷層強(qiáng)地面運(yùn)動(dòng)初步估計(jì)*

孟令媛 周龍泉 劉 杰

（中國(guó)北京 100045 中國(guó)地震臺(tái)網(wǎng)中心）

引言

2013年4月20日8時(shí)2分在我國(guó)四川省雅安市蘆山縣發(fā)生了MS7．0地震，震中位置為30．3°N、103．0°E．蘆山MS7．0地震距2008年5月12日汶川MS8．0地震初始破裂震中約85km，距汶川地震余震密集分布區(qū)南端約60km．蘆山MS7．0地震震中位于龍門(mén)山斷裂帶南段的前山斷裂附近，主震破裂過(guò)程呈北北東向的擠壓逆沖變形特征，與龍門(mén)山斷裂帶總體走向和運(yùn)動(dòng)性質(zhì)一致（龍門(mén)山斷裂帶走向主要為近北東向）．初步分析認(rèn)為，蘆山MS7．0地震與汶川MS8．0地震均為巴顏喀拉地塊向東運(yùn)動(dòng)遇到華南地塊阻擋，應(yīng)力積累和釋放的結(jié)果（劉杰等，2013）．龍門(mén)山斷裂帶全長(zhǎng)約500km，寬約30—40km．其橫向主要由中央斷裂帶、山后斷裂帶、山前斷裂帶及推覆構(gòu)造帶組成；縱向可以分為3段：北川—寧強(qiáng)、勉縣為東北段；北川—都江堰為中段；都江堰—瀘定、康定附近為西南段．這次蘆山MS7．0地震發(fā)生在西南段，2008年汶川MS8．0地震則發(fā)生在中北段，汶川地震的發(fā)生使得龍門(mén)山斷裂帶整體由震中—東北方向破裂了約300km．

蘆山MS7．0地震震中附近100km范圍內(nèi)，1900年以來(lái)共發(fā)生MS≥5．0地震12次，其中MS6．0—6．9地震3次，最大地震為2008年汶川MS8．0地震①http：∥10．5．160．59／netoffice／module／info／portal／custom／dzj／dzj＿news＿display．jsp？GeneralID＝12007．．截止到2013年4月26日，蘆山MS7．0地震造成196人遇難，21人失蹤，逾萬(wàn)人受傷．其中破壞最嚴(yán)重的發(fā)生在震中附近的蘆山、寶興等地區(qū)．震區(qū)余震不斷，截止到4月26日12時(shí)，蘆山地震余震區(qū)已記錄MS≥3．0余震110次．其中MS≥5．0余震4次，最大余震為4月21日17時(shí)5分MS5．4地震，震中位于蘆山縣與邛崍市交界（圖1）．

本研究旨在分析和討論蘆山MS7．0地震的震源參數(shù)特征，進(jìn)而估計(jì)其近斷層區(qū)域強(qiáng)地面運(yùn)動(dòng)的特征和大小．基于已知的震源參數(shù)計(jì)算地震視應(yīng)力等震源參數(shù)，明確蘆山MS7．0地震為震源動(dòng)態(tài)模型中的應(yīng)力下調(diào)模式；結(jié)合斷層面上滑動(dòng)集中區(qū)的分布，進(jìn)一步構(gòu)建動(dòng)態(tài)復(fù)合震源模型（dynamical composite source model，簡(jiǎn)寫(xiě)為DCSM），選取極震區(qū)寶興和蘆山兩個(gè)特征點(diǎn)模擬計(jì)算其強(qiáng)地面運(yùn)動(dòng)參數(shù)．初步結(jié)果顯示，蘆山MS7．0地震的近斷層區(qū)域強(qiáng)地面運(yùn)動(dòng)受有效應(yīng)力降及斷層面上滑動(dòng)集中區(qū)的影響較為明顯．

1 震源參數(shù)及強(qiáng)地面運(yùn)動(dòng)特征

Wyss和Brune（1968）提出視應(yīng)力σa的概念，視應(yīng)力是與地震輻射能量直接相關(guān)的物理量．其物理意義在于，發(fā)震斷層單位面積發(fā)生單位錯(cuò)動(dòng)輻射的地震波能量的大小，即

式中，M0為地震矩，單位為N·m；μ為介質(zhì)的剪切模量，通常取值為μ＝3×104N·m．ES為地震波輻射能，美國(guó)地質(zhì)調(diào)查局（USGS）最新公布的測(cè)定結(jié)果為（2．3—2．8）×1014J．結(jié)合式（1）和表1，可以求得蘆山MS7．0地震的視應(yīng)力σa約為0．49—0．84MPa．

圖1 蘆山MS7．0地震主震位置（紅色星號(hào)）及MS≥3．0余震（綠色圓）分布（截止到2013年4月26日）Fig．1 Locations of main shock（red star）and MS≥3．0aftershocks（green circles）of the Lushan MS7．0earthquake，where the aftershocks are from 20April 2013to 26April 2013

表1 蘆山MS7．0地震震源及斷層參數(shù)Table 1 Source and fault parameters of the Lushan MS7．0earthquake

陳學(xué)忠等（2003）指出，視應(yīng)力σa與余震強(qiáng)度存在著一定的統(tǒng)計(jì)關(guān)系，即通常一次主震后高視應(yīng)力的地震一般對(duì)應(yīng)后續(xù)較大的余震產(chǎn)生．基于對(duì)造成視應(yīng)力取值大小的震源力學(xué)過(guò)程及斷層破裂過(guò)程中動(dòng)摩擦機(jī)制影響的考慮，以及單一σa取值大小對(duì)余震危險(xiǎn)性進(jìn)行評(píng)估存在一定的不確定性，本研究引入地震發(fā)生時(shí)動(dòng)態(tài)破裂過(guò)程模型，結(jié)合視應(yīng)力σa的取值，對(duì)蘆山MS7．0地震斷層破裂過(guò)程進(jìn)行分析．在不考慮破裂能（Er）的前提下，可以將震源動(dòng)態(tài)模型過(guò)程分為3種主要的模型：完全應(yīng)力降模型（Orowan，1960；Brune，1970）、應(yīng)力上調(diào)模型（Savage，Wood，1971）和應(yīng)力下調(diào)模型（Brune，1970，1976；Smithetal，1991；Zú?iga，1993）．這3種模型分別對(duì)應(yīng)摩擦應(yīng)力等于、大于和小于斷層面上的終止剪切應(yīng)力的情況．結(jié)合斷層破裂過(guò)程中靜態(tài)應(yīng)力降（Δσs），可以表示為：當(dāng)2σa＝Δσs時(shí)，為完全應(yīng)力降模型；當(dāng)2σa＜Δσs時(shí)，為應(yīng)力上調(diào)模型；當(dāng)2σa＞Δσs時(shí)，為應(yīng)力下調(diào)模型．

對(duì)于已知斷層面上滑動(dòng)位移，靜態(tài)應(yīng)力降（Δσs）可表示為（Starr，1928；Knopoff，1958；Keilis-Borok，1959）

式中，D為斷層面上的平均滑動(dòng)位移，W為斷層寬度．結(jié)合表1，取D≈1．0m，W≈30km．參照式（2）和式（3），Δσs約為0．6—0．8MPa．

結(jié)合已知斷層面上滑動(dòng)位移的分布情況來(lái)估計(jì)蘆山MS7．0地震的有效破裂面積．首先應(yīng)用矩震級(jí)與斷層破裂面積對(duì)數(shù)之間的雙線性經(jīng)驗(yàn)關(guān)系，估算發(fā)震斷層破裂面積（Hanks，Bakun，2002）

式中，MW為矩震級(jí)，A為斷層破裂面積．式（5）在Hanks和Bakun（2002）原文中為分段函數(shù)，即以A取537km2為分段點(diǎn)，當(dāng)A＝537km2時(shí)，MW≈6．71．由于本文涉及的蘆山MS7．0地震矩震級(jí)小于6．71（表1），因此我們?cè)隍?yàn)算該地震的破裂面積時(shí)，僅采用式（5）中給出的這一部分，進(jìn)一步可以推算出蘆山MS7．0地震的有效破裂面積約為389—447km2．對(duì)比圖2與表1，圖2a中滑動(dòng)位移大于0．2m的有效破裂面積約為450km2，圖2b中滑動(dòng)位移大于0．3m的有效破裂面積約為400—500km2．事實(shí)上，斷層面上滑動(dòng)位移的反演結(jié)果并不具備唯一性（表1），式（5）的推算結(jié)果與圖2中的結(jié)果具備相對(duì)一致性．由于蘆山地震斷層面上滑動(dòng)分布集中且相對(duì)單一，研究過(guò)程中可以將其近似為圓盤(pán)模型進(jìn)行理論計(jì)算．令圓盤(pán)模型半徑R≈11km，結(jié)合式（4）則Δσs約為0．36MPa．

圖2 蘆山MS7．0地震主斷層面上滑動(dòng)位移分布情況（a）陳運(yùn)泰等的模型反演結(jié)果① http：∥www．cea-igp．a(chǎn)c．cn／upload／Image／mrtp／tpxw／20130420ds／2668196836．jpg．；（b）本文根據(jù)DCSM模型給出的滑動(dòng)位移不均勻分布特征Fig．2 Slip-displacement distribution on the main fault plan for the Lushan MS7．0earthquake（a）The inversion result from Chen et al①；（b）Inhomogeneous distribution of slipon the fault used in our DCSM model

基于式（2）—（4）的計(jì)算結(jié)果，已知Δσs取值范圍約為0．36—0．8MPa，結(jié)合σa的估計(jì)范圍為0．49—0．84MPa，可以得出蘆山地震為2σa＞Δσs情況下的應(yīng)力下調(diào)模式的初步結(jié)論．Brune（1970）認(rèn)為，應(yīng)力下調(diào)模式對(duì)應(yīng)著地震時(shí)斷層錯(cuò)動(dòng)會(huì)突然受阻而可能發(fā)生的被突然鎖住的情況．對(duì)比2008年汶川MS8．0地震，基于USGS給出的震源參數(shù)及式（2），可以計(jì)算得到汶川MS8．0地震的σa僅為0．55MPa，屬于應(yīng)力上調(diào)模式．該模式對(duì)應(yīng)地震時(shí)的斷層錯(cuò)動(dòng)過(guò)頭情況．盡管2008年汶川MS8．0地震與2013年蘆山MS7．0地震同發(fā)生在龍門(mén)山斷裂帶上，但二者σa大小存在較大差異，且斷層破裂過(guò)程中對(duì)應(yīng)的兩種錯(cuò)動(dòng)模式完全不同．

截止到2013年4月19日，汶川MS8．0地震余震區(qū)共發(fā)生MS≥5．0余震50次，最大余震為MS6．4；截止到4月26日，蘆山MS7．0地震余震區(qū)已記錄MS≥5．0余震4次，最大余震為MS5．4．汶川地震與蘆山地震余震的時(shí)空分布特征明顯不同：汶川地震具備明顯的分段性，盡管其主震發(fā)震位置在龍門(mén)山斷裂帶南西段的汶川—映秀一段，但其余震展布則呈明顯的沿該斷裂帶北東段的青川一段分布，長(zhǎng)達(dá)300km；蘆山地震斷層破裂長(zhǎng)度較短，震后余震展布較汶川地震明顯集中程度較高．

1．1 估算近斷層區(qū)域強(qiáng)地面運(yùn)動(dòng)的理論值

Andrews（1986）對(duì)地震強(qiáng)地面運(yùn)動(dòng)的頻譜響應(yīng)進(jìn)行分析時(shí)給出

式中，σB為Brune圓盤(pán)模型下的有效應(yīng)力降；AFS為S波輻射圖型因子，其在震源球面上的均方根為AFS＝．由前面已給出的計(jì)算結(jié)果σa≈0．49—0．84MPa，則σB≈2．1—3．6 MPa．進(jìn)一步采取Brune（1970，1971）圓盤(pán)模型下給出的近場(chǎng)強(qiáng)震質(zhì)點(diǎn)運(yùn)動(dòng)速度的估算方法，即

式中，ü和˙u分別為質(zhì)點(diǎn)加速度和質(zhì)點(diǎn)速度，則有ü≈420—720cm／s2，˙u≈21—36cm／s．

蘆山MS7．0地震發(fā)生后，中國(guó)地震局于4月27日給出了實(shí)測(cè)烈度分布圖，極震區(qū)最大烈度為Ⅸ度．烈度通常是通過(guò)對(duì)震區(qū)實(shí)際調(diào)查獲取的震害資料，由于地震發(fā)生后的影響區(qū)域很大，往往是震后一段時(shí)間才能給出現(xiàn)場(chǎng)震害調(diào)查的烈度分布圖．然而在地震應(yīng)急救援中，為快速可靠地估計(jì)極震區(qū)的烈度值，可以利用烈度與強(qiáng)地面運(yùn)動(dòng)參數(shù)之間的轉(zhuǎn)換關(guān)系（Waldetal，1999）

式中，PGA單位為cm／s2，PGV單位為cm／s，IMM為場(chǎng)地烈度．將基于式（7）得到的近場(chǎng)強(qiáng)地面運(yùn)動(dòng)的理論估計(jì)值，近似為峰值加速度（PGA）和峰值速度（PGV），代入式（8）中，IMM≈7．0—9．9．該結(jié)果與實(shí)測(cè)的極震區(qū)最大烈度Ⅸ較為接近．

1．2 構(gòu)建動(dòng)態(tài)復(fù)合震源模型計(jì)算近斷層區(qū)域強(qiáng)地面運(yùn)動(dòng)

針對(duì)蘆山MS7．0地震構(gòu)建動(dòng)態(tài)復(fù)合震源模型（DCSM）（孟令媛，史保平，2011），模型參數(shù)的設(shè)定綜合考慮前面提到的反演結(jié)果（表1、圖2a）．作者在2011年針對(duì)汶川MS8．0地震的研究中，修正了復(fù)合震源模型，提出可以根據(jù)已有反演數(shù)據(jù)建立動(dòng)態(tài)化的復(fù)合震源模型（DCSM），即斷層面上的走向、傾向和滑動(dòng)角參數(shù)不再單一取值，而是可以動(dòng)態(tài)化地進(jìn)行賦值．此外，DCSM模型修正了原有的試錯(cuò)算法，應(yīng)用地震矩和地震波輻射能（有效應(yīng)力降或視應(yīng)力）兩個(gè)條件同時(shí)約束模型．盡管目前尚無(wú)查到有關(guān)蘆山MS7．0地震的速度結(jié)構(gòu)文獻(xiàn)，由于蘆山MS7．0地震與汶川MS8．0地震同處龍門(mén)山斷裂帶，因此，本文針對(duì)蘆山MS7．0地震建立DCSM的相關(guān)速度結(jié)構(gòu)主要參照汶川MS8．0地震（劉啟元等，2009），其斷層面上滑動(dòng)位移分布見(jiàn)圖2b．

圖3 應(yīng)用DCSM模型給出的質(zhì)點(diǎn)運(yùn)動(dòng)加速度（a）、速度（b）三分量時(shí)程曲線（寶興特征點(diǎn)）Fig．3 Simulated acceleration（a）and velocity（b）time-histories in the NS，EW and UD directions，respectively，at Baoxing

圖4 應(yīng)用DCSM模型給出的質(zhì)點(diǎn)運(yùn)動(dòng)加速度（a）、速度（b）三分量時(shí)程曲線（蘆山特征點(diǎn)）Fig．4 Simulated acceleration（a）and velocity（b）time-histories in the NS，EW and UD directions，respectively，at Lushan

應(yīng)用已構(gòu)建的DCSM模型，對(duì)近斷層區(qū)域的寶興和蘆山兩個(gè)特征點(diǎn)進(jìn)行了模擬計(jì)算，并分別給出了加速度和速度的模擬時(shí)程曲線圖（圖3、圖4）．挑選寶興和蘆山兩個(gè)特征點(diǎn)的主要原因在于這兩個(gè)特征點(diǎn)位處Ⅷ— Ⅸ度烈度的極震區(qū)，距離蘆山地震的震中位置及主斷層均較近．寶興和蘆山兩個(gè)特征點(diǎn)距離主斷層的最短距離分別為7和15km，即均為近斷層區(qū)域（＜30km），二者的震中距分別為19和17km．此外，由于寶興和蘆山兩地為震后人員傷亡的集中地區(qū)，因此本研究確定兩個(gè)特征點(diǎn)的經(jīng)緯度位置時(shí)主要參考了兩個(gè)地區(qū)人口集中點(diǎn)的經(jīng)緯度．國(guó)家強(qiáng)震動(dòng)臺(tái)網(wǎng)中心針對(duì)蘆山MS7．0地震公布的強(qiáng)震數(shù)據(jù)統(tǒng)計(jì)情況，其初步獲取的95組三分量加速度記錄中，震中距范圍為27．7—769km．獲取的PGA最大為400．7cm／s2，對(duì)應(yīng)震中距為27．7km．已獲取記錄中PGA大于200cm／s2的記錄共8條．

圖3給出了寶興特征點(diǎn)地表質(zhì)點(diǎn)加速度、速度三分量數(shù)值模擬計(jì)算結(jié)果．圖中給出的分別為地表質(zhì)點(diǎn)運(yùn)動(dòng)NS、EW及UD三分量加速度和速度模擬結(jié)果的時(shí)程曲線圖．由圖3可以看出，模擬結(jié)果的加速度和速度的各自三分量時(shí)程曲線圖均呈現(xiàn)波形集中和高頻成分偏多的特點(diǎn)．波形集中即時(shí)程曲線圖的持續(xù)時(shí)間相對(duì)較短，僅約為5—10s，這一點(diǎn)與寶興臺(tái)站距離震中位置較近而接收到的地震波成分中高頻成分較多是趨于一致的．寶興特征點(diǎn)模擬的NS、EW及UD三分量的PGA分別為389，500和331cm／s2，PGV模擬結(jié)果分別為25，43和31cm／s．由于現(xiàn)有DCSM模型近似采用了汶川MS8．0地震的速度結(jié)構(gòu)（劉啟元等，2009），針對(duì)寶興特征點(diǎn)的模擬結(jié)果是基于基巖上獲得的，并未添加土層的速度數(shù)據(jù)，因此圖3中強(qiáng)地面運(yùn)動(dòng)參數(shù)模擬結(jié)果中高頻的成分相對(duì)較多，尤其在加速度時(shí)程曲線圖中更為明顯．

圖4給出了蘆山特征點(diǎn)地表質(zhì)點(diǎn)加速度、速度NS、EW及UD三分量模擬結(jié)果的時(shí)程曲線圖．蘆山特征點(diǎn)模擬的NS、EW及UD三分量PGA分別為194，205和289cm／s2，PGV模擬結(jié)果分別為20，26和16cm／s．圖4給出的蘆山特征點(diǎn)加速度和速度的三分量時(shí)程曲線的模擬結(jié)果相較于圖3呈現(xiàn)出部分低頻成分，尤其是速度三分量的時(shí)程曲線，且持續(xù)時(shí)間約為10s．蘆山特征點(diǎn)PGA和PGV均比寶興特征點(diǎn)?。畤?guó)家強(qiáng)震動(dòng)臺(tái)網(wǎng)在已公布PGA的三分量加速度記錄中，PGA第二高值出現(xiàn)在距離震中32．6km的蘆山飛仙臺(tái)，該臺(tái)站實(shí)測(cè)的NS、EW及UD三分量PGA分別為357．0，387．4和267．4cm／s2．由于目前尚未公布相關(guān)PGV實(shí)測(cè)記錄及相關(guān)時(shí)程曲線圖，且本研究選取的蘆山特征點(diǎn)與蘆山飛仙臺(tái)站位置有所差別，因此，圖4中給出的模擬結(jié)果與蘆山飛仙臺(tái)的實(shí)測(cè)記錄并不具有完全的可比性．盡管如此，仍可以依據(jù)蘆山飛仙臺(tái)的PGA實(shí)測(cè)記錄，判斷本研究基于DCSM模型針對(duì)所選取特征點(diǎn)強(qiáng)地面運(yùn)動(dòng)的模擬結(jié)果具有較好的合理性．

本研究給出的強(qiáng)地面運(yùn)動(dòng)模擬結(jié)果不包含儀器響應(yīng)差異的影響，相關(guān)程序計(jì)算地表運(yùn)動(dòng)三分量的過(guò)程主要為選定速度結(jié)構(gòu)、給出滑動(dòng)模型、輸入相關(guān)計(jì)算參數(shù)后，DCSM模型即為確定．實(shí)際上，模型一旦確定，針對(duì)不同特征點(diǎn)（某地點(diǎn)）計(jì)算得到三分量的時(shí)程曲線圖即為一次性計(jì)算得到．由于格林函數(shù)的計(jì)算結(jié)果需要在球坐標(biāo)系與笛卡爾坐標(biāo)系中進(jìn)行轉(zhuǎn)換，程序設(shè)定的默認(rèn)方向即為NS、EW和UD三個(gè)方向，因此一次性的計(jì)算結(jié)果被分解為這三個(gè)方向．DCSM的優(yōu)勢(shì)之一即可以給出含有高頻成分的地震到模擬結(jié)果．與同樣能給出含高頻成分的有限斷層隨機(jī)振動(dòng)模型不同，DCSM模型考慮應(yīng)力降、地震波輻射能的影響，并且能夠給出垂直分量的時(shí)程曲線圖．

2 討論與結(jié)論

2013年4月20日蘆山地震發(fā)生后，關(guān)于蘆山MS7．0地震與汶川MS8．0地震的關(guān)聯(lián)性問(wèn)題立即引起了學(xué)術(shù)界的廣泛討論．從構(gòu)造背景的角度考慮，蘆山地震與汶川地震的發(fā)震位置均為巴顏喀拉塊體的邊界上，蘆山地震與汶川地震的震中位置距離約85km，且蘆山地震與汶川地震均發(fā)生在龍門(mén)山斷裂帶上，受到龍門(mén)山斷裂帶總體走向和運(yùn)動(dòng)特征的影響，兩次地震的震源機(jī)制解也較為一致，均為逆沖型的地震．前面也提到，蘆山MS7．0地震為2σa＞Δσs情況下的應(yīng)力下調(diào)模式，對(duì)應(yīng)著地震時(shí)斷層錯(cuò)動(dòng)會(huì)突然受阻而可能發(fā)生被突然鎖住的情況；而汶川地震則屬于應(yīng)力上調(diào)模式，對(duì)應(yīng)地震時(shí)斷層錯(cuò)動(dòng)過(guò)頭的情況．由于二者σa大小存在較大差異，且斷層破裂過(guò)程中對(duì)應(yīng)兩種完全不同的錯(cuò)動(dòng)模式，不同的錯(cuò)動(dòng)模式反映了斷層面上應(yīng)力變化的不同情況，對(duì)強(qiáng)地面運(yùn)動(dòng)地震學(xué)意義重大，因其決定了高頻輻射能的相對(duì)數(shù)量，進(jìn)而決定了給定應(yīng)力降之后的峰值加速度和峰值速度模擬結(jié)果的差異性．

蘆山MS7．0地震具備一個(gè)特征，即其震源過(guò)程中有效應(yīng)力降與汶川MS8．0地震的差異性并不大．主要原因在于兩次地震的σa差異較小，蘆山MS7．0地震σa為0．49—0．84 MPa，汶川MS8．0地震σa為0．55MPa，有效應(yīng)力降σB的差異性也不明顯．由此基于Brune圓盤(pán)模型對(duì)蘆山地震近場(chǎng)質(zhì)點(diǎn)加速度和速度進(jìn)行估算可以發(fā)現(xiàn)，盡管蘆山地震斷層破裂尺度相對(duì)于汶川地震明顯偏小，但其近場(chǎng)強(qiáng)地面運(yùn)動(dòng)的理論估算值仍相對(duì)較高，這一點(diǎn)可以利用極震區(qū)烈度值的大小得到驗(yàn)證．

關(guān)于發(fā)震斷層破裂過(guò)程的模擬，本文基于斷層面上的滑動(dòng)位移分布的反演結(jié)果約束斷層面上滑動(dòng)集中區(qū)的分布位置（圖2a），基于大小不同地震破裂過(guò)程具有自相似的假定（Frankel，1991），即主震斷層面可由多個(gè)隨機(jī)分布的大小尺度不同的子源疊加而成，隨機(jī)疊加而成的斷層面上的滑移分布遵從k－2模型（Zeng，Anderson，1996），進(jìn)而完成對(duì)斷層破裂過(guò)程的運(yùn)動(dòng)學(xué)描述（圖2b）．現(xiàn)階段DCSM的改進(jìn)程度仍然存在一定的局限性，主要在于對(duì)最大一二個(gè)子源位置的約束．實(shí)際上，如果DCSM模型構(gòu)建過(guò)程中不約束子源位置，完全隨機(jī)分配，對(duì)模擬結(jié)果整體PGA和PGV的平均水平不會(huì)有大的影響，但對(duì)近斷層區(qū)域強(qiáng)地面運(yùn)動(dòng)的模擬會(huì)影響很大，尤其是對(duì)最大子源位置附近區(qū)域特征點(diǎn)模擬結(jié)果的PGA和PGV影響很大．

關(guān)于DCSM模型相關(guān)程序計(jì)算三分量的具體過(guò)程，主要為選定速度結(jié)構(gòu)、給出滑動(dòng)模型和輸入相關(guān)計(jì)算參數(shù)三個(gè)步驟．模型一旦確定，針對(duì)不同特征點(diǎn)的計(jì)算需要調(diào)整的參數(shù)只有特征點(diǎn)的經(jīng)緯度參數(shù)．根據(jù)不同的經(jīng)緯度參數(shù)，會(huì)針對(duì)該特征點(diǎn)重新計(jì)算一次格林函數(shù)．同樣地，一旦確定某一個(gè)特征點(diǎn)，計(jì)算得到三分量的時(shí)程曲線圖即為一次性計(jì)算得到．DCSM模型計(jì)算程序設(shè)定的默認(rèn)方向即為NS、EW和UD三個(gè)方向，因此一次性的計(jì)算結(jié)果被分解為這三個(gè)方向．

DCSM的前身為復(fù)合震源模型，創(chuàng)新性地采用了隨機(jī)子源分配的方式，解決了地震波時(shí)程曲線圖模擬結(jié)果中的高頻成分問(wèn)題．同樣地，DCSM的優(yōu)勢(shì)之一即可以給出含有高頻成分的地震動(dòng)模擬結(jié)果．DCSM模型考慮應(yīng)力降、地震波輻射能的影響，能夠分別給出水平和垂直分量的時(shí)程曲線圖．原始的復(fù)合震源模型通常針對(duì)小尺度的走滑型斷層進(jìn)行模擬，而DCSM模型改進(jìn)后不僅可以針對(duì)走滑型斷層進(jìn)行模擬，也可以針對(duì)逆沖型斷層模型進(jìn)行模擬．除考慮應(yīng)力降、地震波輻射能的影響外，還實(shí)現(xiàn)了分段斷層走向、傾向的滑動(dòng)角的動(dòng)態(tài)化設(shè)定．但現(xiàn)有DCSM模型所給出的強(qiáng)地面運(yùn)動(dòng)模擬結(jié)果并未包含儀器響應(yīng)的差異．隨著蘆山MS7．0地震強(qiáng)震觀測(cè)數(shù)據(jù)和震區(qū)速度結(jié)構(gòu)的研究的不斷完善，未來(lái)還可以對(duì)已構(gòu)建的DCSM模型進(jìn)行逐步地修正和改進(jìn)，以求獲取適用于龍門(mén)山斷裂帶強(qiáng)震地表運(yùn)動(dòng)預(yù)測(cè)的有效模型．

2011年2月21日發(fā)生在新西蘭克萊斯特徹奇的MW6．1地震，造成200多人遇難．相對(duì)該地區(qū)半年前發(fā)生的MW7．0地震，這次地震給克萊斯特徹奇城造成了更為嚴(yán)重的建筑物損毀和人員傷亡．一個(gè)重要的原因就在于這次MW6．1地震造成的近斷層強(qiáng)地面運(yùn)動(dòng)強(qiáng)且震中位處城市中心附近，加之半年前MW7．0地震對(duì)城市建筑物的損害，克萊斯特徹奇城內(nèi)很多建筑物仍在重建，尚未完全修復(fù)，經(jīng)不起二次打擊（孟令媛，史保平，2012）．因此，未來(lái)的研究應(yīng)對(duì)可能發(fā)生的、距離城市較近的大中型地震給予足夠的重視，避免城市直下型地震的加強(qiáng)型破壞；需要加強(qiáng)地震發(fā)生后對(duì)城市內(nèi)建筑物可能產(chǎn)生強(qiáng)地面運(yùn)動(dòng)影響的判斷，尤其是人口密集的城市以及曾經(jīng)受過(guò)地震破壞的城市建筑，重視震后及時(shí)、快速、高標(biāo)準(zhǔn)修復(fù)，以便應(yīng)對(duì)未來(lái)地震可能對(duì)人口密集城市造成的破壞．

本文在撰稿過(guò)程中，史保平教授、蔣海昆研究員和張永仙研究員給予了有益的指導(dǎo)與幫助，作者在此一并表示誠(chéng)摯的謝意．

陳學(xué)忠，王小平，王林瑛，張?zhí)熘校?003．地震視應(yīng)力用于震后趨勢(shì)快速判定的可能性［J］．國(guó)際地震動(dòng)態(tài)，（7）：1-3．

劉杰，易桂喜，張致偉，官致君，阮祥，龍鋒，杜方．2013．2013年4月20日四川蘆山M7．0級(jí)地震介紹［J］．地球物理學(xué)報(bào)，56（4）：1404-1407．

劉啟元，李昱，陳九輝，郭彪，李順成．2009．汶川MS8．0地震：地殼上地幔S波速度結(jié)構(gòu)的初步研究［J］．地球物理學(xué)報(bào)，52（2）：309-319．

孟令媛，史保平．2011．應(yīng)用動(dòng)態(tài)復(fù)合震源模型模擬汶川MW7．9地震強(qiáng)地面運(yùn)動(dòng)［J］．地球物理學(xué)報(bào)，54（4）：1010-1027．

孟令媛，史保平．2012．2011年新西蘭MW6．1地震震源過(guò)程及強(qiáng)地面運(yùn)動(dòng)特征初步分析［J］．地球物理學(xué)報(bào)，55（5）：1601-1612．

王衛(wèi)民，郝金來(lái)，姚振興．2013．2013年4月20日四川蘆山地震震源破裂過(guò)程反演初步結(jié)果［J］．地球物理學(xué)報(bào)，56（4）：1412-1417．

張勇，許力生，陳運(yùn)泰．2013．蘆山4·20地震破裂過(guò)程及其致災(zāi)特征初步分析［J］．地球物理學(xué)報(bào)，56（4）：1408-1411．

Andrews D J．1986．Objective determination of source parameters and similarity of earthquakes of different size［C］∥Das S，Boatwright J，Scholz C H eds．EarthquakeSourceMechanics．Washington：AGU：259-267．

Brune J N．1970．Tectonic stress and spectra of seismic shear waves from earthquakes［J］．JGeophysRes，75（26）：4997-5009．Brune J N．1971．Correction of"tectonic stress and spectra of seismic shear waves from earthquakes"［J］．JGeophys Res，76（20）：5002．

Brune J N．1976．The physics of earthquake strong motion［C］∥Lomnitz C，Rosenblueth E eds．SeismicRiskandEngineeringDecisions．New York：Elsevier Sci Publ Co：141-177．

Frankel A．1991．High-frequency spectral falloff for earthquakes，fractal dimension of strength on faults［J］．JGeophys Res，96（B4）：6291-6302．

Hanks T C，Bakun W H．2002．A bilinear source-scaling model forM-logAobservations of continental earthquakes［J］．BullSeismolSocAm，92（5）：1841-1846．

Keilis-Borok V I．1959．On estimation of the displacement in an earthquake source and source dimension［J］．Annalidi Geofisica，12（2）：205-214．

Knopoff L．1958．Energy release in earthquake［J］．GeophysJInter，1（1）：44-52．

Orowan E．1960．Mechanism of seismic faulting in rock deformation：A symposium［J］．GeolSocAmMem，79：323-345．

Savage J C，Wood M D．1971．The relation between apparent stress and stress drop［J］．BullSeismolSocAm，61（5）：1381-1388．

Smith K D，Brune J N，Priestly K F．1991．The seismic spectrum，radiated energy，and the Savage and Wood inequality for complex earthquake［J］．Tectonophysics，188（3／4）：303-320．

Starr A T．1928．Slip in a crystal and rupture in a solid due to shear［J］．ProcCambPhilSoc，24（4）：489-500．

Wald D J，Quitoriano V，Heaton T H，Thomas H，Kanamori H．1999．Relationships between peak ground acceleration，peak ground velocity，and modified Mercalli intensity in California［J］．EarthquakeSpectra，15（3）：557-564．

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

Zeng Y，Anderson J G．1996．A composite source modeling of the 1994Northridge earthquake using genetic algorithm［J］．BullSeismolSocAm，86（1B）：71-83．

Zú?iga F R．1993．Fractional overshoot and partial stress drop．Which one？［J］．BullSeismolSocAm，83（3）：939-944．