基于廣義反演法的不同位置地震引起盆地放大效應(yīng)的研究

王棪雋 劉啟方

摘要：基于強(qiáng)震觀測記錄，采用廣義反演法和HVSR法分析日本仙臺盆地48個強(qiáng)震臺站以及場地vS30和盆地深度對放大效應(yīng)的影響。研究表明，相比于廣義反演法，HVSR法能夠較好地給出場地主頻，但會明顯低估放大效應(yīng)的幅值；處于盆地外不同位置的淺源地震引起的盆地放大效應(yīng)差異明顯，仙臺盆地南部海域地震引起的放大效應(yīng)最大，盆地北部陸地地震引起的放大效應(yīng)最小;盆地南部和東部地震引起的S波的放大效應(yīng)與vS30的相關(guān)性較強(qiáng)，北部地震的放大效應(yīng)與vS30基本不相關(guān);盆地S波的放大效應(yīng)與盆地深度在0.5～5 Hz頻段內(nèi)的相關(guān)性較強(qiáng)，在0.25～0.5 Hz和5～10 Hz頻段內(nèi)基本不相關(guān)。

關(guān)鍵詞：仙臺盆地; 廣義反演; 放大倍數(shù); 盆地深度; S波

中圖分類號： TU4????? 文獻(xiàn)標(biāo)志碼：A?? 文章編號： 1000-0844（2024）03-0614-11

DOI：10.20000/j.1000-0844.20220820002

Basin amplification effect caused by earthquakes at different locations based on the generalized inversion method：a case study of Sendai Basin

WANG Yanjun， LIU Qifang

（Key Laboratory of Structure Engineering of Jiangsu Province， Suzhou Universityof Science and Technology， Suzhou 215009， Jiangsu， China）

Abstract：?Based on the strong motion records， the amplification effects of 48 stations in Sendai Basin， Japan， were analyzed by using the generalized inversion and HVSR methods. Then， the impacts of site vS30 and basin depth on the amplification effect were analyzed. Results show that compared with the generalized inversion method， the HVSR method can provide the main frequency of a site， but it underestimates the amplitude of the amplification effect. The amplification effects of the basin induced by shallow earthquakes in different locations outside the basin considerably vary. The amplification effect caused by ocean earthquakes in the southern part of Sendai Basin is the largest， whereas that caused by land earthquakes in the northern part of Sendai Basin is the smallest. The amplification effect of S-wave caused by earthquakes in the southern and eastern parts of the basin shows a strong correlation with vS30， whereas that in the northern part shows little correlation with vS30. The amplification effect of the S-wave in the basin is strongly correlated with the depth of the basin in the frequency band of 0.5-5 Hz but is basically irrelevant in the frequency band of 0.25-0.5 and 5-10 Hz.

Keywords：Sendai Basin; generalized inversion method; amplification factor; basin depth; S-wave

0 引言

影響盆地放大效應(yīng)的因素，除了眾所周知的震源（震級、震源深度和斷層類型等）、距離（幾何擴(kuò)散和非彈性衰減）和場地（vS30和覆蓋層厚度）三個因素外，還包括臺站相對于盆地的位置、震源相對于盆地的位置和震源到達(dá)盆地的入射角等因素。

廣義反演法（GIT）由Andrews［1］首先提出，可以分離傅里葉振幅譜中依賴頻率傳播的路徑效應(yīng)、場地效應(yīng)和震源效應(yīng)。一些學(xué)者已通過這種方法來處理全球范圍內(nèi)的地震動數(shù)據(jù)［2-7］。

由于不同位置震源到達(dá)盆地的入射角和路徑不同，所以不同震源引起的盆地地震動放大也不同。已有的研究表明，沉積盆地中存在復(fù)雜的路徑衰減和強(qiáng)烈的場地效應(yīng)［8-13］。強(qiáng)烈地震動很可能是由于盆地產(chǎn)生的面波對體波存在較強(qiáng)的影響，這些影響在盆地邊緣局部存在，已被證實(shí)是地震損害增強(qiáng)的主要原因，如1995年神戶地震發(fā)生的神戶損害帶［11］，以及1994年北嶺地震圣莫尼卡地區(qū)的特別損害帶［8］。盆地的不規(guī)則形狀導(dǎo)致地震傳播的路徑不同，從而產(chǎn)生的聚焦效應(yīng)也不同，因此分析不同位置震源作用下盆地的放大差異有重要價值。

