海底震后形變多基站聯(lián)合估計(jì)模型

薛樹(shù)強(qiáng), 肖圳,*, 董杰, 韓保民, 孫悅,, 楊文龍

1 中國(guó)測(cè)繪科學(xué)研究院地理信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100830

2 山東理工大學(xué)建筑工程與空間信息學(xué)院, 山東淄博 255049

0 引言

隨著以全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System, GNSS)為代表的空間大地測(cè)量觀(guān)測(cè)技術(shù)的發(fā)展和成熟,大地測(cè)量監(jiān)測(cè)網(wǎng)絡(luò)已實(shí)現(xiàn)全球陸域有效覆蓋,且實(shí)現(xiàn)了長(zhǎng)期連續(xù)觀(guān)測(cè)(楊元喜等,2011,2020).然而,占地球表面積71%的海洋還存在大量海底基準(zhǔn)空白和導(dǎo)航定位服務(wù)盲區(qū)(楊元喜等,2017).為此,國(guó)際大地測(cè)量協(xié)會(huì)專(zhuān)門(mén)成立了海洋大地測(cè)量交叉委員會(huì)(Poutanen and Rózsa,2020);同時(shí),我國(guó)在“十三五”期間也開(kāi)展了海底大地測(cè)量基準(zhǔn)試驗(yàn)網(wǎng)建設(shè)(中國(guó)科學(xué)院,2022).新一代大地測(cè)量觀(guān)測(cè)系統(tǒng)(Global Geodetic Observing System, GGOS)致力于地球變化監(jiān)測(cè)并為地震等災(zāi)害過(guò)程研究提供高精度觀(guān)測(cè),但由于絕大部分地震發(fā)生在板塊交接帶,且分布于海底(Bürgmann and Chadwell,2014),導(dǎo)致海底構(gòu)造探測(cè)和地震等災(zāi)害過(guò)程研究在廣闊海域缺少有效約束(Liu等,2023).因此,海底大地測(cè)量對(duì)于全球海底板塊監(jiān)測(cè)、地球動(dòng)力學(xué)和地震等災(zāi)害過(guò)程研究具有重要意義.

最近20多年來(lái),震后形變理論和算法得到了較大的發(fā)展(Sun and Okubo,1993;喬學(xué)軍等,2021;劉泰等,2022),國(guó)內(nèi)外許多學(xué)者基于大地測(cè)量觀(guān)測(cè)對(duì)震后形變開(kāi)展了大量研究工作,例如有學(xué)者利用震后黏彈性模型研究地震后對(duì)周邊區(qū)域地殼形變的影響(余建勝等,2018;劉泰等,2017;陳飛等,2020);還有學(xué)者從Sentinel-1和ALOS-2圖像中獲得同震形變場(chǎng)并進(jìn)行一系列分析(姜衛(wèi)平等,2022);也有學(xué)者利用水準(zhǔn)監(jiān)測(cè)數(shù)據(jù)建模推斷震后形變?cè)春驼鸷笮巫儥C(jī)制(Hao et al.,2012;Iinuma et al.,2012);Tobita結(jié)合日本東北部GNSS時(shí)間序列和指數(shù)、對(duì)數(shù)模型構(gòu)建了震后形變模型(Tobita,2016).在GNSS坐標(biāo)時(shí)間序列震后形變估計(jì)方面,不少學(xué)者研究采用分步估計(jì)法,即先利用震前數(shù)據(jù)估計(jì)站速度,然后通過(guò)去除震后數(shù)據(jù)中站速度引起的位移和同震形變得到震后形變(Gonzalez-Ortega et al.,2014;Tobita,2016).對(duì)GNSS時(shí)間序列進(jìn)行整體建模的模型參數(shù)估計(jì)具有更為嚴(yán)密的理論基礎(chǔ)(Altamimi et al.,2016).此外,GNSS測(cè)站震后形變具有相同的驅(qū)動(dòng)機(jī)制,所以相近測(cè)站具有相似的震后形變,從而導(dǎo)致GNSS時(shí)間序列中存在共模形變,因此已有學(xué)者采用主成分分析法(Principal component analysis, PCA)提取共模形變,且采用兩步法估計(jì)或迭代估計(jì)得到更為合理的震后形變參數(shù)(Savage and Svarc,1997;蘇利娜等,2019).

