黃浩，袁伏全，趙玉紅

(青海省地震局，青海 西寧 810001)

0 引 言

地震孕育過程研究的一個基本問題是震源區(qū)周邊的應力調(diào)整過程，即對其空間分布和時間演化特征進行分析，以及對其作為前兆現(xiàn)象的特征進行提取。但由于多方面的原因，目前難以對應力分布情況和應力水平做出一個直接的測量或估計[1]。近年來，地震視應力逐漸被應用于斷裂帶應力狀態(tài)分析與地震趨勢判定[2-9]。

唐古拉地震帶中強地震活躍，近年來發(fā)生了2013年雜多5.1級地震、2015年瑪多5.2級地震、2016年丁青5.5級地震、2016年10月17日雜多6.2級地震、2017年6月唐古拉4級震群、2018年5月6日稱多5.3級地震、2018年8月3日治多5.1級地震、2020年1月25日丁青5.1級地震和2020年4月1日石渠5.6級地震。唐古拉地區(qū)的監(jiān)測能力有限，前兆方面僅玉樹地震臺架設(shè)了觀測手段，在震情跟蹤時常面臨預測指標較少、數(shù)字地震學提取指標研究不深入等問題。

為了提取有效預測指標滿足震情跟蹤需要，根據(jù)監(jiān)測能力篩選唐古拉地震帶東段2012年1月至2020年4月間2.5≤ML≤3.5地震計算地震視應力，研究地震視應力演化特征與5.0級以上地震之間的關(guān)系，提取唐古拉地震帶中強地震視應力預測指標，為震情跟蹤提供科學依據(jù)。

1 資料選取

選取2012年1月至2020年4月青海區(qū)域數(shù)字地震臺網(wǎng)記錄的發(fā)生在唐古拉地震帶東段2.5≤ML≤3.5地震波形(圖1)。在近震源條件下，選用震中距在250 km以內(nèi)的臺站。根據(jù)各臺站記錄的波形情況，選取信噪比較高、記錄較清晰的波形數(shù)據(jù)，將這些臺站記錄到的波形數(shù)據(jù)進行去傾和儀器響應校正。

2 方 法

地震視應力σapp[10]是表征震源區(qū)應力水平的物理量，地震視應力可作為區(qū)域絕對應力水平的下限估計[11]，其定義為：

(1)

其中ES為地震波輻射能量；M0為地震矩；μ為震源區(qū)介質(zhì)剪切模量，通常取3.0×104MPa。

利用波形數(shù)據(jù)計算地震能量和地震矩時，首先將波形數(shù)據(jù)進行去傾處理，然后進行傅氏變換得到觀測譜。在近震源條件下，震源距較小，可以忽略非彈性衰減的影響，通過儀器響應和幾何擴散校正得到震源譜。若忽略非彈性衰減，震源譜[12]可以表示為:

(2)

圖1 唐古拉地震帶中小地震與臺站分布

上式中Ω0為震源譜零頻極限值，fc為拐角頻率。當Ω0和fc給定時可由式(2)確定震源譜[2]。

地震矩M0[13]可以根據(jù)下式求得：

(3)

其中ρ為地殼介質(zhì)密度(取2.71 g/cm3)；ν為波速(S 波取3.5 km/s)；d為震源距；Ω0為震源譜零頻極限值；R為輻射因子，可用其均方根代替(S 波取為0.63)[14]。

地震輻射能量可由對速度譜的平方積分求得[15]：

(4)

根據(jù)上述方法得到每個臺站的地震矩M0、地震波輻射能量ES，即可根據(jù)式(1)對各臺站求平均、并消除個別臺站的異常高值對平均值的影響[16]得到地震視應力。

3 計算結(jié)果及時空特征

選取了至少2個臺站參與計算的事件，最終獲得了278次2.5≤ML≤3.5地震的視應力。由于選取的地震震級相差較小、在一定程度上減小了視應力與震級相關(guān)的影響，因此采用視應力時序變化分析區(qū)域應力波動情況。在分析視應力演化特征時，定義視應力超過1倍均方差為視應力高值異常。

2012年1月至2020年4月，唐古拉地震帶東段共發(fā)生了7次5級以上地震，即2013年雜多5.1級地震、2015年瑪多5.2級地震、2016年丁青5.5級地震、2016年雜多6.2級地震、2018年稱多5.3級地震、2020年丁青5.1級地震和2020年石渠5.6級地震。圖2展示了278次2.5≤ML≤3.5地震的震級和視應力時序變化，結(jié)果顯示視應力數(shù)值在中強地震前均存在明顯異常波動。雖然我們選擇的地震震級相差較小，但震級與視應力仍然存在一定相關(guān)性。因此采用多項式擬合震級與視應力之間的關(guān)系，獲得的擬合關(guān)系為：σapp=0.202 2×ML-0.406 1。為了在一定程度上減小震級與視應力相關(guān)的影響，將震級與視應力擬合關(guān)系作為背景值，再將278次2.5≤ML≤3.5地震的視應力減去背景值，得到規(guī)準化視應力。