本文將廣義反演法應(yīng)用于仙臺盆地105次地震的3 327組記錄的大量地震數(shù)據(jù)集，以分離震源、場地和路徑效應(yīng);將衰減函數(shù)進(jìn)行參數(shù)化，得到了與頻率相關(guān)的S波品質(zhì)因子;利用ω-2模型對分離的震源譜進(jìn)行解釋［14-15］;比較并分析了利用GIT反演的場地放大倍數(shù)與利用HVSR法計算的場地放大倍數(shù)之間的關(guān)系；重點(diǎn)分析了不同位置震源引起的場地放大倍數(shù)與vS30和盆地深度的關(guān)系。

1 分析方法

1.1 參數(shù)化廣義反演法

參數(shù)化廣義反演法可以分離觀測到的傅里葉振幅譜中依賴頻率的傳播路徑效應(yīng)、場地效應(yīng)和震源效應(yīng)三項因素。觀測記錄的地震動傅里葉振幅譜在頻域中可以表示為以下三項的乘積：

Oij（f，Mi，Rij）=Si（f，Mi）·Pij（f，Rij）·Gj（f）（1）

式中：Oij（f，Mi，Rij）是i震源在j臺站觀測到的震級為Mi，頻率為f的S波傅里葉振幅譜，Rij為i震源到j(luò)臺站的震源距;Si（f，Mi）表示i震源在頻率為f時的震源譜;Pij（f，Rij）為i震源到j(luò)臺站的路徑衰減;Gj（f）為j臺站的場地放大效應(yīng)。對式（1）的兩邊取對數(shù)，可將式（1）線性化為式（2）：

lnOij（f，Mi，Rij）=lnSi（f，Mi）+lnPij（f，Rij）+lnGj（f）（2）

1.2 單步非參數(shù)化廣義反演法

Oth等［5］在研究日本地震動特性時提出了單步非參數(shù)化廣義譜反演方法，以便更穩(wěn)定地分離日本地震動的傳播路徑效應(yīng)、場地效應(yīng)和震源效應(yīng)三項因素。式（2）可以表示為Ax=b的線性方程組，其中b為所有臺站記錄到的S波的傅里葉振幅譜的一維向量;x為震源譜、場地放大倍數(shù)、路徑衰減這些未知項的一維向量;A為表示兩者關(guān)系的稀疏矩陣。矩陣形式如下：

100·…-100·…-100·…

010·…-100·…-100·…

····…-····…-····…

100·…-010·…-010·…

w100·…-000·…-····…

000·…-w300·…-····…

-w22w2-w22·…-000·…-····…

0-w22w2-w22…-····…-····…

····…-····…-····…×

lnP（f，RBin1）

lnP（f，RBinND）…

lnG1（f）

lnGNS（f）…

lnS1（f，M1）

lnSNE（f，MNE）=lnO1（f，M1，R1）lnOij（f，Mi，Rij）00（3）

式（3）中，路徑衰減項（P）、場地放大項（G）、震源項（S）垂直疊加，并以水平虛線分隔;在稀疏矩陣中，由兩條垂直虛線分隔的左、中、右子矩陣，分別與路徑衰減項（P）、場地放大項（G）、震源項（S）相對應(yīng);w1是路徑衰減的權(quán)重系數(shù)，使lnP（f，RBin1）=0，其中RBin1為最小震源距;權(quán)重系數(shù)w2用來約束路徑衰減曲線，使其變得平滑，盡可能排除場地相關(guān)的影響因素;權(quán)重系數(shù)w3是一個限制條件，用來消除場地放大系數(shù)和震源項之間存在的權(quán)衡，對應(yīng)參考場地的場地放大系數(shù)。

2 數(shù)據(jù)篩選及數(shù)據(jù)處理

2.1 地震動記錄篩選

本研究中使用的數(shù)據(jù)均來自日本K-NET和KiK-net（井上），通過下列原則從仙臺盆地周邊篩選數(shù)據(jù)：（1）震級MJMA>4;（2）震源距在5～250 km;（3）震源深度不大于30 km;（4）每個臺站至少有4次地震記錄，每次地震至少被4個臺站記錄。最終數(shù)據(jù)庫包含了105次地震，48個臺站的3 327組記錄。圖1是地震動記錄的震級和震源距的關(guān)系。