利用海底大地測(cè)量開(kāi)展地震過(guò)程研究,已成為最近十多年來(lái)的熱點(diǎn)研究課題(Fujiwara et al.,2022;Iinuma et al.,2012;Ozawa et al.,2012;Watanabe et al.,2014;孫悅等,2023).目前,海底大地測(cè)量、海底板塊運(yùn)動(dòng)和擴(kuò)張監(jiān)測(cè)等通?；谌?qū)Ш叫l(wèi)星系統(tǒng)-聲學(xué)(Global Navigation Satellite System-Acoustics, GNSS-A)組合觀(guān)測(cè)方法實(shí)現(xiàn),該方法最早由美國(guó)斯克利普斯海洋學(xué)研究所提出(Spiess,1980).水下定位受復(fù)雜海洋溫度、鹽度等環(huán)境參數(shù)影響,特別是海洋聲速場(chǎng)隨時(shí)間和空間都存在顯著的差異性,給水下精密定位帶來(lái)了極大挑戰(zhàn)(趙建虎和梁文彪,2019;薛樹(shù)強(qiáng)等,2023).為了提高水下定位精度和可靠性,通常需要在等水深的半徑上均勻布設(shè)多個(gè)海底基準(zhǔn)站構(gòu)成區(qū)域海底基準(zhǔn)網(wǎng)(Yokota et al.,2016;Matsumoto and Fujita,2008),例如,美國(guó)學(xué)者提出的在海底布設(shè)3～4個(gè)海底基準(zhǔn)站的方法一直被日本等沿用至今(Spiess,1985a,1985b).該方法經(jīng)過(guò)最近20多年的發(fā)展和完善,其平面定位精度已經(jīng)達(dá)到2～3 cm(Fujita et al.,2006;Watanabe et al.,2020;Yang and Qin,2021),而Sato等利用日本東北大地震前的數(shù)據(jù)對(duì)日本海底地殼運(yùn)動(dòng)觀(guān)測(cè)系統(tǒng)進(jìn)行了評(píng)估,結(jié)果表明,海底基準(zhǔn)站高程定位精度也可達(dá)到10 cm(Sato et al.,2013b).海底基準(zhǔn)站觀(guān)測(cè)設(shè)備由于供電等問(wèn)題需要新舊站更替或替換電池,由于水下作業(yè)的固有成本和技術(shù)難度,通常采用新舊站更替方式保證時(shí)間序列的連續(xù)性,且每年一般只能復(fù)測(cè)3～5次,導(dǎo)致海底基準(zhǔn)站坐標(biāo)時(shí)序樣本相對(duì)較小,無(wú)法與GNSS等實(shí)時(shí)連續(xù)大地測(cè)量觀(guān)測(cè)相比(Yokota et al.,2018).除上述原因之外,水下定位還會(huì)受到測(cè)量船相對(duì)于海底控制點(diǎn)的航跡、聲速剖面誤差的影響(馬越原等,2022;肖圳等,2023;Xue et al., 2023).

日本近20年建立了世界上最為密集的海底地殼運(yùn)動(dòng)觀(guān)測(cè)網(wǎng)絡(luò)(Ozawa et al.,2002;Ozawa et al.,2012;Yokota et al.,2016),該網(wǎng)絡(luò)清晰記錄了一些大型地震引起的形變等,許多學(xué)者在此基礎(chǔ)上對(duì)2011年日本東北部地震(MW9.0)進(jìn)行了研究,例如,有學(xué)者實(shí)現(xiàn)了災(zāi)害過(guò)程模擬(Iinuma et al.,2012; Sun et al.,2014),研究了日本東北部地震導(dǎo)致的同震滑動(dòng)分布和板塊間耦合(Iinuma et al.,2012;Sato et al.,2013a).日本學(xué)者利用海底大地測(cè)量也開(kāi)展了大量地震過(guò)程研究,并在2020年公開(kāi)了7個(gè)海底站的十年尺度監(jiān)測(cè)資料,這為本文研究提供了數(shù)據(jù)條件.本文借鑒國(guó)際地球參考框架2020(ITRF 2020)震后形變研究思想(Altamimi et al.,2016),結(jié)合海底基準(zhǔn)站坐標(biāo)時(shí)序樣本小、存在中斷等實(shí)際問(wèn)題,構(gòu)建了以站速度和震后形變作為公共參數(shù)的多站聯(lián)合時(shí)序分析模型;后又提出了將N-E方向震后形變弛豫因子作為公共參數(shù)的震后形變參數(shù)聯(lián)合估計(jì)模型,并提出了求解站坐標(biāo)、站速率等線(xiàn)性參數(shù)和震后形變非線(xiàn)性參數(shù)混合模型的“一步法”;最后采用日本海底基準(zhǔn)站網(wǎng)觀(guān)測(cè)時(shí)序驗(yàn)證了本文的主要研究結(jié)果.

