陳學(xué)忠 李艷娥 陳麗娟
1)中國(guó)地震局地球物理研究所,北京 100081
2)重慶市地震局,重慶 401147
(1)預(yù)測(cè)區(qū)域:位于菲律賓板塊南邊界的帛琉群島(Palau)地區(qū)(圖1)。
圖1 預(yù)測(cè)區(qū)域示意圖Fig.1 The area enclosed by the red polygon is the prediction region
(2)預(yù)測(cè)震級(jí):MS≥8.0。
(3)預(yù)測(cè)發(fā)震時(shí)間:2023年12月31日前。
在震源參數(shù)中,視應(yīng)力可作為地殼應(yīng)力水平的間接估計(jì)。所以視應(yīng)力也常被用來分析強(qiáng)震前應(yīng)力變化。很多震例認(rèn)為,強(qiáng)震前震源區(qū)附近視應(yīng)力會(huì)顯著升高[1-4]。
視應(yīng)力的定義為[5]:
式中,μ是剪切模量(對(duì)地殼介質(zhì),μ取3×104MPa),ES和M0分別為地震波輻射能量和地震矩。具體計(jì)算方法詳見文獻(xiàn)[3-4]。
在圖2紅色實(shí)線所圍成的區(qū)域內(nèi),視應(yīng)力大約從2008年初開始從0.4上升,截止到2021年6月,達(dá)到2.2的高值,相對(duì)上升幅度高達(dá)450%。到2021年底持續(xù)時(shí)間達(dá)13年(圖3)。如此長(zhǎng)的持續(xù)時(shí)間,對(duì)應(yīng)的地震震級(jí)可以超過MS8.0。
強(qiáng)震發(fā)生前震源系統(tǒng)會(huì)變得極不穩(wěn)定,平時(shí)對(duì)震源不起作用的外因可能會(huì)起顯著作用,或者說震源系統(tǒng)對(duì)外因敏感了。作用在震源系統(tǒng)上的外力(如,引潮力等)是外因之一。當(dāng)?shù)厍蜃赞D(zhuǎn)速率發(fā)生變化時(shí),其在震源斷層面上引起的應(yīng)力也相應(yīng)發(fā)生變化。由于地球自轉(zhuǎn)速率變化非常小,在震源斷層面上引起的應(yīng)力變化也非常弱[6]。如果震源區(qū)處于極不穩(wěn)定的狀態(tài),這種微弱的應(yīng)力變化或許會(huì)激發(fā)一些地震活動(dòng)。這些被激發(fā)的地震活動(dòng)將會(huì)表現(xiàn)出與地球自轉(zhuǎn)速率變化的顯著相關(guān)性。我們?cè)治隽?976年唐山MS7.8地震、2008年汶川MS8.0地震、2011年日本MW9.1地震及2019年興文MS5.7等強(qiáng)震前破裂區(qū)及附近地區(qū)中小地震與地球自轉(zhuǎn)的相關(guān)性,發(fā)震前小地震均存在低P值異常[7-10]。隨著應(yīng)力的增加,帛琉群島地區(qū)已變得很不穩(wěn)定,其表現(xiàn)是地震活動(dòng)與地球自轉(zhuǎn)速率變化之間具有顯著相關(guān)性。這個(gè)相關(guān)性是否顯著可以用P值來描述,當(dāng)P≤5%時(shí),地震活動(dòng)非隨機(jī)發(fā)生,與地球自轉(zhuǎn)顯著相關(guān)。具體計(jì)算方法詳見文獻(xiàn)[7,11]。圖4a給出了與地球自轉(zhuǎn)存在顯著相關(guān)性的區(qū)域(藍(lán)色實(shí)線圍成的區(qū)域),圖4b給出了P值隨時(shí)間的變化。從2018年初起P值出現(xiàn)快速下降,于年底達(dá)到最大低值~0.05%。
圖2 視應(yīng)力異常區(qū)(圖中紅色實(shí)線所圍成的區(qū)域)Fig.2 Apparent stress anomaly area (solid red rectangle)
圖3 視應(yīng)力時(shí)間變化曲線(圖2中紅色實(shí)線所圍成的區(qū)域)Fig.3 Apparent stress as a function of time (the solid red rectangular area in Figure 2)
圖4 (a)2018年發(fā)生在地球自轉(zhuǎn)加速或減速期間的M≥5.0地震震中分布。加速期間地震集中發(fā)生在藍(lán)色實(shí)線圍成的區(qū)域內(nèi)。(b)P值隨時(shí)間的變化(紅色區(qū)域表示95%置信區(qū)間)Fig.4(a) Spatial distribution of event count for earthquakes (M≥5.0),and earthquakes that occurred during the earth’s rotation are concentrated in the blue rectangle.(b)P-value as a function of time for earthquakes(M≥5.2)that occurred in the blue rectangle(the red area indicates the 95% confidence range)