為了研究地震震源位置對盆地S波放大效應(yīng)的影響，我們依據(jù)震源相對盆地的空間分布將地震記錄劃分為3個研究區(qū)域，如圖2所示。區(qū)域Ⅰ包含48個臺站在27次M4.0～5.9地震中記錄到的1 101組數(shù)據(jù)，用紅色圓形表示;區(qū)域Ⅱ包含47個臺站在10次M4.0～5.8地震中記錄到的214組數(shù)據(jù)，用黃色圓形表示;區(qū)域Ⅲ包含48個臺站在68次M4.0～6.4地震中記錄到的2 012組數(shù)據(jù)，用藍(lán)色圓形表示。本文以地下剪切波速為vS=1.5 km/s的等值面為盆地底部［16］，定義為Z1.5，深度等值線如圖2（b）所示，其中標(biāo)記為黃色的上三角形為位于盆地外的參考基巖臺，黑色的上三角形為選取的所有臺站。

2.2 地震動記錄的數(shù)據(jù)處理

（1） 對原始記錄進(jìn)行零線校準(zhǔn)，再以0.2～30 Hz為濾波頻帶進(jìn)行濾波，以減少噪聲對記錄的影響。

（2） 采用Husid Plot函數(shù)［17］識別S波的起始時間，計算如下：

Hn（T）=∫T0［a（t）］2dt∫∞0［a（t）］2dt （4）

式中：Hn（T）為Husid Plot函數(shù)，反映了整個加速度時程能量分布;a（t）為加速度時程，以函數(shù)值急劇增加的點(diǎn)為S波的起始時間。

S波的結(jié)束時間，采用累積均方根函數(shù)（CRMS）來識別［18］，計算如下：

CRMS（T）=1T∫T0‖a（t）‖2dt （5）

式中：CRMS（T）為累積均方根函數(shù);a（t）為加速度時程，以函數(shù)值開始下降的點(diǎn)為S波的結(jié)束時間。

為了保證記錄對頻率的分辨率要求，提取到的S波記錄長度至少是4 s，使傅里葉譜的分辨率達(dá)到0.25 Hz，同時為了避免面波對S波的干擾，提取到的S波記錄至多是12 s。在確定S波的起止時間后，在S波前后各多截取10%的S波長度數(shù)據(jù)，并且乘以邊瓣余弦函數(shù)窗，來減少截斷誤差。圖3給出了基于上述方法截取的一條S波地震動記錄。

（3） 直接通過傅里葉變化得到的曲線有很多毛刺，在進(jìn)行反演時這些尖刺峰值可能會影響反演結(jié)果，因此采用0.5 Hz的矩形窗對S波記錄進(jìn)行平滑處理。

2.3 參考臺站的選取

為了消除震源項和場地放大項之間存在的影響，必須將一個合適的臺站作為參考場地，將這個臺站的放大倍數(shù)在所有頻率范圍內(nèi)設(shè)置為2。一般按照下列條件選取作為參考場地的基巖臺站：（1）vS30>600 m/s的場地鉆孔剪切波速;（2）H/V譜比曲線較為平坦且幅值較小。

本研究選取臺站MYGH03為參考臺站，在計算范圍內(nèi)將這個臺站的場地放大系數(shù)假設(shè)為2，即GMYGH03（f）=2。該臺站覆蓋層30 m處的平均剪切波速vS30=934 m/s，可近似認(rèn)為是巖石場地;同時該臺站的H/V譜比曲線較為平坦且幅值較小，滿足參考臺站的選取要求（圖4）。

3 結(jié)果分析和討論

3.1 震源譜及震源參數(shù)的計算

圖5給出了通過廣義反演計算得到的105次地震的位移震源譜。本文計算了這些地震的地震矩密度函數(shù)M0i（f）˙，計算如下［19］：

Si（f）=（2πf）2·RθφVF4πρV3sR0·M0i（f）˙（6）

式中：Rθφ為輻射花樣系數(shù)，通常取平均值0.63［20］;V表示橫波能量分為兩個水平分量，取12;ρ為地殼介質(zhì)的密度，本文取2.7 g/cm3;vS為地殼介質(zhì)平均剪切波速，本文取3.6 km/s;R0為參考距離，通常取1 km;M0i（f）表示第i次地震在頻率f處的地震矩，可按下式計算：

