藥曉東 章文波

（中國北京100049中國科學(xué)院大學(xué)地球科學(xué)學(xué)院）

?

2013年四川蘆山MS7.0地震強(qiáng)地面運(yùn)動模擬

（中國北京100049中國科學(xué)院大學(xué)地球科學(xué)學(xué)院）

運(yùn)用經(jīng)驗(yàn)格林函數(shù)法模擬2013年4月20日蘆山MS7.0地震的近場強(qiáng)地面運(yùn)動． 在擬合過程中， 首先參考前人遠(yuǎn)場反演結(jié)果給出的滑動量分布特征和主震波形的包絡(luò)線特征， 確定強(qiáng)震動生成區(qū)的大致范圍和數(shù)量； 然后利用Somerville等提出的地震矩與凹凸體面積的經(jīng)驗(yàn)關(guān)系式確定強(qiáng)震動生成區(qū)細(xì)小劃分的初值， 繼而利用遺傳優(yōu)化算法確定以上二者的最優(yōu)值及其它震源參數(shù)． 將數(shù)值模擬波形與實(shí)際地震觀測記錄在時間域和頻率域分別進(jìn)行比較， 結(jié)果顯示， 在所選取的30個觀測臺站中， 多數(shù)臺站的數(shù)值模擬結(jié)果與實(shí)際觀測結(jié)果符合得很好， 特別是大于1 Hz的高頻部分． 斷層面上有兩個強(qiáng)震動生成區(qū)， 其位置與前人反演的滑動量集中分布區(qū)相一致， 而且強(qiáng)震動生成區(qū)規(guī)模比Somerville等獲得的標(biāo)度率估計值要?。?/p>

強(qiáng)地面運(yùn)動 蘆山MS7.0地震 經(jīng)驗(yàn)格林函數(shù)法 震源參數(shù) 波形對比

引言

2013年4月20日在我國四川省雅安市蘆山縣境內(nèi)發(fā)生了MS7.0地震， 這是繼2008年汶川大地震后發(fā)生在龍門山斷裂帶的又一次顯著的災(zāi)害性事件． 地震發(fā)生后， 我國地震工作者先后根據(jù)現(xiàn)場應(yīng)急科考、 GPS觀測數(shù)據(jù)、 遠(yuǎn)場地震數(shù)據(jù)和近震強(qiáng)震動數(shù)據(jù)發(fā)布了此次地震的發(fā)震構(gòu)造模型、 震源機(jī)制解、 地震震源破裂過程、 精定位及其與汶川地震的關(guān)系等研究結(jié)果（陳運(yùn)泰等， 2013； 王衛(wèi)民等， 2013； 武艷強(qiáng)等， 2013； 謝祖軍等， 2013； 徐錫偉等， 2013； 張勇等， 2013； 金明培等， 2014）． 此外， 孟令媛等（2014）利用隨機(jī)震源模型模擬了蘆山地震的烈度分布， 張冬麗等（2013）使用有限差分法模擬了其強(qiáng)地面運(yùn)動特征． 但是這些工作都沒有直接模擬近場強(qiáng)地面運(yùn)動的具體波形和頻譜， 而波形峰值、 持時及頻譜特征是地震動的三要素．

近年來發(fā)生的破壞性大地震, 例如我國1999年臺灣集集地震、 2008年汶川地震， 日本1995年神戶地震、 2003年十勝（Tokachi）地震、 2011年東北（Tohoku）地震， 1999年土耳其Kocaeli地震, 2010年海地地震以及2010年智利大地震， 都展示給我們： 大地震產(chǎn)生的強(qiáng)地面運(yùn)動是工程建筑物動力災(zāi)變的最主要原因， 是造成各類建筑物損毀甚至倒塌的直接原因， 是人員傷亡和各類次生災(zāi)害發(fā)生的主要誘因． 目前地震學(xué)和地震工程學(xué)研究人員已經(jīng)認(rèn)識到， 近場強(qiáng)地面運(yùn)動預(yù)測對于減輕未來大地震所造成的災(zāi)害、 重大工程的抗震設(shè)計以及地震危險性分析均有極其重要的作用， 所以如何有效地預(yù)測破壞性大地震的強(qiáng)地面運(yùn)動成為地震學(xué)和地震工程學(xué)領(lǐng)域中最為重要的課題之一．

目前預(yù)測未來大地震強(qiáng)地面運(yùn)動的方法主要包括： ① 確定性的數(shù)值計算方法， 即利用計算機(jī)通過數(shù)值計算來實(shí)現(xiàn)仿真模擬， 包括運(yùn)動學(xué)模擬和動力學(xué)模擬（Zhangetal, 2010； Ruiz, Madariaga, 2013）； ② 隨機(jī)方法， 即基于地震點(diǎn)源ω2模型（Brune, 1970, 1971）， 利用隨機(jī)震源理論來模擬地震的高頻地面地震動， 或者采用改進(jìn)的隨機(jī)有限斷層方法（Boore, 1983； Motazedian, Atkinson, 2005）； ③ 理論方法， 即以位移表示定理為基礎(chǔ)， 以理論格林函數(shù)計算為核心的方法（Eringen, Suhubi, 1975； Aki, Richards， 1980）． 雖然以上3種方法中理論方法具有嚴(yán)格的數(shù)學(xué)物理意義， 但是目前只能計算水平成層的理論格林函數(shù)， 因此運(yùn)用理論格林函數(shù)計算理論地震圖還存在很大的局限． 這3種方法均有各自不足： 數(shù)值方法對計算資源的消耗比較大， 費(fèi)時較多； 隨機(jī)方法缺乏明確的物理意義； 理論方法在非水平成層和非均勻介質(zhì)的計算中沒有解析解， 如果借助數(shù)值計算則耗時較多． 為了克服這些不足， Hartzell（1978） 提出一種經(jīng)驗(yàn)格林函數(shù)法----選取主震的小震（余震或前震）記錄作為格林函數(shù)， 取代理論格林函數(shù)的計算， 通過一定的破裂方式進(jìn)行疊加獲得主震的地面運(yùn)動記錄． 由于小震記錄已經(jīng)包含了震源、 傳播路徑和局部場地效應(yīng)的影響， 因此避免了計算理論格林函數(shù)的復(fù)雜性． 針對這一方法， 研究人員進(jìn)行了深入研究并加以改進(jìn)， 其中最重要的包括： Kanamori（1979）引入距離和輻射花樣校正； Irikura（1983）根據(jù)大小地震的相似性， 將小震的選取轉(zhuǎn)化為如何對小震進(jìn)行較正， 并運(yùn)用到實(shí)際大地震的強(qiáng)震動實(shí)例分析中， 由此該方法得到了廣泛應(yīng)用（Joyner, Boore, 1986; 羅奇峰, 1989; Danetal, 1990; 張有兵， 章文波， 2010； 藥曉東等， 2015）． 本文利用這一方法對2013年蘆山MS7.0地震的強(qiáng)地面運(yùn)動進(jìn)行模擬， 旨在為未來大地震的地震動預(yù)測和該地區(qū)的工程抗震設(shè)計提供可靠依據(jù)．

