彭宇洋, 王 墩,2

(1. 中國(guó)地質(zhì)大學(xué)(武漢)地球科學(xué)學(xué)院, 湖北 武漢 430074;2. 中國(guó)地質(zhì)大學(xué)湖北巴東地質(zhì)災(zāi)害國(guó)家野外科學(xué)觀測(cè)研究站, 湖北 武漢 430074)

0 引言

一般而言,在大震發(fā)生后,因數(shù)據(jù)及相關(guān)參數(shù)選擇的不確定性,地震震源參數(shù)快速測(cè)定過(guò)程中往往力求簡(jiǎn)單、準(zhǔn)確測(cè)定地震基本參數(shù),來(lái)服務(wù)地震減災(zāi)及快速響應(yīng)。因此,在震后地震參數(shù)反演中往往將地震震源當(dāng)作一個(gè)點(diǎn)源。但是在這種情況下,通常只能得到基本的地震破裂參數(shù)(如震級(jí)、震中/質(zhì)心位置)。如果地震是“多段式”破裂,那么不同破裂階段的地震信息將在單個(gè)點(diǎn)源設(shè)定的反演結(jié)果中不能區(qū)分。這雖然對(duì)震級(jí)大小測(cè)定沒有什么太大影響,但可能使得有些信息被掩蓋。如2008年汶川大地震,該地震面波震級(jí)8.0,破裂尺度超過(guò)300 km,在地震發(fā)生初期,地震應(yīng)急救援力量主要聚焦于震中位置的汶川附近,忽視了其他地區(qū)。震后地震調(diào)查發(fā)現(xiàn),同在一個(gè)斷層面的北川破壞同樣嚴(yán)重。究其原因主要是在地震應(yīng)急中往往把地震設(shè)定為點(diǎn)源,默認(rèn)震中破壞最大。但實(shí)際上地震是沿線或者沿面破裂的,最開始破裂的地方不一定是地震災(zāi)害最嚴(yán)重的地方。

現(xiàn)在發(fā)展起來(lái)的地震波形擬合方法可以獲得精細(xì)的斷層面位錯(cuò)分布,但其往往需要更多的人工設(shè)定參數(shù)。目前較為成熟的獲得斷層面滑移分布的方法為有限斷層反演方法,但其準(zhǔn)確結(jié)果往往在地震震后數(shù)天至一個(gè)月產(chǎn)生,遠(yuǎn)不能滿足大震應(yīng)急響應(yīng)的需要。那么如何在時(shí)效性和結(jié)果精細(xì)準(zhǔn)確之間實(shí)現(xiàn)一個(gè)“兩全其美”的方法就顯得尤為重要。

1 震源破裂反演方法簡(jiǎn)介

實(shí)時(shí)地震研究中使用的方法有傳統(tǒng)地震學(xué)、測(cè)地學(xué)和最新提出的重力波[1-2]方法。這里主要介紹與本研究相關(guān)的地震學(xué)方法和地震波觀測(cè)數(shù)據(jù)。按使用地震觀測(cè)臺(tái)站遠(yuǎn)近不同,實(shí)時(shí)地震研究可大致分為基于近場(chǎng)地震的觀測(cè)數(shù)據(jù)和全球地震波觀測(cè)數(shù)據(jù)兩個(gè)主要研究范疇。基于近場(chǎng)地震觀測(cè)數(shù)據(jù)的實(shí)時(shí)地震學(xué)研究主要包括利用P波和S波的地震預(yù)警技術(shù)及系統(tǒng)架構(gòu)等[3-9]或利用近場(chǎng)或區(qū)域地震波數(shù)據(jù)反演矩張量[10-11]或搜索與之最匹配的理論波形從而求得矩張量[12]?；谌虻卣鸩ㄓ^測(cè)數(shù)據(jù)的研究主要包括W-phase反演、反投影技術(shù)、有限斷層反演等[13-23]。因數(shù)據(jù)易獲取性等原因,這里我們主要聚焦基于全球地震波觀測(cè)數(shù)據(jù)的實(shí)時(shí)地震學(xué)方法。

基于全球地震波觀測(cè)數(shù)據(jù)的實(shí)時(shí)地震學(xué)研究方法主要有基于W-phase的矩張量反演、反投影成像、有限斷層反演技術(shù)和全波形擬合反演技術(shù)。這些方法所需數(shù)據(jù)類型、時(shí)效性、優(yōu)缺點(diǎn)列于表1。因全波形擬合反演費(fèi)時(shí)較多(如GCMT等),常不用于實(shí)時(shí)地震學(xué)研究,故這里略去。