M0i（f）˙=M0i1+ffc2 （7）

式中：M0i為第i次地震的地震矩；fc為拐角頻率。采用Andrews［1］提出的利用震源譜估計拐角頻率fc和震源譜平臺值Ω0的方法，定義了兩個積分常數(shù)SD2和SV2，由以下公式計算得出：

SD2=2∫∞0D2（f）df （8）

SV2=2∫∞0V2（f）df （9）

式中：D（f）為位移譜，可由加速度譜D（f）=A（f）（2πf）2得出;V（f）為速度譜，可由加速度譜V（f）=A（f）2πf得出。

假設(shè)震源譜符合Brune的ω-2震源模型，那么D（f）可表示為：

D（f）=Ω01+ffc2 （10）

聯(lián)立式（8）～（10），可得：

fc=12πSD2SV2 （11）

Ω0=2SD2πfc （12）

同時，可得地震矩M0：

M0=4πρV3sΩ0Rθφ （13）

采用該方法計算了各地震的M0和fc，以最小二乘法擬合數(shù)據(jù)，結(jié)果如圖6所示。

lg（M0）=24.74-2.26lg（fc） （14）

lg（M0）=22.77-2.37lg（fc） （15）

通過與任葉飛［21］得出的結(jié)果［式（15）］進(jìn)行比較，發(fā)現(xiàn)兩者的各項系數(shù)比較接近，說明本文得到的震源譜較為可靠。

3.2 反演的品質(zhì)因子分析

反演的過程中計算了與距離有關(guān)的路徑衰減函數(shù)Pij（f，Rij）。從6～210 km的距離被分成102個2 km寬的分段;在R=6 km處，Pij（f，Rij）在所有頻率范圍設(shè)定為1;衰減函數(shù)Pij（f，Rij）包含著由幾何擴(kuò)散引起的衰減和由Q值引起的衰減，Q值可以通過假設(shè)幾何擴(kuò)散模型通過計算得到。 本文假設(shè)幾何傳播的函數(shù)為：

G（R）=R0Rα （16）

式中：R為震源距;R0為參考距離，由于觀測記錄中的最小震源距為6 km，因此將參考距離R0取為6 km;α為衰減指數(shù)，是S波的幾何擴(kuò)散。利用品質(zhì)因子Q和幾何擴(kuò)展G（R）可以估計路徑衰減函數(shù)。頻率相關(guān)的路徑衰減函數(shù)Pij（f，Rij）可以通過下式計算：

lnPij（f，Rij）-lnG（Rij）=-πfRijQ（f）Vs（17）

式中：Rij為第i次地震到第j個臺站的震源距;Q（f）為品質(zhì)因子;vS為S波平均剪切波速，取3.6 km/s。

Q（f）作為表征地震動衰減特性的函數(shù)，通常表示為指數(shù)形式Q（f）=Q0fα（Q0、α為常數(shù)）。通過對衰減函數(shù)所得結(jié)果采用最小二乘法回歸擬合，得到區(qū)域Ⅰ、Ⅱ、Ⅲ地區(qū)經(jīng)驗(yàn)關(guān)系分別為117f1.0、75.2f0.7、120.2f1.1。如圖7所示，圖中藍(lán)色直線和散點(diǎn)表示區(qū)域Ⅰ的品質(zhì)因子，紅色直線和散點(diǎn)表示區(qū)域Ⅱ的品質(zhì)因子，黑色直線和散點(diǎn)表示區(qū)域Ⅲ的品質(zhì)因子。圖7中還給出了Oth等［5］基于日本東部地區(qū)地震反演得到的Q（f）結(jié)果（黑色虛線）。圖7中可以看出，本文反演得到的區(qū)域Ⅰ、Ⅲ地區(qū)的Q（f）要略大于Oth等［5］計算結(jié)果，區(qū)域Ⅱ地區(qū)的Q（f）要略小于Oth等［5］計算結(jié)果。考慮到Oth等［5］選取的是整個日本火山帶東側(cè)的數(shù)據(jù)，也就是區(qū)域Ⅱ所在位置，這就意味著區(qū)域Ⅰ、Ⅲ地區(qū)S波非彈性衰減要快于東部地區(qū)。另外，本文選取的是震源距范圍6～210 km的強(qiáng)震動記錄，而Oth［5］選取的數(shù)據(jù)震源距均小于150 km。