1 經(jīng)驗(yàn)格林函數(shù)法原理

利用理論方法計算地震圖的核心是計算格林函數(shù)， 雖然目前格林函數(shù)的計算取得了一些進(jìn)展（Kennett, 1980; Bouchon, 1981； Chen, 1995, 1996）， 但真實(shí)地球的復(fù)雜性， 使得三維復(fù)雜構(gòu)造的格林函數(shù)一直很難確定． 為了克服這一難題， Hartzell（1978）提出了利用大地震后余震記錄作為一種經(jīng)驗(yàn)格林函數(shù)來模擬大震記錄的方法． 其基本思路是： 一個小地震記錄（前震或者余震）， 由于已經(jīng)包含了震源、 傳播路徑和局部場地信息， 所以近似視為一對“源點(diǎn)-接收點(diǎn)”之間的格林函數(shù)， 通過這種格林函數(shù)與變化的震源空間時間函數(shù)作卷積， 可以合成大地震所引起的地面運(yùn)動. 這種合成方法是半經(jīng)驗(yàn)性的， 所以被稱為“經(jīng)驗(yàn)格林函數(shù)法”． 該方法的優(yōu)勢在于格林函數(shù)取自小震記錄， 自然包括了真實(shí)地球介質(zhì)從震源、 傳播路徑到接受點(diǎn)的場地條件影響的主要信息， 從而避免了理論格林函數(shù)的計算， 減小了計算量， 加快了計算速度．

經(jīng)驗(yàn)格林函數(shù)法最初的假定是一個臺站記錄到覆蓋整個破裂面的可視為點(diǎn)源的小震記錄， 則可以通過疊加模擬出該臺站的大震記錄（Hartzell, 1978）． 然而實(shí)際情況是臺站上的記錄遠(yuǎn)遠(yuǎn)不足， 所以該方法難以實(shí)現(xiàn)． 為了使該方法在實(shí)踐中得到更廣泛的應(yīng)用， Kanamori和Irikura先后作出了重要的改進(jìn)． Kanamori（1979）考慮到目標(biāo)地震具有一定斷層規(guī)模且由許多小震源組成， 所以即使個別小震記錄可以選作斷層面上的點(diǎn)源記錄， 也要考慮到位置因素， 因此必須進(jìn)行距離和輻射圖案校正． 此外， 為了使小震記錄作為格林函數(shù)， 要求其拐角頻率高于地震動的目標(biāo)最高頻率， 即要求所選擇的小震破裂尺度足夠小， 而且記錄的信噪比足夠高． 由于小震破裂尺度足夠小和記錄信噪比足夠高之間存在著固有矛盾， 所以此條件不容易滿足． Irikura（1983）依據(jù)對大小地震震源參數(shù)的統(tǒng)計研究結(jié)果， 提出了大小地震滿足相似律的假設(shè)， 并將大地震中的子源和所利用的小震的位錯函數(shù)均視為斜坡函數(shù)， 則大地震產(chǎn)生的地面運(yùn)動可由對子源在斷層上的二重求和轉(zhuǎn)換為對小震在斷層上和時間域內(nèi)的三重求和． 基于這些假設(shè)， 該方法放寬了對小震記錄的限制； 而且把問題從選擇什么樣的小震作為格林函數(shù)轉(zhuǎn)換為如何合理地對小震記錄進(jìn)行修正以得到大震中的子源貢獻(xiàn)， 這使得這一方法得到了廣泛的應(yīng)用（Joyner, Boore, 1986; 羅奇峰, 1989; Danetal, 1990; 張有兵， 章文波， 2010； 藥曉東等， 2015）．

本文主要采用Irikura（1983）所改進(jìn)的方法， 并針對當(dāng)所選取的“小震”較大和強(qiáng)震動生成區(qū)不唯一時， 如何劃分主斷層上的子斷層作了改進(jìn)． 圖1a為經(jīng)驗(yàn)格林函數(shù)法示意圖． 該方法基本原理的數(shù)學(xué)表達(dá)為

圖1 經(jīng)驗(yàn)格林函數(shù)法示意圖

（1）

（2）

（3）

式中：U（t）和u（t）分別為主震記錄和小震記錄， 這里將u（t）近似看作格林函數(shù)來擬合U（t）；r，r0，rij和ξij分別表示小震到臺站的距離、 主震初始破裂點(diǎn)到臺站的距離、 主震斷層上子斷層（i，j）到臺站的距離， 以及主震初始破裂點(diǎn)到子斷層（i，j）的距離；β為剪切波速度；VR為破裂傳播速度；N為主震斷層劃分為細(xì)小子斷層的特征數(shù)， 通常主震斷層劃分為小斷層的數(shù)目用N×N表示；C為地震尺寸和應(yīng)力降的參數(shù)， 與N聯(lián)合在一起表示大小地震之間規(guī)模與應(yīng)力降之比； 由于主震和小震在震源時間函數(shù)和頻率上存在差異， 所以在模擬時需要對小震記錄的頻譜進(jìn)行調(diào)整， 式中F（t）是一個濾波窗函數(shù)， 通過與小震記錄作卷積起到調(diào)整小震頻譜的作用， 其形狀參見圖1b；n′為對頻率上限進(jìn)行調(diào)整的參數(shù)．

從經(jīng)驗(yàn)格林函數(shù)法的原理和數(shù)學(xué)表達(dá)中可以看出， 該方法的關(guān)鍵在于經(jīng)驗(yàn)格林函數(shù)小震的選?。?如果小震記錄能夠在震源過程、 傳播路徑和場地條件等3方面與主震一致， 那么通過一定的疊加則可以很好地合成主震記錄． 但在實(shí)際應(yīng)用中， 對于同一個臺站的記錄， 通常只有場地條件最容易滿足， 而傳播路徑和震源信息則無法實(shí)現(xiàn)一致， 所以須綜合考慮小震與主震的差異， 通過距離和頻率校正來生成可靠的記錄（張有兵， 章文波， 2010； 藥曉東等， 2015）．