1 海底震后形變模型

海底震后形變會(huì)影響海底基準(zhǔn)站坐標(biāo)時(shí)序分析,對(duì)海底基準(zhǔn)維持產(chǎn)生不利影響.為此,本文嘗試引入海底基準(zhǔn)參考?xì)v元,并構(gòu)建海底基準(zhǔn)站坐標(biāo)、速度、震后形變等數(shù)聯(lián)合估計(jì)模型.由于海底基準(zhǔn)站局部網(wǎng)通常布設(shè)在地殼構(gòu)造穩(wěn)定區(qū)域,且分布以水深為半徑的圓上,可假設(shè)有K個(gè)海底基站,第i個(gè)海底基站在參考?xì)v元t0時(shí)刻的平面坐標(biāo)為xi0(i=1,2,…,K),所有海底基準(zhǔn)站具有相同的速度v0和相同的震后形變,則t時(shí)刻的海底基準(zhǔn)站網(wǎng)坐標(biāo)可表示為

(1)

其中,v0稱(chēng)為海底基準(zhǔn)站在參考?xì)v元t0的速度,δPSD為震后形變,εi(t)為觀(guān)測(cè)誤差.上述模型將鄰近站點(diǎn)的海底震后形變作為公共參數(shù),因此我們稱(chēng)其為海底震后形變網(wǎng)解模型.考慮海底基準(zhǔn)站定位精度相對(duì)較低,且海底復(fù)測(cè)周期長(zhǎng)導(dǎo)致觀(guān)測(cè)樣本有限,所以本文建議采用上述多站聯(lián)合的網(wǎng)解模型.相對(duì)于GNSS測(cè)站相距幾十千米,海底基準(zhǔn)站間相距僅有幾百米到幾千米,故而將站速度和震后形變參數(shù)作為公共參數(shù)是合理的.相比于基于主成分分析方法的震后共模形變提取方法(蘇利娜等,2019),本文認(rèn)為上述公共參數(shù)法對(duì)于平滑觀(guān)測(cè)噪聲可具有更強(qiáng)的約束性.

海底基準(zhǔn)站隨板塊存在近似線(xiàn)性運(yùn)動(dòng)的同時(shí),還會(huì)受地震等地質(zhì)事件影響,存在非線(xiàn)性運(yùn)動(dòng).震后形變通常采用指數(shù)或?qū)?shù)函數(shù)模擬,如指數(shù)函數(shù)模型可表示為(Altamimi et al.,2016,2018)

(2)

其中,y=[KEKNτ]T,KE,KN分別為E方向和N方向的振幅因子,τ分別為震后形變弛豫因子(我們假設(shè)N方向和E方向具有相同的弛豫因子),T0為地震發(fā)生時(shí)刻.已有學(xué)者利用指數(shù)和對(duì)數(shù)組合模型開(kāi)展震后形變分析,并給出了相應(yīng)參數(shù)的物理意義(Tobita,2016).在實(shí)踐中,分別對(duì)測(cè)站坐標(biāo)時(shí)序分量構(gòu)建震后形變模型,會(huì)得到不同的弛豫時(shí)間因子.然而,對(duì)于同一個(gè)地震事件,不同方向的震后弛豫因子應(yīng)該是相同的.因此,本文采用N-E時(shí)序聯(lián)合估計(jì)地震形變弛豫時(shí)間因子.需要指出,雖然理論上地震會(huì)影響全球測(cè)站坐標(biāo)時(shí)序,但無(wú)論是從計(jì)算成本還是從實(shí)際應(yīng)用方便角度,對(duì)震后形變參數(shù)實(shí)施假設(shè)檢驗(yàn),當(dāng)震后形變參數(shù)顯著時(shí)方可視為模型誤差予以參數(shù)化估計(jì),即上述模型應(yīng)建立在震后形變識(shí)別基礎(chǔ)上.當(dāng)采用假設(shè)檢驗(yàn)方法對(duì)附加參數(shù)進(jìn)行識(shí)別時(shí),需要充分考慮檢驗(yàn)參數(shù)的相關(guān)性問(wèn)題.對(duì)于長(zhǎng)時(shí)序分析問(wèn)題,參數(shù)間的相關(guān)性一般很低,此時(shí)建議采用t檢驗(yàn)對(duì)對(duì)逐個(gè)參數(shù)施加假設(shè)檢驗(yàn)(陶本藻,1986).此外,對(duì)于多維模型參數(shù)選取問(wèn)題,亦可采用更為高效的AIC或BIC信息熵準(zhǔn)則(Akaike,1973;薛樹(shù)強(qiáng)和楊元喜,2013;薛樹(shù)強(qiáng)等,2022).