3.3 反演的場地放大效應(yīng)分析

3.3.1 HVSR法和GIT法得到的場地反應(yīng)比較

在0.1～20 Hz的頻率范圍內(nèi)，以48個臺站中位于盆地中心的15個臺站的場地放大系數(shù)Gj（f）為例與HVSR法比較，如圖8所示。其中實(shí)線為各臺站廣義反演估算的站點(diǎn)放大系數(shù)，藍(lán)色表示區(qū)域Ⅰ，紅色表示區(qū)域Ⅱ，黑色表示區(qū)域Ⅲ;圖中虛線為各區(qū)域通過HVSR法計算的結(jié)果。HVSR法假設(shè)局部場地條件對地面運(yùn)動的垂直分量沒有顯著影響，HVSR法可以用作場地放大的指標(biāo)。以往的研究表明［21］，HVSR法與場地的地質(zhì)條件一致，能夠揭示主要的頻率峰值。然而，通常的結(jié)論是，該方法在振幅水平上有衰減，特別是對高階頻率，且水平和垂直分量都受到散射和吸收造成的能量損失的影響。因此，HVSR除了場地的放大效應(yīng)外，也反映了高頻能量的一些衰減。

從圖8中可以看出，對于部分臺站（如IWTH26、MYG003、MYGH08），GIT法和HVSR法獲得的場地反應(yīng)較一致;但對于大部分臺站，GIT得到的場地反應(yīng)都要大于HVSR得到的結(jié)果，如圖9所示。兩種方法得到場地反應(yīng)差異的主要原因是：Parolai［22］通過數(shù)值模擬及分離不同震源對HVSR的影響，發(fā)現(xiàn)在波阻抗變化大的基巖-土層交界面，S波和P波轉(zhuǎn)換引起垂直方向的場地放大，從而導(dǎo)致HVSR估計的場地反應(yīng)偏小。對于不同地理位置的震源，部分臺站如IWTH26、MYG012、MYG013、MYG014及MYGH08三個區(qū)域所得到的場地放大均較為相似;其余大部分臺站三個區(qū)域通過GIT法所得到的場地放大均有較大差異，對于大部分臺站，區(qū)域Ⅲ產(chǎn)生的場地反應(yīng)最大，區(qū)域Ⅰ其次，區(qū)域Ⅱ產(chǎn)生的場地反應(yīng)最小;而HVSR法得到的結(jié)果并不能看出明顯的差異。

3.3.2 場地放大系數(shù)與場地地質(zhì)條件的關(guān)系

地表土層30 m處的剪切波平均速度（vS30）是用于衡量強(qiáng)運(yùn)動振幅和土壤或巖石分類的重要現(xiàn)場參數(shù)。對于有鉆孔資料的臺站，本文根據(jù)已有的數(shù)據(jù)如圖10 ，通過延拓法（Wang等［23］）計算這些臺站的vS30。在對數(shù)刻度下分別對三個區(qū)域在0.25～0.5 Hz、0.5～1 Hz、1～5 Hz和5～10 Hz的平均場地放大因子通過最小二乘法擬合線性關(guān)系，擬合結(jié)果如圖11所示。擬合公式為：

ln（y）=b1·vS30+b2 （18）

式中：y為放大倍數(shù)；vS30為地表土層30 m處的剪切波平均速度。

結(jié)果顯示，區(qū)域Ⅰ和Ⅲ在4個頻段內(nèi)的平均場地放大與vS30都存在較強(qiáng)的相關(guān)性，其中，0.5～1 Hz和1～5 Hz內(nèi)的相關(guān)性會更強(qiáng)，斜率分別為-3.96×10-3、-4.34×10-3、-4.18×10-3、-4.06×10-3，相關(guān)系數(shù)也都達(dá)到了0.65以上;區(qū)域Ⅱ僅在0.5～1 Hz和1～5 Hz內(nèi)存在較強(qiáng)的相關(guān)性，直線斜率為-3.53×10-3和-2.79×10-3，相關(guān)系數(shù)分別為0.63和0.60，而0.25～0.5 Hz和5～10 Hz內(nèi)基本不相關(guān)，在這兩個頻段內(nèi)的相關(guān)系數(shù)都小于0.5。

3.3.3 場地放大系數(shù)與盆地深度的關(guān)系