2 蘆山MS7.0地震分析

2.1 經(jīng)驗(yàn)格林函數(shù)小震的選取

如前文所述， 經(jīng)驗(yàn)格林函數(shù)法的優(yōu)勢在于包含了震源過程、 傳播路徑和場地條件上的主要信息， 因此擬合的關(guān)鍵在于小震的選取是否恰當(dāng)． 因?yàn)樾≌鹋c主震的空間位置和震源機(jī)制是否相似決定了擬合的優(yōu)劣， 所以我們將震源位置和震源機(jī)制作為選擇小震的標(biāo)準(zhǔn)． 此外擬合臺站的數(shù)量和覆蓋范圍也是我們考慮的因素， 擬合臺站數(shù)量過少或者不能覆蓋主震斷層的不同方位也是需要避免的． 為了滿足擬合臺站數(shù)量和記錄信噪比， 本文只考慮M≥3.0的小震．

蘆山地震后， 不同研究小組利用全球測震資料發(fā)布了此次地震的震源機(jī)制． 哈佛大學(xué)全球質(zhì)心矩張量目錄（Dziewonskietal, 1981; Ekstr?metal, 2012）給出的震源機(jī)制為： 斷層走向212°， 傾角42°， 滑動角100°， 深度21.9 km， 地震矩1.02×1019N·m，MW6.6 （GCMT， 2013）； 美國地質(zhì)調(diào)查局地震情報中心給出的震源機(jī)制為： 斷層走向198°， 傾角33°， 滑動角71°， 深度12 km， 地震矩1.0×1019N·m，MW6.6 （USGS， 2013）； 中國地震局地球物理研究所陳運(yùn)泰院士研究組給出的震源機(jī)制為： 斷層走向220°， 傾角35°， 滑動角95°， 深度12 km， 地震矩1.6×1019N·m，MW6.7 （中國地震臺網(wǎng)中心， 2013）． 不同機(jī)構(gòu)給出的震源機(jī)制結(jié)果均反映了震源的逆沖特征， 本文將主要參考國內(nèi)機(jī)構(gòu)給出的結(jié)果， 即斷層走向220°， 傾角35°， 滑動角95°—101°， 深度12—13 km， 地震矩（1.4—1.6）×1019N·m． 國家強(qiáng)震動臺網(wǎng)中心獲得了此次地震的主、 余震記錄， 其中主震記錄363條， 余震記錄3300多條， 余震震級范圍為MS2.3—5.4， 所有地震記錄均為三分向加速度記錄， 采樣頻率均為200 Hz．

表1 主震和作為經(jīng)驗(yàn)格林函數(shù)的余震的震源位置和震源機(jī)制Table 1 Focal mechanisms and locations of the main shock and the aftershock used as EGF in this study

表2 強(qiáng)震動臺站位置信息Table 2 Locations of the 30 strong motion stations used in this study

根據(jù)以上信息， 并結(jié)合哈佛大學(xué)全球質(zhì)心矩張量目錄以及呂堅(jiān)等（2013）給出的余震震源機(jī)制和定位信息， 本文篩選出一次MW4.8余震作為經(jīng)驗(yàn)格林函數(shù)小震. 主震和所選取余震的震源參數(shù)見表1， 記錄地震數(shù)據(jù)的臺站信息見表2， 主震、 余震和臺站的空間分布如圖2所示． 此次模擬涉及的臺站為30個． 盡管另外兩個臺站（051LSL樂山凌云臺和051MCL沐川利店臺）也獲取了主、 余震記錄， 但由于這兩個臺站的主震記錄異常而被剔除． 圖3給出了典型的余震記錄加速度時程及其傅里葉譜（051BXD和051CDZ臺站）． 我們按以下步驟對數(shù)據(jù)進(jìn)行統(tǒng)一處理： 首先進(jìn)行基線校正， 然后利用凱澤（Kaiser）窗帶通濾波器對記錄進(jìn)行濾波， 帶通頻率范圍為0.2—20 Hz．

圖2 主震、 余震和臺站空間分布圖

2.2 模型參數(shù)估計

蘆山地震發(fā)生后， 研究人員利用遠(yuǎn)場地震數(shù)據(jù)（王衛(wèi)民等， 2013； 張勇等， 2013； 趙翠萍等， 2013）以及GPS數(shù)據(jù)與近場地震數(shù)據(jù)聯(lián)合（金明培等， 2014）反演了汶川地震的斷層面滑動量．

王衛(wèi)民等（2013）的斷層模型為： 主斷層走向205°， 傾角38.5°， 斷層規(guī)模66 km×35 km． 反演結(jié)果顯示此次地震為逆沖事件， 最大滑動量為1.59 m， 最大滑動量所在的凹凸體基本位于斷層中部， 在震源深度偏上位置， 在凹凸體周圍同時還分布著3個較小的滑動量集中區(qū)（圖4a）． 張勇等（2013）的斷層模型為： 主斷層走向219°， 傾角33°， 斷層規(guī)模60 km×45 km． 反演結(jié)果顯示這是一個以逆沖為主、 但兼具比汶川地震還小的右旋走滑分量的地震， 最大滑動量為1.3 m， 最大滑動量所在的凹凸體基本位于斷層中部， 包括了震源所在位置， 在凹凸體西南方、 斷層面下緣同時還有另一個較小的滑動量集中區(qū)（圖4b）． 趙翠萍等（2013）的斷層模型為： 主斷層走向220°， 傾角35°， 斷層規(guī)模100 km×45 km． 反演結(jié)果顯示此次地震破裂自起始破裂點(diǎn)向地表快速發(fā)展, 呈雙側(cè)破裂的逆沖事件， 凹凸體在起始破裂點(diǎn)之上40 km×30 km 區(qū)域, 最大滑動量為1.8 m， 在震中地表的北東方向有一個滑動量為0.4 m左右的區(qū)域（圖4c）． 金明培等（2014）的斷層模型為： 主斷層走向212°， 傾角35°—54°， 斷層規(guī)模60 km×39 km． 反演結(jié)果顯示此次地震以逆沖為主， 最大滑動量為1.1 m， 最大滑動量所在的凹凸體基本位于斷層中部， 包含了震源所在位置， 在凹凸體周圍同時還分布著3個較小的滑動量集中區(qū)（圖4d）． 無論是遠(yuǎn)震臺站反演， 還是GPS與近場強(qiáng)震聯(lián)合反演， 其結(jié)果都有共同的特征， 即在斷層面中心位置附近有唯一的凹凸體， 其最大滑動量為1.1—1.8 m．