從空間分布特征來看，7次5級以上地震均發(fā)生在規(guī)準化視應力異常區(qū)邊緣，地震并沒有發(fā)生在規(guī)準化視應力異常最顯著的地方(由于圖件較多不作逐一展示)。以下分別或成組對7次5級以上地震前后的規(guī)準化視應力時間演化特征進行分析。

圖2 選取資料M-T(a)和視應力時序圖(b)

3.1 2013年雜多5.1級地震

2012年1月7日～2012年2月10日視應力在2個月內(nèi)出現(xiàn)2次高于1倍均方差的高值異常，其后131天規(guī)準化視應力持續(xù)正常，2012年6月22日～2012年7月29日規(guī)準化視應力再次出現(xiàn)高值異常。規(guī)準化視應力演化經(jīng)歷了異常出現(xiàn)—恢復—異常再次出現(xiàn)的過程、持續(xù)202天，異常結(jié)束后6個月發(fā)生2013年1月30日雜多5.1級地震(圖3)。

圖3 選取資料M-T(a)和規(guī)準化視應力時序圖(b)

3.2 2015年瑪多5.2、2016年丁青5.5和雜多6.2級地震

2014年7月20日～2014年8月17日視應力在1個月內(nèi)出現(xiàn)2次高于1倍均方差的高值異常，其后226天規(guī)準化視應力持續(xù)正常，2015年4月2日～2015年5月17日規(guī)準化視應力再次出現(xiàn)高值異常。規(guī)準化視應力演化經(jīng)歷了異常出現(xiàn)—恢復—異常再次出現(xiàn)的過程、持續(xù)301天，異常結(jié)束后5個月、12個月和17個月分別發(fā)生了2015年10月12日瑪多5.2級地震、2016年5月11日丁青5.5級地震和2016年10月17日雜多6.2級地震(圖4)。

圖4 選取資料M-T(a)和規(guī)準化視應力時序圖(b)

3.3 2018年稱多5.3級地震

2018年1月3日～2018年2月8日規(guī)準化視應力在2個月內(nèi)出現(xiàn)2次高于1倍均方差的高值異常，異常結(jié)束后3個月發(fā)生2018年5月6日稱多5.3級地震(圖5)。

圖5 選取資料M-T(a)和規(guī)準化視應力時序圖(b)

3.4 2020年丁青5.1級地震和石渠5.6級地震

2019年2月16日～2019年4月12日規(guī)準化視應力在2個月內(nèi)出現(xiàn)2次高于1倍均方差的高值異常，其后227天規(guī)準化視應力持續(xù)正常，2019年11月27日～2020年1月3日規(guī)準化視應力再次出現(xiàn)高值異常。規(guī)準化視應力演化經(jīng)歷了異常出現(xiàn)—恢復—異常再次出現(xiàn)的過程、持續(xù)321天，異常結(jié)束后1個月和3個月分別發(fā)生了2020年1月25日丁青5.1級地震和2020年4月1日石渠5.6級地震(圖6)。

圖6 選取資料M-T(a)和規(guī)準化視應力時序圖(b)

4 結(jié)論與討論

基于青海區(qū)域數(shù)字地震臺網(wǎng)觀測資料，獲得了2012年以來唐古拉地震帶東部278次2.5≤ML≤3.5地震的視應力。對比分析研究區(qū)域規(guī)準化視應力不同時段的時空特征，獲得以下結(jié)論：

(1)唐古拉地震帶東段2.5≤ML≤3.5地震規(guī)準化視應力，在中強地震前表現(xiàn)出“規(guī)準化視應力異?！謴汀俅纬霈F(xiàn)異?！钡难莼^程，表明中強地震前震源區(qū)域及附近的應力水平相對較高，可以作為唐古拉地震帶中強地震預測指標。規(guī)準化視應力異常演化時間與中強地震對應情況表明，規(guī)準化視應力異常演化過程在200～300天存在發(fā)生5級地震的可能，超過300天形成中強地震成組的可能性較大。

(2)7次中強地震均發(fā)生在規(guī)準化視應力異常區(qū)域邊緣，表明中強地震震中可能不發(fā)生在規(guī)準化視應力異常最高的位置，對中強地震地點預測應關(guān)注整個異常區(qū)及邊緣。

(3)目前，唐古拉地震帶東段視應力演化過程和持續(xù)時間與2014～2015年類似，表明唐古拉地震帶東段近期應力水平相對較高，應密切跟蹤視應力演化過程、注意中強地震的可能。