表1 基于全球地震波數(shù)據(jù)的實(shí)時(shí)地震學(xué)常用研究方法簡(jiǎn)介

W-phase因其到時(shí)快(介于P、S震相之間)、不易限幅、結(jié)果相對(duì)穩(wěn)定、反演速度快等特點(diǎn),在大震震后應(yīng)急響應(yīng)及次生災(zāi)害(特別是海嘯風(fēng)險(xiǎn)評(píng)估)中得到了廣泛應(yīng)用[15,24]。該技術(shù)在反演中利用滑動(dòng)時(shí)窗和格點(diǎn)搜索來(lái)獲得矩張量質(zhì)心的最佳位置和時(shí)刻,所以能較為準(zhǔn)確地獲得基于點(diǎn)源模型的矩張量解,為海嘯風(fēng)險(xiǎn)模擬計(jì)算提供科學(xué)依據(jù)[25]。然而,這種基于簡(jiǎn)單點(diǎn)源模型的假定,在提高結(jié)果穩(wěn)定性的同時(shí),也失去了對(duì)復(fù)雜地震事件震源破裂時(shí)空過(guò)程的解析度,使反演結(jié)果產(chǎn)生偏差[26]。

反投影方法自2005年Ishii 等[13]和Kruger等[17]在《Nature》雜志上同期發(fā)表利用該方法成功成像2004年蘇門答臘島MW9.1大地震震源破裂后,在地震學(xué)界掀起了應(yīng)用熱潮,獲得了對(duì)震源破裂參數(shù)及其他特征的進(jìn)一步認(rèn)識(shí)[27-31]。該方法無(wú)需設(shè)定斷層模型,只用通過(guò)簡(jiǎn)單的波形偏移疊加,就可獲得穩(wěn)定可靠的震源空間遷移分布圖。但是該方法不使用格林函數(shù)和反演,得到的結(jié)果只能反映破裂空間遷移信息,其物理意義尚未明晰。

有限斷層反演技術(shù)是三種方法中最成熟的一種。其在20世紀(jì)70—80年代就開始被用來(lái)反演斷層面位錯(cuò)及子事件[16,32-33],并逐漸作為地震震源研究的一種標(biāo)準(zhǔn)手段得到廣泛應(yīng)用[18-22,34-35]。通過(guò)有限斷層反演可以獲得矩震級(jí)、斷層面位錯(cuò)分布、破裂速度等重要震源信息,為海嘯預(yù)警、震后應(yīng)急響應(yīng)和地震物理等研究提供了基礎(chǔ)數(shù)據(jù)。在震后快速(或自動(dòng))有限斷層反演中,由于借助尺度律和初始震級(jí)來(lái)確定模型參數(shù),常常存在模型設(shè)置不合理導(dǎo)致結(jié)果偏差較大或者穩(wěn)定性較差。因此,本工作在現(xiàn)有的W-phase反演方法基礎(chǔ)上[15,22,25],按照有限斷層原理發(fā)展了多點(diǎn)源W-phase方法。通過(guò)余震空間展布獲得粗略的斷層空間特征,在震源區(qū)域設(shè)定多點(diǎn)源,用W-phase反演方法反演各點(diǎn)源的震源機(jī)制及地震矩張量。并將其應(yīng)用到2004年蘇門答臘島MW9.1大地震中。其結(jié)果與該處俯沖帶海溝滑移方向以及歷史地震震源機(jī)制特征都有較好的一致性。

2 基于多點(diǎn)源的W-phase反演方法

在介紹本文所使用的W-phase多點(diǎn)源方法之前,首先介紹一下傳統(tǒng)的多點(diǎn)源方法。Kikuchi等[16,36-37]開發(fā)了一種迭代算法來(lái)反演復(fù)雜破裂過(guò)程的多點(diǎn)源解,其算法的基本操作總結(jié)如下:

(1) 依據(jù)先驗(yàn)信息進(jìn)行準(zhǔn)備操作。對(duì)已知的地震斷層面布設(shè)網(wǎng)格,并計(jì)算每個(gè)網(wǎng)格到臺(tái)站的格林函數(shù),并將格林函數(shù)與假定的震源時(shí)間函數(shù)進(jìn)行卷積。