不同的反演結(jié)果均顯示蘆山地震的震源機(jī)制比較一致， 同時反演給出的震源過程也比較簡單． 此外由記錄到的加速度波形可以看出， 大多數(shù)臺站的記錄只體現(xiàn)出一個波包的整體特征， 這也說明蘆山地震的發(fā)震過程比較簡單． 但對所有的主震記錄進(jìn)行觀察和分析時， 發(fā)現(xiàn)部分臺站的記錄體現(xiàn)出兩個明顯的波包特征． 以051LSF臺站記錄為例（圖5）， 從包絡(luò)線上可以看到EW和UD分量存在間隔為7—8 s的兩個明顯波包， 在EW和UD分量對應(yīng)的第二個波包處， NS分量也有個明顯的波包出現(xiàn)， 但在EW和UD分量的第一個波包位置， NS分量分裂為兩個較小的波包． 因此可以初步判斷蘆山地震的強(qiáng)震動生成不是一次單一事件， 而是由至少兩次事件， 即至少兩個強(qiáng)震動生成區(qū)的發(fā)震所引起的． 由前文遠(yuǎn)場反演的滑動量分布可以看出這兩個強(qiáng)震動生成區(qū)的震源機(jī)制相互一致， 因此考慮選擇一個小地震記錄作為經(jīng)驗(yàn)格林函數(shù)． 已有的反演結(jié)果也顯示， 在蘆山主震的斷層面上有一個主要的凹凸體和1—3個較小的局部滑動顯著區(qū)（王衛(wèi)民等， 2013； 張勇等， 2013； 趙翠萍等， 2013； 金明培等， 2014）， 因此本文構(gòu)建的基本模型為兩個強(qiáng)震動生成區(qū)， 其中較大的生成區(qū)位于震源附近， 較小的在其周圍． 對于是否需要增加強(qiáng)震動生成區(qū)的數(shù)量則應(yīng)根據(jù)實(shí)際模擬來確定， 同時參考前文的反演斷層模型， 本文將事件發(fā)生的斷層設(shè)定為66 km×35 km． 由于蘆山地震震源機(jī)制變化較小， 所以模型中僅選用了同一個經(jīng)驗(yàn)格林函數(shù)余震． 模型中還涉及到強(qiáng)震動生成區(qū)的位置， 本文主要參考王衛(wèi)民等（2013）的反演結(jié)果， 同時也參考了其它反演結(jié)果， 最終的確切位置由遺傳算法優(yōu)化搜索獲取．

圖3 051BXD（a）和051CDZ（b）臺站記錄到的余震加速度時程（上）及其傅里葉譜與噪聲譜（下）

圖4 蘆山地震斷層面上的滑動量分布

圖5 051LSF臺站記錄的蘆山主震加速度時程. 圖中數(shù)字為記錄到的加速度峰值

確定強(qiáng)震動生成區(qū)的可能位置后， 還要確定強(qiáng)震動生成區(qū)劃分成子塊的數(shù)量N． 以往主要根據(jù)地震的相似性， 包括幾何形狀、 應(yīng)力降、 平均位移和上升時間等， 推導(dǎo)出N的確定方法以及改進(jìn)的震源譜比法（Irikura, 1983; Irikura, Kamae, 1994）． 但是由于這些方法在強(qiáng)震動生成區(qū)不唯一， 無法分配不同強(qiáng)震動生成區(qū)的大小， 因此在蘆山地震模擬的實(shí)例中， 我們首先借鑒凹凸體與地震矩之間的標(biāo)度關(guān)系進(jìn)行初值估計， 然后再利用優(yōu)化算法最終確定N值．

Somerville等（1999）根據(jù)15個MW5.6—7.2地震的研究結(jié)果給出了地震矩與凹凸體總面積的經(jīng)驗(yàn)關(guān)系為

（4）

式中：M0為地震矩， 單位為N·m；Aa為凹凸體總面積， 單位為km2． 根據(jù)哈佛大學(xué)全球質(zhì)心矩張量目錄給出的結(jié)果： 蘆山地震主震標(biāo)量地震矩為1.02×1019N·m， 小震的標(biāo)量地震矩為2.15×1016N·m， 兩者凹凸體總面積之比為60.8． 據(jù)此假設(shè)主震與小震的強(qiáng)震動生成區(qū)的比值為61， 此即為兩個強(qiáng)震動生成區(qū)細(xì)小劃分總和的初始估計值． 由于不同研究小組給出的地震矩存在一些差別， 其中較大的一個是中國地震局地球物理研究所陳運(yùn)泰院士研究組給出的1.6×1019N·m （中國地震臺網(wǎng)中心， 2013）， 求出的相應(yīng)比值為82． 為此， 我們將主震與小震的凹凸體面積之比的初始范圍設(shè)定為61—82， 具體劃分則要借助遺傳算法進(jìn)行優(yōu)化搜索． 這里之所以選擇遺傳算法是因?yàn)楸疚牡哪M問題是多極值問題， 智能優(yōu)化算法是相對更合理的解決策略． 根據(jù)以往將強(qiáng)震動生成區(qū)劃分為N×N塊的方法， 我們也將每個強(qiáng)震動生成區(qū)的子網(wǎng)格在走向方向和傾向方向按相同數(shù)量進(jìn)行劃分．

在確定強(qiáng)震動生成區(qū)的大小時， 我們首先設(shè)定強(qiáng)震動生成區(qū)劃分的子斷層長寬相等， 然后采用遺傳算法（Carroll, 2001）在0.1—5 km范圍內(nèi)進(jìn)行搜索． 在確定斷層破裂速度時參考相關(guān)研究， 將剪切波速度設(shè)定為3.5 km/s， 同時將斷層破裂速度設(shè)定為2.1—3.6 km/s．這是因?yàn)槟壳把芯拷Y(jié)果表明， 破裂過程的平均速度一般為剪切波速的60%—90%， 但對于個別大地震的破裂過程， 由于其局部某個時段有可能發(fā)生超剪切破裂（張海明， 2006）， 因此放寬了搜索范圍． 對于其它震源參數(shù)， 依然借助遺傳算法搜索得到．

在利用遺傳算法優(yōu)化搜索參數(shù)時， 衡量參數(shù)最優(yōu)的標(biāo)準(zhǔn)是使模擬波形與真實(shí)記錄波形之間的適配函數(shù)達(dá)到最?。?本文使用的適配函數(shù)是對加速度、 速度和位移波形的三分量全部進(jìn)行比較， 且采用等權(quán)重相加． 由于時間域波形比較復(fù)雜， 所以我們采用比較包絡(luò)線的方式， 具體形式為