盆地的場地放大與震源的位置息息相關(guān)，不同位置震源入射到盆地的角度和盆地形狀的不同，導(dǎo)致盆地的放大效應(yīng)也不同［24-27］。為了分析不同地震位置S波放大與盆地深度的關(guān)系，本文以地下剪切波速vS=1.5 km/s的等值面為盆地底部，定義為Z1.5，在對數(shù)刻度下擬合了Z1.5與不同頻率段內(nèi)三個區(qū)域平均場地放大因子的線性關(guān)系，如圖12所示，擬合公式為：

ln（y）=b1·h+b2 （19）

式中：y為放大倍數(shù);h為盆地深度。

圖12中，每個散點(diǎn)代表一個臺站所在位置的盆地深度和GIT反演的平均場地放大倍數(shù)，直線表示相應(yīng)區(qū)域擬合的線性關(guān)系，藍(lán)色表示區(qū)域Ⅰ，紅色表示區(qū)域Ⅱ，黑色表示區(qū)域Ⅲ。

從總體上看，區(qū)域Ⅲ在各頻段內(nèi)的放大倍數(shù)（斜率均值為1.75×10-3）均大于區(qū)域Ⅰ（1.58×10-3）、Ⅱ（1.262×10-3）的斜率均值;而區(qū)域Ⅰ在0.25～0.5 Hz和1～10 Hz頻段內(nèi)盆地深度小于200 m時放大倍數(shù)略小于區(qū)域Ⅱ，之后隨著盆地深度的增加逐漸大于區(qū)域Ⅱ。

從圖12中可以看出，在0.25～5 Hz頻率段三個區(qū)域的放大倍數(shù)與盆地深度的相關(guān)性均較弱，斜率均值為0.764×10-3;在0.5～1 Hz頻率段，隨著盆地深度的增加，放大倍數(shù)明顯上升，放大倍數(shù)與盆地深度有較強(qiáng)的相關(guān)性，斜率均值為2.92×10-3是0.25～5 Hz時的約4倍。當(dāng)頻率處于1～5 Hz時，放大倍數(shù)與盆地深度依然有較強(qiáng)的相關(guān)性，斜率均值為1.93×10-3是0.25～5 Hz時的2.5倍左右。但是，當(dāng)頻率超過5 Hz時放大倍數(shù)與盆地深度的相關(guān)性開始減弱，在5～10 Hz時，放大倍數(shù)與盆地深度基本不相關(guān)，甚至出現(xiàn)輕微的負(fù)相關(guān)，即放大倍數(shù)隨著盆地深度增加逐漸變小。

4 結(jié)論

本文利用日本K-NET和KiK-net記錄的大量地震動數(shù)據(jù)，選取了仙臺盆地周邊48個臺站在105次地震中觀測的3 327組記錄，通過廣義反演法（GIT）分析了地震波的震源譜、路徑衰減和臺站的場地放大效應(yīng)。通過對震源譜的研究，得出了拐角頻率與地震矩的函數(shù)關(guān)系，并通過與其他研究的結(jié)果比較，證明了震源譜的可靠性。擬合所得的品質(zhì)因子表明仙臺盆地南部和東部震源S波衰減較快且與頻率存在很強(qiáng)的相關(guān)性，北部震源S波衰減較慢。在對得到的場地效應(yīng)進(jìn)行重點(diǎn)分析后，獲得了如下結(jié)論：

（1） 將本文得到的場地放大效應(yīng)數(shù)據(jù)與HVSR法的結(jié)果進(jìn)行了比較，發(fā)現(xiàn)HVSR法與廣義反演技術(shù)得到的函數(shù)形狀吻合良好?？梢钥闯?，兩種方法得到臺站放大倍數(shù)的形狀基本上是相同的，但大部分站點(diǎn)HVSR法的反應(yīng)幅值明顯低于廣義反演的結(jié)果。說明HVSR法能很好地估計場地放大的峰值頻率，但是會低估場地反應(yīng)幅值。