當(dāng)存在n個(gè)觀(guān)測(cè)歷元時(shí),可構(gòu)建觀(guān)測(cè)方程：

j=1,2,…,n

(3)

其中,xij=xi(tj)為海底基準(zhǔn)站在觀(guān)測(cè)歷元tj的坐標(biāo),可由GNSS-A觀(guān)測(cè)技術(shù)獲取.上述多站聯(lián)合網(wǎng)解模型及N-E時(shí)序聯(lián)合估計(jì)模型,也是解決目前海底基準(zhǔn)站小樣本坐標(biāo)時(shí)序觀(guān)測(cè)和觀(guān)測(cè)精度相對(duì)較低問(wèn)題的有效途徑.

2 海底基準(zhǔn)站坐標(biāo)與速度聯(lián)合估計(jì)模型

Eij=Ei0+vE0(tj-t0)+δEPSD(tj)+εEij,Nij=Ni0+vN0(tj-t0)+δNPSD(tj)+εNij,

(4)

線(xiàn)性化可得

dLi,j=ai,jdxi+bi,jdy,

(5)

其中,dxi=[dEidNidvEdvN]T,dy=[dKEdKNdτ]T,

(6)

(7)

(8)

當(dāng)存在K個(gè)海底基準(zhǔn)站和n個(gè)觀(guān)測(cè)歷元時(shí),即可構(gòu)成如下矩陣形式的觀(guān)測(cè)模型：

dL=CdX+ε,

(9)

(10)

其中,

(11)

為K個(gè)測(cè)站在n個(gè)觀(guān)測(cè)歷元的累計(jì)設(shè)計(jì)矩陣.考慮當(dāng)非線(xiàn)性迭代廣泛存在的不適定問(wèn)題,建議考慮附加一定先驗(yàn)約束.

基于高斯-牛頓法可迭代計(jì)算線(xiàn)性化觀(guān)測(cè)模型(9)的非線(xiàn)性最小二乘解,即

(12)

w(dLi)=

(13)

為等價(jià)權(quán)因子,σ0為觀(guān)測(cè)單位權(quán)方差.IGGIII方案擁有正常權(quán)段、可疑降權(quán)段、以及淘汰權(quán)段,可以充分利用觀(guān)測(cè)數(shù)據(jù),具有較強(qiáng)的抗差性,其中,k1=1.5,k2=2.5為常用推薦值(楊元喜,1996).

采用類(lèi)似的方法,可以分別獲取E方向和N方向的震后形變振幅和弛豫時(shí)間因子,下文不再贅述.此外,擬合震后形變模型也可采用以下“兩步法”,即

(1)首先,忽略模型(5)中的震后形變項(xiàng),獲取觀(guān)測(cè)時(shí)序殘差V,即

V=[I-A(ATPA)ATP]dL;

(14)

(2)將上述殘差向量作為虛擬觀(guān)測(cè)量,采用高斯-牛頓法求解震后形變參數(shù),即

(15)

其中,dSk=V-δPSD(y,t),δPSD(y,t)為全部觀(guān)測(cè)的非線(xiàn)性函數(shù)模型項(xiàng).

需要指出,上述兩步法需要考慮殘差之間的相關(guān)性,并需要迭代計(jì)算以消除分步估計(jì)的影響.對(duì)于小樣本觀(guān)測(cè),忽略殘差的相關(guān)性可能會(huì)對(duì)震后形變參數(shù)估計(jì)產(chǎn)生較大影響.考慮到海底基準(zhǔn)站一般每年復(fù)測(cè)3～5次,觀(guān)測(cè)樣本較小會(huì)導(dǎo)致殘差相關(guān)性較大,本文建議采用理論更為嚴(yán)密的“一步法”,即采用模型(12)一次獲取站坐標(biāo)及震后形變參數(shù).