Emis=

（5）

式中，eobs-acc，eobs-vel和eobs-disp分別表示觀測的加速度、 速度和位移的包絡(luò)線；esyn-acc，esyn-vel和esyn-disp分別表示擬合的加速度、 速度和位移的包絡(luò)線．

3 結(jié)果與討論

通過模型參數(shù)的估計和遺傳算法優(yōu)化搜索確定了最優(yōu)模型， 模型參數(shù)見表3， 其中N1=8，N2=4， 子斷層劃分的總和為80， 介于初始估計值61—82之間． 強(qiáng)震動生成區(qū)1和強(qiáng)震動生成區(qū)2的面積分別為64 km2和16 km2， 具體位置和面積大小見圖6， 強(qiáng)震動生成區(qū)的總面積為80 km2， 比用式（4）估計的109 km2， 即Somerville等（1999）獲得的標(biāo)度率的估計值要小． 繼Somerville等（1999）研究標(biāo)度率后， Irikura和Miyake（2011）也總結(jié)了類似的標(biāo)度率， 但其總結(jié)的主要是破裂面積與地震矩之間的標(biāo)度關(guān)系． 我們借助Somerville等（1999）獲得的凹凸體與地震矩之間的標(biāo)度關(guān)系式（4）， 以及破裂面與地震矩之間的關(guān)系得到

（6）

表3 最優(yōu)模型參數(shù)和遺傳算法搜索范圍Table 3 The optimal parameters of determined model and their search scopes for genetic algorithm

圖6 基于王衛(wèi)民等（2013）的反演結(jié)果構(gòu)建的強(qiáng)震動震源模型

式中：M0為地震矩， 單位為N·m；A為破裂面積， 單位為km2． 這樣可以估計出破裂面積與凹凸體面積之比為4.46， 借助該比值， 我們將Irikura和Miyake（2011）總結(jié)的破裂面積與地震矩之間的標(biāo)度關(guān)系轉(zhuǎn)化為凹凸體總面積與地震矩之間的標(biāo)度關(guān)系, 即

（7）

（8）

Aa=1.19×10-18M0，M0≥7.5×1020,

（9）

式中：M0為地震矩， 單位為N·m；Aa為凹凸體總面積， 單位為km2． 該標(biāo)度率為三段形式， 見圖7a的黑色實(shí)線． 可以看出， 式（4）與式（7）表達(dá)相同， 意味著中等震級地震在Somerville等（1999）和Irikura和Miyake（2011）中服從相同的標(biāo)度關(guān)系， 但從式（8）和式（9）可以看出， Irikura和Miyake（2011）對大地震和巨大地震所總結(jié)的標(biāo)度關(guān)系與Somerville等（1999）出現(xiàn)了差異， 這是因?yàn)樵诳偨Y(jié)標(biāo)度關(guān)系時Irikura和Miyake（2011）加入了更多的地震數(shù)據(jù)， 尤其是大地震和巨大地震的數(shù)據(jù)． 雖然此次模擬比式（4）中對應(yīng)的面積要小， 但式（8）同樣出現(xiàn)了比式（4）面積小的趨勢． 因此此次模擬結(jié)果佐證了式（8）的趨勢， 也為大地震的標(biāo)度率總結(jié)提供了一則實(shí)例依據(jù)．

圖7 強(qiáng)震動生成區(qū)面積（a）和上升時間（b）與地震矩的標(biāo)度關(guān)系

由圖6可以看出， 最優(yōu)模型的強(qiáng)震動生成區(qū)1的初始破裂點(diǎn)幾乎位于正中心， 然后向兩側(cè)破裂， 這也驗(yàn)證了蘆山地震由震源發(fā)震， 向雙側(cè)破裂的特征； 同時與應(yīng)力降有關(guān)的參數(shù)C1和C2分別為1.6和0.7， 反映了強(qiáng)震動生成區(qū)1的應(yīng)力降比強(qiáng)震動生成區(qū)2的要大； 而兩個強(qiáng)震動生成區(qū)的上升時間均為1.0 s， 與Somerville等（1999）獲得的標(biāo)度率的估計值0.95 s非常接近（圖7b）．

從圖6可以看出， 強(qiáng)震動生成區(qū)1與凹凸體位置一致， 體現(xiàn)了高頻能量與低頻能量的主要生成區(qū)域是相似的； 但強(qiáng)震動生成區(qū)2卻與其它反演結(jié)果存在差別， 這源于本文的模擬頻率范圍與確定性方法的頻率范圍不同， 本文的頻率范圍為1—20 Hz， 而確定性方法由于其數(shù)值計算的頻率受計算能力和耗時的影響， 往往低于1 Hz， 所以出現(xiàn)兩者不同的可能性非常大．

圖8給出了加速度波形的模擬結(jié)果． 可以看出， 051BXD臺站的EW分量達(dá)到了1005.4 cm/s2， 是我國內(nèi)陸地區(qū)觀測到的首個超過1g（1g=9.8×102m/s2）的記錄． 該臺站無論是時間域還是頻率域均擬合得很好． 水平加速度峰值超過100 cm/s2的臺站共有12個， 其中10個臺站（051BXM， 051BXY， 051BXZ， 051HYY， 051LDL， 051LSF， 051TQL， 051YAD， 051YAL和051YAM）在時間域和頻率域均擬合得很好， 只有2個臺站（051PJD和051QLY）的估計值明顯偏?。?這種偏小估計值也出現(xiàn)在051CDZ和051DJZ臺站上. 這4個臺站剛好位于蘆山地震主震斷層的北東向， 而且這種明顯偏小的估計值并未出現(xiàn)在其它方位的臺站上． 結(jié)合哈佛大學(xué)質(zhì)心矩張量目錄給出的蘆山地震主、 余震的滑動角分別為100°和74°（表1）， 可以看出， 主震體現(xiàn)為逆沖為主兼微小右滑機(jī)制， 而余震體現(xiàn)為逆沖為主兼微小左滑機(jī)制. 這體現(xiàn)了余震震源機(jī)制對擬合的影響直接關(guān)系到加速度峰值和幅值譜的大?。?擬合結(jié)果中最差的是051QLY臺站的EW分量， 無論在時間域或頻率域均偏離較大． 其原因一方面由于震源機(jī)制的影響， 另一方面與該記錄的信噪比有關(guān)（圖9）． 通過051QLY臺站EW分量的信噪比與其它記錄的信噪比對比（如051YAM臺站的EW分量）可以看出， 051QLY臺站的EW分量在1 Hz左右的信噪比很低， 因此在擬合過程中會影響結(jié)果． 在分析過程中我們也考慮過在斷層靠近北東向的位置增加一個強(qiáng)震動生成區(qū)， 這樣可以彌補(bǔ)051PJD臺站估計值偏小的問題， 但增加一個強(qiáng)震動生成區(qū)會增大對其它臺站的貢獻(xiàn)， 導(dǎo)致其它臺站模擬結(jié)果偏大， 因此我們對兩個強(qiáng)震動生成區(qū)的模型并未作調(diào)整．