（2） 通過對不同震源位置GIT法獲得場地放大效應(yīng)的分析，發(fā)現(xiàn)盆地對不同震源位置地震波的放大不同。因此，我們在對數(shù)刻度下分別回歸分析三個區(qū)域在0.25～0.5 Hz、0.5～1 Hz、1～5 Hz和5～10 Hz頻段內(nèi)平均場地放大倍數(shù)與盆地深度的線性關(guān)系，發(fā)現(xiàn)盆地南部海域地震的場地放大效應(yīng)明顯大于盆地北部和東部，而盆地東部地震場地反應(yīng)略大于盆地北部。另外，我們發(fā)現(xiàn)在0.5～5 Hz頻段，放大倍數(shù)隨著盆地深度逐漸增大，但是在0.25～0.5 Hz和5～10 Hz頻段內(nèi)兩者的相關(guān)性變得很弱。盆地東部和南部地震所引起的場地反應(yīng)與盆地深度的相關(guān)性較好，覆蓋層越厚，盆地引起的場地反應(yīng)越大。

（3） 本文進(jìn)一步分析了不同位置震源對同一臺站引起的場地放大效應(yīng)與vS30的關(guān)系，發(fā)現(xiàn)三個區(qū)域在0.5～1 Hz和1～5 Hz頻段內(nèi)的平均場地放大倍數(shù)與vS30的相關(guān)性較強(qiáng);在0.25～0.5 Hz和5～10 Hz頻段內(nèi)的相關(guān)性較弱。其盆地北部地震所引起的場地放大效應(yīng)在0.25～0.5 Hz和5～10 Hz頻段內(nèi)無明顯相關(guān)性。

本文僅考慮了不同位置震源對盆地放大效應(yīng)的影響，進(jìn)一步工作可考慮不同深度震源對盆地放大效應(yīng)的影響。

參考文獻(xiàn)（References）

［1］ ANDREWS D J.Objective determination of source parameters and similarity of earthquakes of different size［M］//Earthquake Source Mechanics.Washington D.C.：American Geophysical Union，2013：259-267.

［2］ CASTRO R R，ANDERSON J G，SINGH S K.Site response，attenuation and source spectra of S waves along the Guerrero，Mexico，subduction zone［J］.Bulletin of the Seismological Society of America，1990，80：1481-1503.

［3］ CASTRO R R，PACOR F，SALA A，et al.S wave attenuation and site effects in the region of Friuli，Italy［J］.Journal of Geophysical Research：Solid Earth，1996，101（B10）：22355-22369.

［4］ OTH A，PAROLAI S，BINDI D，et al.Source spectra and site response from S waves of intermediate-depth Vrancea，Romania，earthquakes［J］.Bulletin of the Seismological Society of America，2009，99（1）：235-254.

［5］ OTH A，BINDI D，PAROLAI S，et al.Spectral analysis of K-NET and KiK-net data in Japan，part II：on attenuation characteristics，source spectra，and site response of borehole and surface stations［J］.Bulletin of the Seismological Society of America，2011，101（2）：667-687.

［6］ OTH A.On the characteristics of earthquake stress release variations in Japan［J］.Earth and Planetary Science Letters，2013，377-378：132-141.

［7］ WANG H W，REN Y F，WEN R Z.Source parameters，path attenuation and site effects from strong-motion recordings of the Wenchuan aftershocks （2008—2013） using a non-parametric generalized inversion technique［J］.Geophysical Journal International，2018，212（2）：872-890.

［8］ GRAVES R W.Modeling three-dimensional site response effects in the Marina District Basin，San Francisco，California［J］.Bulletin of the Seismological Society of America，1993，83（4）：1042-1063.

［9］ HANKS T C.Strong ground motion of the San Fernando，California，earthquake：ground displacements［J］.Bulletin of the Seismological Society of America，1975，65（1）：193-225.

［10］ LIU H L，HEATON T.Array analysis of the ground velocities and accelerations from the 1971 San Fernando，California，earthquake［J］.Bulletin of the Seismological Society of America，1984，74（5）：1951-1968.

［11］ KAWASE H.The cause of the damage belt in Kobe：“the basin-edge effect” constructive interference of the direct S-wave with the basin-induced diffracted/rayleigh waves［J］.Seismological Research Letters，1996，67（5）：25-34.

［12］ PITARKA A，IRIKURA K，IWATA T，et al.Three-dimensional simulation of the near-fault ground motion for the 1995 Hyogo-Ken Nanbu （Kobe），Japan，earthquake［J］.Bulletin of the Seismological Society of America，1998，88（2）：428-440.

［13］ SCHNABEL P B.SHAKE-a computer program for earthquake response analysis of horizontally layered sites［J］.Eerc，1972，70（18）：41278.