3 試驗(yàn)結(jié)果分析

2011年3月11日,日本海底發(fā)生了9.0級(jí)大地震,國(guó)內(nèi)外學(xué)者圍繞日本震后地球物理反演開(kāi)展了大量研究工作(Ozawa et al.,2012;Sun et al.,2014;陳飛等,2020).本文利用日本公布的7個(gè)觀(guān)測(cè)站在2011—2020年的GNSS-A實(shí)測(cè)數(shù)據(jù)(空間分布如圖1所示)驗(yàn)證本文方法的有效性.這7個(gè)站均受到2011年日本大地震的影響,以FUKU站為例,E方向的中心點(diǎn)站坐標(biāo)時(shí)序如圖2所示,可以看到在2011年3月日本大地震發(fā)生后坐標(biāo)受到地震影響,隨著時(shí)間的推移逐漸恢復(fù)線(xiàn)性趨勢(shì).

圖1 日本海底基準(zhǔn)站分布

圖2 受震后影響的FUKU站中心點(diǎn)坐標(biāo)時(shí)序

使用本文提出的海底震后形變網(wǎng)解模型以及兩步法分別對(duì)日本公開(kāi)的7個(gè)測(cè)站坐標(biāo)時(shí)序求解各測(cè)站的站速度和震后形變等參數(shù),水平方向的震后形變參數(shù) (振幅、弛豫因子)如表1所示.

表1 兩步法結(jié)果與本文模型解對(duì)比Table 1 Comparison between two-step method and the proposed model solution

由于弛豫因子主要與地殼黏滯系數(shù)和彈性厚度有關(guān),結(jié)合區(qū)域構(gòu)造特征可知(見(jiàn)表1),KAMN、KAMS、MYGI三個(gè)站所在塊體具有相同運(yùn)動(dòng)趨勢(shì),而MYGW、FUKU、CHOS三個(gè)站又與上述三個(gè)站具有反向運(yùn)動(dòng)趨勢(shì),而TOS2距離震中較遠(yuǎn),又表現(xiàn)出于上述三個(gè)站不同的運(yùn)動(dòng)趨勢(shì),然而,兩步法得到的震后弛豫因子多數(shù)均在0.2年左右,他們之間的區(qū)分度很小,這意味著利用小樣本時(shí)序殘差震后形變參數(shù)容易低估弛豫因子的數(shù)值.相比之下,一步法得到的弛豫因子結(jié)果區(qū)分度較好.另外,因?yàn)椴煌５谆鶞?zhǔn)站相對(duì)于地震發(fā)生時(shí)刻的復(fù)測(cè)時(shí)間不同,所以基于不同海底基準(zhǔn)站坐標(biāo)時(shí)序分析得到的震后形變振幅和弛豫因子并不在統(tǒng)一的時(shí)間參考下.為此,我們統(tǒng)計(jì)了日本海底基準(zhǔn)站相對(duì)于地震發(fā)生時(shí)刻的復(fù)測(cè)時(shí)間信息,如圖3所示.

圖3給出的7個(gè)海底基準(zhǔn)站復(fù)測(cè)時(shí)間分布圖表明,TOS2站捕獲的地震形變弛豫因子會(huì)相對(duì)較小,其他站也會(huì)存在微小差異,我們稱(chēng)之為弛豫因子偏差.為此,我們將計(jì)算得到弛豫因子偏差(如表2所示)補(bǔ)償?shù)礁鳒y(cè)站時(shí)序提取的弛豫因子上,得到統(tǒng)一到地震發(fā)生起始時(shí)刻的震后形變弛豫因子.

圖3 7個(gè)海底基準(zhǔn)站復(fù)測(cè)時(shí)間分布

表2 各測(cè)站弛豫因子偏差Table 2 Relaxation factor deviation of each station

將各測(cè)站弛豫因子偏差補(bǔ)償?shù)接?jì)算出來(lái)的弛豫因子參數(shù)后的兩步法與海底震后形變網(wǎng)解模型一步法(分為聯(lián)合估計(jì)N-E方向弛豫因子和分別估計(jì)N、E方向弛豫因子)結(jié)果對(duì)比如表3所示.