圖8 各臺站加速度時程（a）和加速度傅里葉譜（b）的觀測記錄（黑色）與最優(yōu)模型的擬合結(jié)果（灰色）對比（Ⅰ）

圖8 各臺站加速度時程（a）和加速度傅里葉譜（b）的觀測記錄（黑色）與最優(yōu)模型的擬合結(jié)果（灰色）對比（Ⅱ）

圖8 各臺站加速度時程（a）和加速度傅里葉譜（b）的觀測記錄（黑色）與最優(yōu)模型的擬合結(jié)果（灰色）對比（Ⅲ）

結(jié)合模擬過程中的計算， 我們發(fā)現(xiàn)擬合對小震的定位非常敏感， 而對主余震斷層的幾何長寬的敏感度較差． 這一點(diǎn)可能是由于主、 余震斷層的幾何形狀比較相似， 而且多數(shù)臺站的震中距要大于斷層幾何尺度的緣故．

4 結(jié)論

本文利用經(jīng)驗(yàn)格林函數(shù)法對蘆山地震的30個臺站進(jìn)行模擬研究， 選擇一個余震記錄作為經(jīng)驗(yàn)格林函數(shù)， 采用兩個強(qiáng)震生成區(qū)進(jìn)行模擬， 結(jié)果顯示在時間域和頻率域， 尤其是大于1 Hz的高頻部分， 多數(shù)臺站均擬合得很好． 對于加速度峰值超過100 cm/s2的臺站， 特別是051BXD臺站， 其擬合結(jié)果非常好， 體現(xiàn)出經(jīng)驗(yàn)格林函數(shù)法的優(yōu)勢所在， 也說明了經(jīng)驗(yàn)格林函數(shù)法的可靠性．

圖9 051QLY（a）和051YAM（b）臺站記錄到的余震加速度時程（上）及其傅里葉譜與噪聲譜（下）

小震震源機(jī)制對模擬結(jié)果的影響非常明顯， 其與目標(biāo)地震的震源機(jī)制越接近, 則模擬結(jié)果越好． 我們還發(fā)現(xiàn)擬合結(jié)果對小震的定位非常敏感， 而我們所掌握的小震資料是有限的， 因此關(guān)于小震震源機(jī)制和精確定位的相關(guān)研究， 對強(qiáng)震動模擬非常重要． 同時小震記錄的信噪比也是影響模擬的重要因素．

本文確定的強(qiáng)震動生成區(qū)的模型參考了遠(yuǎn)場反演的模型， 但最終結(jié)果與其又有一定的差別， 因此總結(jié)兩者之間的異同， 尋找兩者之間的關(guān)聯(lián)對震源研究和強(qiáng)震動模擬具有重要的意義．

雖然經(jīng)驗(yàn)格林函數(shù)法對強(qiáng)地面運(yùn)動的模擬非常有效， 但是對于預(yù)測未來大地震的強(qiáng)地面運(yùn)動， 還需要對大量的地震實(shí)例和模擬進(jìn)行總結(jié)， 對小震的震級、 定位、 震源機(jī)制以及主余震震源的差異等都需要進(jìn)一步研究． 同時地震標(biāo)度率的總結(jié)是未來災(zāi)害性地震動預(yù)測的基礎(chǔ)， 是特征化模型的重要組成部分， 因此地震標(biāo)度率尤其是大地震的標(biāo)度關(guān)系， 需要加入更多的地震實(shí)例， 進(jìn)行更深入的研究．

陳運(yùn)泰, 楊智嫻, 張勇, 劉超. 2013. 從汶川地震到蘆山地震[J]. 中國科學(xué): 地球科學(xué), 43（6）: 1064--1072.

Chen Y T, Yang Z X, Zhang Y, Liu C. 2013. From 2008 Wenchuan earthquake to 2013 Lushan earthquake[J].ScientiaSinicaTerrae, 43（6）: 1064--1072 （in Chinese）.

金明培, 汪榮江, 屠泓為. 2014. 蘆山7級地震的同震位移估計和震源滑動模型反演嘗試[J]. 地球物理學(xué)報, 57（1）: 129--137.

Jin M P, Wang R J, Tu H W. 2014. Slip model and co-seismic displacement field derived from near-source strong motion records of the LushanMS7.0 earthquake on 20 April 2013[J].ChineseJournalofGeophysics, 57（1）: 129--137 （in Chinese）.

羅奇峰. 1989. 近場加速度的半經(jīng)驗(yàn)合成[D]. 哈爾濱: 中國地震局工程力學(xué)研究所： 8--14.