(2) 通過(guò)已知的格林函數(shù)與臺(tái)站波形,對(duì)每個(gè)網(wǎng)格點(diǎn)進(jìn)行反演求解,得到誤差最小的點(diǎn)源結(jié)果。

(3) 從觀測(cè)值中減去步驟(2)點(diǎn)源結(jié)果與格林函數(shù)的合成值,然后將剩余波形作為觀測(cè)數(shù)據(jù)輸入到下一次迭代。重復(fù)這個(gè)過(guò)程,得到N個(gè)點(diǎn)源,直到達(dá)到可接受的整體波形擬合或沒有進(jìn)一步合理的點(diǎn)源解。

如果將傳統(tǒng)方法運(yùn)用到W-phase反演中,相當(dāng)于對(duì)式(1)的反復(fù)迭代運(yùn)算,每一次得到一組最佳地震矩張量后都從臺(tái)站中將這個(gè)地震矩張量正演形成的波形剔除掉,剩余波形再進(jìn)行一次反演。這樣做有許多優(yōu)點(diǎn),比如能夠按照能量釋放從大到小依次提取相關(guān)震源。但是也有不足之處:先提取出來(lái)的震源可能為了盡可能地減小誤差而包含了后續(xù)震源的地震矩張量、越到后續(xù)提取出的震源結(jié)果可能越來(lái)越不可靠、整體擬合的誤差與點(diǎn)源解的個(gè)數(shù)難以平衡等等。

(1)

(2)

與傳統(tǒng)方法不同的是,本文中離散的震源點(diǎn)是一次性選好,而不是在網(wǎng)格中搜索,反演結(jié)果是一步到位,不是迭代得到。這樣做的好處是不僅僅節(jié)省了運(yùn)算時(shí)間,也避免了傳統(tǒng)迭代法中后續(xù)提取出來(lái)的震源結(jié)果可靠性不佳的影響。

雖然這樣做有許多優(yōu)點(diǎn),但其缺點(diǎn)是需要一次性配置好所有的震源點(diǎn)。這里我們首先依據(jù)余震空間展布信息獲得的粗略斷層空間特征來(lái)等緯度布設(shè)多點(diǎn)源。這樣就可以通過(guò)W-phase反演獲得設(shè)定多點(diǎn)源的地震矩張量,從而得到其震源機(jī)制及相對(duì)釋放能量大小,為定量判定破裂空間展布特征及破裂尺度提供科學(xué)數(shù)據(jù)。

3 2004年MW9.1蘇門答臘大地震結(jié)果及分析

3.1 地震基本信息

2004年12月26日蘇門答臘MW9.1地震是一次毀滅性的地震,引發(fā)了巨大的海嘯,其破裂時(shí)間持續(xù)了8～10 min,破裂長(zhǎng)度延伸了1 300～1 500 km,矩震級(jí) (MW) 為9.1。這使它成為21世紀(jì)目前最大的地震,僅比1960年的智利和1964年的阿拉斯加地震小,它是現(xiàn)代寬帶地震儀和全球定位系統(tǒng) (GPS) 設(shè)備記錄的第一次大型逆沖斷層破裂。在其主震破裂后的三個(gè)小時(shí)內(nèi),發(fā)生了45次大型余震,這些余震分布在整個(gè)大型逆沖斷層面和西安達(dá)曼斷層沿線。其中7次余震在2 h內(nèi)從主震處由南向北連續(xù)單發(fā)[38]。海底的垂向運(yùn)動(dòng)引發(fā)了一場(chǎng)巨大的海嘯,淹沒了蘇門答臘島北部和孟加拉灣周圍的沿海地區(qū),造成了20多萬(wàn)人死亡。

該地震發(fā)生在大洋巖石圈逆沖到大陸邊緣和島弧下方的位置,大型逆沖斷層型地震基本都發(fā)生在這里。根據(jù)古登堡-里氏震級(jí)的頻率關(guān)系,此處小地震比大地震多得多,很少有超過(guò)MW8以上的地震。然而在極少數(shù)情況下,部分地震事件接近或超過(guò)MW9,這樣的大型逆沖斷層地震涉及到斷層寬度(150～300 km)、長(zhǎng)度(500～1 500 km)很大的運(yùn)動(dòng),釋放了俯沖帶中相對(duì)板塊運(yùn)動(dòng)的大部分應(yīng)變能,經(jīng)常產(chǎn)生巨大的海嘯,這大大放大了它們的破壞潛力。在2004年之前,人們認(rèn)為大型逆沖區(qū)地震最可能發(fā)生在相對(duì)年輕的、以高收斂速度俯沖的洋殼區(qū)域[39];1960年智利MW9.5地震和1964年阿拉斯加MW9.2地震都發(fā)生了這種情況。然而,蘇門答臘地震迫使人們重新評(píng)估這些理論[40-42],因?yàn)樗l(fā)生在一個(gè)具有較老洋殼和相對(duì)收斂緩慢的地區(qū),在其破裂帶北部體現(xiàn)的尤為明顯。