［14］ BRUNE J N.Tectonic stress and the spectra of seismic shear waves from earthquakes［J］.Journal of Geophysical Research，1970，75（26）：4997-5009.

［15］ BRUNE J N .Journal of geophysical research［J］.American Geophysical Union，1971，36（3）：1766-1181.

［16］ 李春果.基于強(qiáng)震動記錄殘差分析的關(guān)東盆地放大效應(yīng)研究［D］.哈爾濱：中國地震局工程力學(xué)研究所，2019.

LI Chunguo.Research of kanto basin amplification effect based on residual analysis of ground-motion recordings［D］.Harbin：Institute of Engineering Mechanics，China Earthquake Administration，2019.

［17］ RA L.Gravity effects on the earthquake response of yielding? structures［J］.California Institute of Technology，1967，71（16）：31589.

［18］ MCCANN M，SHAH H.Determining strong-motion duration of earthquakes［J］.Bulletin of the Seismological Society of America，1979，69（4）：1816-1825.

［19］ KANAMORI H.Mechanism of tsunami earthquakes［J］.Physics of the Earth and Planetary Interiors，1972，6（5）：346-359.

［20］ 裴順平，許忠淮，汪素云.中國大陸及鄰近地區(qū)上地幔頂部Sn波速度層析成像［J］.地球物理學(xué)報，2004，47（2）：250-256.

PEI Shunping，XU Zhonghuai，WANG Suyun.Sn wave tomography in the uppermost mantle beneath the China continent and adjacent regions［J］.Chinese Journal of Geophysics，2004，47（2）：250-256.

［21］ 任葉飛.基于強(qiáng)震動記錄的汶川地震場地效應(yīng)研究［D］.哈爾濱：中國地震局工程力學(xué)研究所，2014.

REN Yefei.Study on site effect in the Wenchuan earthquake using strong-motion recordings［D］.Harbin：Institute of Engineering Mechanics，China Earthquake Administration，2014.

［22］ PAROLAI S.Comparison of different site response estimation techniques using aftershocks of the 1999 Izmit earthquake［J］.Bulletin of the Seismological Society of America，2004，94（3）：1096-1108.

［23］ WANG H Y，WANG S Y.A new method for estimating vS（30） from a shallow shear-wave velocity profile （depth 30 m）［J］.Bulletin of the Seismological Society of America，2015，105（3）：1359-1370.

［24］ DAVIS P M，RUBINSTEIN J L，LIU K H，et al.Northridge earthquake damage caused by geologic focusing of seismic waves［J］.Science，2000，289（5485）：1746-1750.

［25］ CHOI Y.Empirical model for basin effects accounts for basin depth and source location［J］.Bulletin of the Seismological Society of America，2005，95（4）：1412-1427.

［26］ IWAKI A，IWATA T.Simulation of long-period ground motion in the Osaka sedimentary basin：performance estimation and the basin structure effects［J］.Geophysical Journal International，2010，181（2）：1062-1076.

［27］ 劉啟方，劉曉光.基于強(qiáng)震記錄的大阪盆地S波和次生波放大效應(yīng)研究［J］.自然災(zāi)害學(xué)報，2021，30（6）：32-42.

LIU Qifang，LIU Xiaoguang.Study on the amplification effect of S wave and basin induced wave of Osaka Basin based on strong motion records［J］.Journal of Natural Disasters，2021，30（6）：32-42.

（本文編輯：任 棟）

基金項目：國家自然科學(xué)基金項目（51978434）

第一作者簡介：王棪雋（1997-），男，碩士研究生，主要從事盆地放大研究。 E-mail：542966504@qq.com。

通信作者：劉啟方（1969-），男，研究員，博士，主要從事地震工程研究。E-mail：qifang_liu@126.com。

王棪雋，劉啟方.基于廣義反演法的不同位置地震引起盆地放大效應(yīng)的研究——以仙臺盆地為例［J］.地震工程學(xué)報，2024，46（3）：614-624.DOI：10.20000/j.1000-0844.20220820002

WANG Yanjun，LIU Qifang.Basin amplification effect caused by earthquakes at different locations based on the generalized inversion method： a case study of Sendai Basin［J］.China Earthquake Engineering Journal，2024，46（3）：614-624.DOI：10.20000/j.1000-0844.20220820002