Luo Q F. 1989.TheSemi-EmpiricalSynthesisofNear-FieldGroundMotions[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 8--14 （in Chinese）.

呂堅(jiān), 王曉山, 蘇金蓉, 潘林山, 李正, 尹利文, 曾新福, 鄧輝. 2013. 蘆山7.0級地震序列的震源位置與震源機(jī)制解特征[J]. 地球物理學(xué)報, 56（5）: 1753--1763.

Lü J, Wang X S, Su J R, Pan L S, Li Z, Yin L W, Zeng X F, Deng H. 2013. Hypocenter location and source mechanism of theMS7.0 Lushan earthquake sequence[J].ChineseJournalofGeophysics, 56（5）: 1753--1763 （in Chinese）.

孟令媛, 周龍泉, 劉杰. 2014. 2013年四川蘆山MS7.0地震近斷層強(qiáng)地面運(yùn)動模擬及烈度分布估計[J]. 地球物理學(xué)報, 57（2）: 441--448.

Meng L Y, Zhou L Q, Liu J. 2014. Estimation of the near-fault strong ground motion and intensity distribution of the 2013 Lushan, Sichuan,MS7.0 earthquake[J].ChineseJournalofGeophysics, 57（2）: 441--448 （in Chinese）.

王衛(wèi)民, 郝金來, 姚振興. 2013. 2013年4月20日四川蘆山地震震源破裂過程反演初步結(jié)果[J]. 地球物理學(xué)報, 56（4）: 1412--1417.

Wang W M, Hao J L, Yao Z X. 2013. Preliminary result for rupture process of Apr. 20, 2013, Lushan earthquake, Sichuan, China[J].ChineseJournalofGeophysics, 56（4）: 1412--1417 （in Chinese）.

武艷強(qiáng), 江在森, 王敏, 車時, 廖華, 李強(qiáng), 李鵬, 楊永林, 向和平, 邵志剛, 王武星, 魏文薪, 劉曉霞. 2013. GPS監(jiān)測的蘆山7.0級地震前應(yīng)變積累及同震位移場初步結(jié)果[J]. 科學(xué)通報, 58（20）: 1910--1916.

Wu Y Q, Jiang Z S, Wang M, Che S, Liao H, Li Q, Li P, Yang Y L, Xiang H P, Shao Z G, Wang W X, Wei W X, Liu X X. 2013. Preliminary results of the co-seismic displacement and pre-seismic strain accumulation of the LushanMS7.0 earthquake reflected by the GPS surveying[J].ChineseScienceBulletin, 58（20）: 1910--1916 （in Chinese）.

謝祖軍, 金筆凱, 鄭勇, 葛粲, 熊熊, 熊誠, 許厚澤. 2013. 近遠(yuǎn)震波形反演2013年蘆山地震震源參數(shù)[J]. 中國科學(xué): 地球科學(xué), 43（6）: 1010--1019.

Xie Z J, Jin B K, Zheng Y, Ge C, Xiong X, Xiong C, Hsu H T. 2013. Source parameters inversion of the 2013 Lushan earthquake by combining teleseismic waveforms and local seismograms[J].ScientiaSinicaTerrae, 43（6）: 1010--1019 （in Chinese）.

徐錫偉, 聞學(xué)澤, 韓竹軍, 陳桂華, 李傳友, 鄭文俊, 張世民, 任治坤, 許沖, 譚錫斌, 魏占玉, 王明明, 任俊杰, 何仲, 梁明劍. 2013. 四川蘆山7.0級強(qiáng)震: 一次典型的盲逆斷層型地震[J]. 科學(xué)通報, 58（20）: 1887--1893.

Xu X W, Wen X Z, Han Z J, Chen G H, Li C Y, Zheng W J, Zhang S M, Ren Z K, Xu C, Tan X B, Wei Z Y, Wang M M, Ren J J, He Z, Liang M J. 2013. LushanMS7.0 earthquake: A blind reserve-fault earthquake[J].ChineseScienceBulletin, 58（20）: 1887--1893 （in Chinese）.

藥曉東, 章文波, 于湘?zhèn)? 2015. 2008年汶川8.0級大地震近場強(qiáng)地面運(yùn)動的模擬[J]. 地球物理學(xué)報, 58（3）: 886--903.

Yao X D, Zhang W B, Yu X W. 2015. Simulation of near-field strong ground motion caused by the 2008MS8.0 Wenchuan earthquake[J].ChineseJournalofGeophysics, 58（3）: 886--903 （in Chinese）.

張冬麗, 黃蓓, 張獻(xiàn)兵, 徐錫偉, 鄭文俊. 2013. 基于有限斷層模型的蘆山“4·20”7.0級強(qiáng)烈地震強(qiáng)地面運(yùn)動特征[J]. 地震地質(zhì), 35（2）: 423--435.

Zhang D L, Huang B, Zhang X B, Xu X W, Zheng W J. 2013. Strong ground motion distribution and simulation based on finite fault model of Lushan 7.0 earthquake on April 20, 2013[J].SeismologyandGeology, 35（2）: 423--435 （in Chinese）.

張海明. 2006. 復(fù)雜斷層系統(tǒng)動力學(xué)破裂的理論研究和地表影響下的超剪切破裂問題的初步研究[D]. 北京: 北京大學(xué)地球科學(xué)與空間學(xué)院： 1--5.

Zhang H M. 2006.TheoreticalStudyonDynamicsRupturesonComplexFaultSystemsandPreliminaryStudyonSupershearRupturesUndertheEffectofGround[D]. Beijing: School of Earth and Space Sciences, Peking University: 1--5 （in Chinese）.

張勇, 許力生, 陳運(yùn)泰. 2013. 蘆山4·20地震破裂過程及其致災(zāi)特征初步分析[J]. 地球物理學(xué)報, 56（4）: 1408--1411.

Zhang Y, Xu L S, Chen Y T. 2013. Rupture process of the Lushan 4·20 earthquake and preliminary analysis on the disaster causing mechanism[J].ChineseJournalofGeophysics, 56（4）: 1408--1411 （in Chinese）.

張有兵, 章文波. 2010. 用經(jīng)驗(yàn)格林函數(shù)法模擬2008年日本巖手—宮城地震的強(qiáng)地面運(yùn)動[J]. 地震學(xué)報, 32（3）: 320--331.

Zhang Y B, Zhang W B. 2010. Strong ground motion simulation of the 2008 Iwate-Miyagi, Japan, earthquake using empirical Green’s function method[J].ActaSeismologicaSinica, 32（3）: 320--331 （in Chinese）.

趙翠萍, 周連慶, 陳章立. 2013. 2013年四川蘆山MS7.0級地震震源破裂過程及其構(gòu)造意義[J]. 科學(xué)通報, 58（20）: 1894--1900.

Zhao C P, Zhou L Q, Chen Z L. 2013. Source rupture process of LushanMS7.0 earthquake, Sichuan, China and its tectonic implication[J].ChineseScienceBulletin, 58（20）: 1894--1900 （in Chinese）.

中國地震臺網(wǎng)中心. 2013. 四川省雅安市蘆山縣7.0級地震[EB/OL]. [2015-06-10]. http:∥news.ceic.ac.cn/CC20130420080246.html.

China Earthquake Networks Center. 2013.MS7.0 earthquake occurred in Lushan, Ya’an city, Sichuan Province[EB/OL]. [2015-06-10]. http:∥news.ceic.ac.cn/CC20130420080246.html （in Chinese）.

Aki K, Richards P G. 1980.QuantitativeSeismology:TheoryandMethods[M]. San Francisco: W H Freeman and Company: 63--113.

Boore D M. 1983. Stochastic simulation of high-frequency ground motions based on seismological model of the radiated spectra[J].BullSeismolSocAm, 73（6A）: 1865--1894.

Bouchon M. 1981. A simple method to calculate Green’s functions for elastic layered media[J].BullSeismolSocAm, 74（4）: 959--971.

Brune J R. 1970. Tectonic stress and the spectra of seismic shear waves from earthquakes[J].JGeophysRes, 75（26）: 4997--5009.