本次地震發(fā)生的具體位置在蘇門答臘斷層的北部邊緣,印度洋板塊正在俯沖到蘇門答臘大陸架下面[43]。震中附近的相對(duì)板塊運(yùn)動(dòng)約為50 mm/a,但不垂直于海溝,沿蘇門答臘逆沖斷層的印度—澳大利亞板塊邊界變形帶在10°S～5°N之間廣泛發(fā)育,造成此處板塊構(gòu)造環(huán)境格外復(fù)雜[44-45]。

3.2 地震震源機(jī)制解反演結(jié)果

基于前文中所提到的多點(diǎn)源方法,此次地震多點(diǎn)源W-phase反演依據(jù)余震空間展布特征(圖1)布設(shè)多點(diǎn)源,使用了震中距5°～85°范圍內(nèi)全球地震臺(tái)網(wǎng)(GSN)記錄的地震波形(圖2)來(lái)反演本次地震的震源機(jī)制解,為了測(cè)試和驗(yàn)證方法的可靠性,本研究分別設(shè)定1、2、3、4、5、6個(gè)點(diǎn)源?；?個(gè)點(diǎn)源模型的波形擬合結(jié)果見圖3。

紅色五角星:2004年蘇門答臘MW9.1大地震震中;黑色震源機(jī)制:震后3個(gè)月5級(jí)以上地震GCMT結(jié)果;紅色小圓圈:震后3個(gè)月4級(jí)以上余震分布位置,共計(jì)2 456個(gè)(數(shù)據(jù)源自USGS)。圖1 余震分布圖Fig.1 Distribution of aftershokes

圖2 本次地震所選取的臺(tái)站分布Fig.2 Distribution of stations selected for the study of Sumatra earthquake

多點(diǎn)源的選取按照余震空間展布得到的粗略斷層空間特征進(jìn)行等緯度選擇:在斷層點(diǎn)為6個(gè)時(shí)按照2°間隔進(jìn)行布設(shè);在斷層點(diǎn)為5個(gè)時(shí)按照2.5°間隔進(jìn)行布設(shè);在斷層點(diǎn)為4個(gè)時(shí)按照2.8°間隔進(jìn)行布設(shè);在斷層點(diǎn)為3個(gè)或2個(gè)時(shí)按照4°間隔進(jìn)行布設(shè)。其時(shí)窗的設(shè)置見圖4,從上至下依次為1點(diǎn)源到6點(diǎn)源的時(shí)窗設(shè)置?？v坐標(biāo)為地震矩釋放率(即對(duì)地震矩Mo求取關(guān)于時(shí)間的導(dǎo)數(shù)),橫坐標(biāo)為時(shí)間。

反演結(jié)果如圖5所示。當(dāng)設(shè)定為單點(diǎn)源時(shí)[圖5(a)],震源機(jī)制為低傾角斷層,與哥倫比亞大學(xué)GCMT結(jié)果吻合性較好。 當(dāng)為2個(gè)點(diǎn)源時(shí)[圖5(b)],震源機(jī)制斷層節(jié)面走向在震中附近點(diǎn)源時(shí)為北西方向,遠(yuǎn)離震中點(diǎn)源時(shí)為北北西,與該處中強(qiáng)余震震源機(jī)制走向一致。 3～5點(diǎn)源結(jié)果[圖5(c)～(e)]展現(xiàn)的特征比2個(gè)點(diǎn)源時(shí)更加細(xì)致,總體沒有太大變化,但波形擬合誤差更小。6個(gè)點(diǎn)源時(shí)[圖5(f)],其震源機(jī)制斷層節(jié)面走向的偏轉(zhuǎn)與俯沖帶海溝滑移方向表現(xiàn)出高度一致性,也與周緣的中強(qiáng)震震源機(jī)制一致,證實(shí)了該方法的可靠性及準(zhǔn)確性。