表3 補(bǔ)償后兩步法與海底震后形變網(wǎng)解模型參數(shù)結(jié)果對(duì)比Table 3 Comparison of model parameters between compensated two-step method and post earthquake deformation network

由表3可知,兩步法得到的弛豫因子依然沒(méi)有較好的區(qū)分度.需要指出,對(duì)于GNSS基準(zhǔn)站連續(xù)坐標(biāo)時(shí)序的大樣本情形,殘差的相關(guān)性可以忽略,但上述試驗(yàn)結(jié)果表明,對(duì)于海底基準(zhǔn)站復(fù)測(cè)小樣本觀(guān)測(cè)情形,殘差的相關(guān)性可對(duì)震后形變估計(jì)產(chǎn)生很大影響.一步法(7參數(shù)法)即N-E聯(lián)合估計(jì)得到的弛豫因子明顯可以根據(jù)數(shù)據(jù)的大小分為兩組(前三個(gè)站為一組,后四個(gè)站為一組),這主要是因?yàn)?兩組觀(guān)測(cè)站空間分布差異性決定的,即處于兩個(gè)不同的塊體上,這一結(jié)論與實(shí)際站速度觀(guān)測(cè)整體也是基本符合的(如圖4),因此本文認(rèn)為研究海底震后形變參數(shù)時(shí)使用一步法(7參數(shù)法)得到的結(jié)果相比于兩步法要更準(zhǔn)確.

圖4 海底基準(zhǔn)站運(yùn)動(dòng)

4 結(jié)論

對(duì)海底局域網(wǎng)內(nèi)的多個(gè)基準(zhǔn)站坐標(biāo)時(shí)序聯(lián)合處理,不僅可以提高海底板塊運(yùn)動(dòng)的監(jiān)測(cè)精度,也可以提高海底震后形變信息提取的敏感性.受限于觀(guān)測(cè)成本和觀(guān)測(cè)技術(shù)所限,一般每年只能對(duì)海底基準(zhǔn)開(kāi)展3～5次復(fù)測(cè),相對(duì)GNSS而言其復(fù)測(cè)精度也偏低,為此,本文建議將多個(gè)測(cè)站的站速度和震后形變參數(shù)作為公共參數(shù)處理,并采用“一步法”對(duì)線(xiàn)性參數(shù)和非線(xiàn)性參數(shù)進(jìn)行聯(lián)合解算.利用海底基準(zhǔn)站復(fù)測(cè)坐標(biāo)時(shí)序研究震后弛豫因子時(shí),需要充分考慮地震發(fā)生時(shí)刻與海底基準(zhǔn)站復(fù)測(cè)時(shí)刻的偏差影響.

當(dāng)采用“兩步法”提取震后形變信息時(shí),受海底基準(zhǔn)坐標(biāo)觀(guān)測(cè)時(shí)序的樣本大小限制,必然導(dǎo)致其在第一步時(shí)間序列線(xiàn)性擬合階段獲取的殘差時(shí)間序列存在較大的相關(guān)性,若忽略這種相關(guān)性影響,則會(huì)影響其在第二步提取震后形變信息.因此,本文建議采用“一步法”同時(shí)對(duì)測(cè)站坐標(biāo)、速度和震后形變參數(shù)做出聯(lián)合估計(jì).我們認(rèn)為,地震引起地殼形變的弛豫時(shí)間不應(yīng)因坐標(biāo)系的定義不同而不同,因此,為了更進(jìn)一步提高震后形變信息提取的可靠性,我們建議將N、E坐標(biāo)時(shí)序中的震后形變弛豫因子作為公共參數(shù)進(jìn)行估計(jì).

研究發(fā)現(xiàn),日本2011年發(fā)生了MW9.0大地震后,位于不同海底板塊上的海底測(cè)站感應(yīng)到的海底震后形變振幅和弛豫時(shí)間存在較為明顯的差異性,同一板塊上的震后弛豫時(shí)間具有較好的一致性,但當(dāng)采用“兩步法”提取震后弛豫因子時(shí),則很難獲取與之一致的分析結(jié)果.因此本文認(rèn)為,研究海底震后形變參數(shù)時(shí)選擇一步法會(huì)得到比“兩步法”更準(zhǔn)確的結(jié)果.