Brune J R. 1971. Correction to “Tectonic stress and the spectra of seismic shear waves from earthquakes”[J].JGeophysRes, 76（20）: 5002.

Carroll D L. 2001. FORTRAN genetic algorithm （GA） driver[EB/OL]. [2015-06-10]. http:∥www.cuaerospace.com/Technology/GeneticAlgorithm/GADriverFreeVersion.aspx.

Chen X F. 1995. Seismogram synthesis for multi-layered media with irregular interfaces by global reflection/transmission matrices method: Ⅱ. Applications for 2D SH case[J].BullSeismolSocAm, 85（4）: 1094--1106.

Chen X F. 1996. Seismogram synthesis for multi-layered media with irregular interfaces by global reflection/transmission matrices method: Ⅲ. Theory for P-SV case[J].BullSeismolSocAm, 86（2）: 389--405.

Dan K, Watanabe T, Tanaka T, Sato R. 1990. Stability of earthquake ground motion synthesized by using different small-event records as empirical Green’s functions[J].BullSeismolSocAm, 80（6A）: 1433--1455.

Dziewonski A M, Chou T A, Woodhouse J H. 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[J].JGeophysRes, 86（B4）: 2825--2852.

Ekstr?m G, Nettles M, Dziewonski A M. 2012. The global CMT project 2004—2010: Centroid-moment tensors for 13017 earthquakes[J].PhysEarthPlanetInter, 200: 1--9.

Eringen A, Suhubi E. 1975.Elastodynamics，VolumeⅡ:LinerTheory[M]. New York: Academic Press, Inc: 717--834.

GCMT. 2013. Global CMT catalog[EB/OL]. [2015-06-10]. http:∥www.globalcmt.org/cgi-bin/globalcmt-cgi-bin/CMT4/form?itype=ymd&yr=2013&mo=4&day=20&otype=ymd&oyr=2013&omo=4&oday=21&jyr=1976&jday=1&ojyr=1976&ojday=1&nday=1&lmw=0&umw=10&lms=0&ums=10&lmb=0&umb=10&llat=-90&ulat=90&llon=-180&ulon=180&lhd=0&uhd=1000<s=-9999&uts=9999&lpe1=0&upe1=90&lpe2=0&upe2=90&list=0.

Hartzell S H. 1978. Earthquake aftershocks as Green’s functions[J].GeophysResLett, 5（1）: 1--4.

Irikura K. 1983. Semi-empirical estimation of strong ground motions during large earthquakes[J].BullDisasPrevResInst, 33（298）: 63--104.

Irikura K, Kamae K. 1994. Estimation of strong ground motion in broad-frequency band based on a seismic source scaling model and an empirical Green’s function technique[J].AnnGeophys, 37（6）: 25--47.

Irikura, K, Miyake H. 2011. Recipe for predicting strong ground motion from crustal earthquake scenarios[J].PureApplGeophys, 168（1/2）: 85--104.

Joyner W B, Boore D M. 1986. On simulating large earthquakes by Green’s function addition of smaller earthquakes[C]∥EarthquakeSourceMechanics. Washington: American Geophysical Union: 269--274.

Kanamori H. 1979. A semi-empirical approach to prediction of long-period ground motions from great earthquakes[J].BullSeismolSocAm, 69（6）: 1645--1670.

Kennett B L N. 1980.Seismic waves in stratified halfspace: Ⅱ. Theoretical seismograms[J].GeophysJInt, 61（1）: 1--10.

Motazedian D, Atkinson G M. 2005. Stochastic finite-fault modeling based on a dynamic corner frequency[J].BullSeismolSocAm, 95（3）: 955--1010.

Ruiz S, Madariaga R. 2013. Kinematic and dynamic inversion of the 2008 northern Iwate earthquake[J].BullSeismolSocAm, 103（2A）: 694--708.

Somerville P, Irikura K, Graves R, Sawada S, Wald D, Abrahamson N, Iwasaki Y, Kagawa T, Smith N, Kowada A. 1999. Characterizing crustal earthquake slip models for the prediction of strong ground motion[J].SeismolResLett, 70（1）: 59--80.

USGS. 2013.M6.6: 56 km WSW of Linqiong, China[EB/OL]. [2015-06-10]. http:∥earthquake.usgs.gov/earthquakes/eventpage/usb000gcdd#general_summary.

Zhang W B, Iwata T, Kojiro K. 2010. Dynamic simulation of the 1999 Chi-Chi, Taiwan, earthquake[J].JGeophysRes, 115（B4）: B04305.

Strong ground motion simulation for the 2013MS7.0 Lushan, China, earthquake

（CollegeofEarthScience,UniversityofChineseAcademyofSciences,Beijing100049,China）

The near source strong ground motions of the 2013MS7.0 Lushan, China, earthquake were simulated using empirical Green’s function （EGF） method. At first, we estimated the amount and location of strong motion gene-ration areas （SMGAs） based on the characteristics of both slip distributions from far-field seismic inversion and the envelopes of recorded acceleration from the main shock, and determined the amount of subfaults on SMGAs referring to the scaling law of asperity areaversusseismic moment introduced by Somervilleetal. Then, we implemented the genetic algorithm searching for the optimized value of above two and other source parameters. Based on the source models, we synthetized the waveforms for the 30 selected stations near the source region. The comparison of the synthetic waveforms with the observed records indicated that they agreed very well with each other, especially for the part of high-frequency larger than 1 Hz. We found that there were two obvious SMGAs on the fault, which take the position that the asperities from far-field seismic inversion take. The combined SMGAs we obtained were smaller than those predicted by extension of the scaling law by Somervilleetal.

strong ground motion; LushanMS7.0 earthquake; empirical Green’s function method; source parameter; waveform comparison

10.11939/jass.2015.04.007.

國家自然科學(xué)基金項(xiàng)目（41274068）和中國科學(xué)院、 國家外國專家局創(chuàng)新團(tuán)隊(duì)國際合作伙伴計劃（KZZD-EW-TZ-19）聯(lián)合資助.

2014-09-13收到初稿， 2015-03-10決定采用修改稿.

e-mail: wenbo@ucas.ac.cn

4.007

P315.01

A

藥曉東， 章文波. 2015. 2013年四川蘆山MS7.0地震強(qiáng)地面運(yùn)動模擬. 地震學(xué)報, 37（4）: 599--616.

Yao X D, Zhang W B. 2015. Strong ground motion simulation for the 2013MS7.0 Lushan, China, earthquake.ActaSeismologicaSinica, 37（4）: 599--616.

doi:10.11939/jass.2015.04.007.