在傳統(tǒng)W-phase反演中,由于使用的是單一點(diǎn)源,反演結(jié)果只能得到一個(gè)震源機(jī)制解,它獲得的是發(fā)震斷層面整體特征。地震震級(jí)越大,破裂尺度往往更長(zhǎng),破裂時(shí)間更久,此時(shí)基于單點(diǎn)源模型反演并不合適,較難獲取復(fù)雜發(fā)震斷層面空間展布形態(tài)的具體特征;在空間上其能量釋放特征也不清晰,震級(jí)測(cè)定相對(duì)多點(diǎn)源往往偏小(表2)。

以左上角臺(tái)站為例:IC:臺(tái)網(wǎng)名,QIZ:臺(tái)站名,LHZ:寬頻帶垂直向,φ:方位角,Δ:震中距,α:臺(tái)站南北分量與正北方向夾角。黑色波形是實(shí)際觀測(cè)波形,紅色波形為理論預(yù)測(cè)波形,兩個(gè)紅點(diǎn)之間是W-phase震相時(shí)窗。通過(guò)計(jì)算實(shí)際觀測(cè)波形與理論預(yù)測(cè)波形的擬合誤差均方根值,當(dāng)設(shè)定為單點(diǎn)源時(shí)值為0.417 94 mm,當(dāng)6點(diǎn)源時(shí)值為0.330 16 mm圖3 臺(tái)站分布和基于6個(gè)點(diǎn)源模型的波形擬合結(jié)果Fig.3 Stations distribution and waveform fitting results based on six point source models

圖4 多點(diǎn)源W-phase時(shí)窗設(shè)置圖Fig.4 Time windows diagram of multiple point W-phase source

從基于6個(gè)點(diǎn)源模型的結(jié)果分析表明,本次蘇門答臘地震的斷層展布從破裂開始(最南邊震源點(diǎn))到破裂結(jié)束(最北邊震源點(diǎn))的斷層面都是低傾角斷層,但是斷層節(jié)面傾向由東北逐漸向東東南方向轉(zhuǎn)變,這與海溝處斷層的滑移方向相吻合(圖5)。本次地震最終得到的震級(jí)為9.0,能量分布在第2個(gè)點(diǎn)達(dá)到最大,后面3個(gè)點(diǎn)能量釋放相對(duì)較少,能量主要集中在9°N以南(表2)。這與現(xiàn)有研究的反演結(jié)果相吻合(圖6、圖7)。

紅色線條:斷層以及俯沖帶運(yùn)動(dòng)方向;紅色震源機(jī)制解:多點(diǎn)源W-phase反演結(jié)果;淺藍(lán)色的震源機(jī)制解:GCMT此次地震結(jié)果;淺綠色的震源機(jī)制解:所有該處大于6級(jí)以上地震并與本次地震成因相同的震源機(jī)制解圖5 多點(diǎn)源W-phase反演結(jié)果圖Fig.5 Results of multiple point W-Phase source inversion

4 討論

對(duì)于W-phase多點(diǎn)源方法,雖然其結(jié)果展現(xiàn)出與其他方法結(jié)果的高度一致性,但是由于其較為依賴人工設(shè)定的點(diǎn)源位置與持續(xù)時(shí)間,尚不能做到完全自動(dòng)、實(shí)時(shí)地運(yùn)行,這是它將來(lái)有待進(jìn)一步發(fā)展的地方。如果我們根據(jù)震后余震序列展示的主震破裂斷層空間展布來(lái)設(shè)置多點(diǎn)源,那么其時(shí)效性將會(huì)從傳統(tǒng)W-phase方法的震后30 min得到結(jié)果延后到震后幾個(gè)小時(shí)甚至幾天才能得到結(jié)果。

表2 多點(diǎn)源W-phase相關(guān)參數(shù)的設(shè)置及其結(jié)果的具體展示Table 2 Setting of relevant parameters of multiple point W-Phase source and specific display of their results

黑色五角星為震源位置;灰色圓圈是GCMT點(diǎn)源結(jié)果,顏色變化表示了滑動(dòng)量的大小[46-51]。結(jié)果源自瑞士地震服務(wù) (SED) (http://www.seismo.ethz.ch/),由Martin Mai匯編得到圖6 有限斷層反演結(jié)果圖Fig.6 Finite-fault inversion results

但是如果基于快速反投影方法所得到的最大破裂點(diǎn)位置來(lái)進(jìn)行點(diǎn)源位置設(shè)置,那么它將和傳統(tǒng)的W-phase方法一樣迅速,因?yàn)榭焖俜赐队胺椒ǖ慕Y(jié)果大概震后13 min加上震源破裂持續(xù)時(shí)間內(nèi)產(chǎn)出[53],而傳統(tǒng)W-phase結(jié)果大概在震后30 min產(chǎn)出。這意味著在等待W-phase波列傳到所需最遠(yuǎn)臺(tái)站之前,我們就能獲得該地震的多點(diǎn)源破裂特征,而且本程序在高性能計(jì)算服務(wù)器環(huán)境下運(yùn)行,運(yùn)算時(shí)間基本在數(shù)十秒內(nèi)完成。多點(diǎn)源W-phase程序?qū)⑼耆梢詫?shí)現(xiàn)實(shí)時(shí)自動(dòng)反演,從而獲得大震破裂空間能量釋放特征。

W-phase多點(diǎn)源反演方法表明:此方法可以快速測(cè)定破裂過(guò)程復(fù)雜的大震,獲得其變化的震源機(jī)制,這是其他基于點(diǎn)源模型的傳統(tǒng)方法所不能達(dá)到的。特別是對(duì)于俯沖帶地震,快速了解地震震源機(jī)制特別重要。因?yàn)樽呋偷卣鹨话悴划a(chǎn)生大海嘯,而傾滑型地震產(chǎn)生強(qiáng)烈海嘯風(fēng)險(xiǎn),如果某俯沖帶附近先發(fā)生了一個(gè)走滑地震,然后隨之觸發(fā)了海溝垂向分量的運(yùn)動(dòng),那么本算法將能快速判定逆沖事件,為海嘯風(fēng)險(xiǎn)評(píng)估及海嘯預(yù)警提供可靠科學(xué)數(shù)據(jù)。

除此以外,多點(diǎn)源W-phase反演的結(jié)果也可以更好協(xié)助快速有限斷層反演工作。在有限斷層反演中需要準(zhǔn)確的震源機(jī)制解(一般由W-phase反演提供)來(lái)布設(shè)斷層面,現(xiàn)在利用多點(diǎn)源W-phase反演能夠給多斷層有限斷層反演提供更加準(zhǔn)確的參數(shù),實(shí)現(xiàn)更精細(xì)的斷層面位錯(cuò)反演。

對(duì)于本算法,除了與快速反投影方法相結(jié)合以外,也可以與空間格點(diǎn)搜索方法相結(jié)合:在獲得快速反投影得到的破裂點(diǎn)后,在每個(gè)破裂點(diǎn)四周布設(shè)加密網(wǎng)格點(diǎn),并對(duì)每個(gè)格點(diǎn)組合都運(yùn)算一遍其不擬合度,找出目標(biāo)范圍內(nèi)的最佳格點(diǎn)組合。這樣將會(huì)對(duì)快速反投影得到的破裂點(diǎn)更加包容,提升本算法的穩(wěn)健性。

5 總結(jié)

基于傳統(tǒng)W-phase方法,本研究開發(fā)了多點(diǎn)源W-phase反演方法。在通過(guò)結(jié)合余震分布結(jié)果來(lái)等緯度設(shè)定多點(diǎn)源后能夠獲得更加精細(xì)的斷層空間展布及能量釋放特征,為深入認(rèn)識(shí)地震物理過(guò)程提供了新手段。

將多點(diǎn)源W-phase反演方法應(yīng)用到2004年蘇門答臘MW9.1大地震后,其結(jié)果較好地反映出了該地震1 200 km長(zhǎng)破裂的分段震源機(jī)制變化及能量釋放特征。從基于6個(gè)點(diǎn)源模型的結(jié)果進(jìn)行分析表明,本次蘇門答臘地震從破裂開始到破裂結(jié)束,斷層面都是低傾角斷層,斷層面傾向由東北逐漸向東東南方向轉(zhuǎn)變,能量釋放在6°N達(dá)到最大,主要集中在9°N以南,這與已知的有限斷層方法、此處歷史地震的結(jié)果以及海溝處斷層的滑移方向都有高度的一致性。

除此以外,本方法還可以基于快速反投影方法所得到的最大破裂點(diǎn)位置來(lái)實(shí)時(shí)配置程序,這將使本文結(jié)果時(shí)效性與傳統(tǒng)W-phase方法時(shí)效性一致,在震后約30 min即可產(chǎn)出可靠的斷層空間展布及能量釋放特征,使多點(diǎn)源W-phase反演方法有應(yīng)用于海嘯預(yù)警及地震應(yīng)急的潛力。