李自芮,馬翀之,崔 瑾,朱鵬濤
(寧夏回族自治區(qū)地震局,寧夏 銀川 750001)
地傾斜觀測是地殼形變測量的手段之一,它能連續(xù)監(jiān)測地殼傾斜速度及方向的變化。在地震孕育的初期,地殼產(chǎn)生彈性形變,隨著應(yīng)力的積累不斷增強,形變加速,介質(zhì)發(fā)生塑性變化,當(dāng)震源體的應(yīng)力達(dá)到其介質(zhì)破裂的臨界值時,震源區(qū)進入非穩(wěn)定狀態(tài),地殼形變也會出現(xiàn)一些轉(zhuǎn)折和加速等特征。總之,在構(gòu)造力的作用下,致使區(qū)域應(yīng)力場增強以及在一些易積累應(yīng)力的構(gòu)造上形成震源體,在地震孕育的整個過程中都伴隨著地殼形變[1]。
銀川臺石英擺傾斜儀于80年代投入觀測,架設(shè)于寧夏賀蘭山東麓斷裂西側(cè)的一山洞里,洞體基巖為花崗片麻巖,覆蓋厚度約50m,進深約100m,無植被。經(jīng)過洞室的不斷改造,觀測資料于1989年之后趨于正常,觀測精度高于10-8量級[2]。
依據(jù)銀川臺臺站位置及其監(jiān)測能力范圍,本文確定研究區(qū)域為銀川臺石英擺傾斜儀觀測資料異常變化與周邊400km范圍內(nèi)發(fā)生的MS5.0以上地震之間的對應(yīng)關(guān)系,并給出地震預(yù)測指標(biāo)。根據(jù)觀測資料精度定義各級地震的監(jiān)測能力范圍為:
本文統(tǒng)計了在銀川臺監(jiān)測能力范圍內(nèi),自1989年以來發(fā)生的MS5.0以上地震震例,并且分析了對應(yīng)時間段內(nèi)石英擺傾斜儀觀測資料的變化情況,結(jié)果表明:在銀川臺周邊400km范圍內(nèi)發(fā)生的MS5.0以上地震前,石英擺傾斜儀觀測資料均有不同程度的異常反應(yīng)。1989年以來銀川臺周邊地區(qū)發(fā)生MS5.0以上地震統(tǒng)計如表1所示。
表1 1989年以來銀川臺周邊地區(qū)發(fā)生MS5.0以上地震統(tǒng)計
值得一提的是:在離銀川臺千里之外的昆侖山口西8.1級強震前,銀川臺石英擺傾斜儀同樣記錄到了信度較高的異常變化[3],2008年四川汶川8.0級地震與銀川臺石英擺傾斜儀觀測資料也存在一定程度的關(guān)聯(lián)。因此,在分析資料異常變化特征與地震對應(yīng)關(guān)系時一并進行了統(tǒng)計分析。
形態(tài)分析是根據(jù)原始資料的變化形態(tài)進行分析的一種方法。通常當(dāng)資料的正常變化規(guī)律被打破而出現(xiàn)異常變化時,若經(jīng)分析斷定非干擾所致,便可認(rèn)為資料出現(xiàn)了異常。
圖1 銀川臺石英擺傾斜儀五日均值曲線Fig.1 The 5-day mean curve of quartz pendulum in Yinchuan Seismic Station
表2 銀川臺石英擺傾斜儀五日均值異常變化特征及對應(yīng)地震情況
銀川臺石英擺傾斜儀觀測曲線存在正弦曲線形態(tài)的年變規(guī)律,兩個方向的變化規(guī)律基本一致:每年NS向先向北傾斜,EW向先向東傾斜,2—4月NS向轉(zhuǎn)為向南傾斜,EW向轉(zhuǎn)為向西傾斜,9—11月NS向又轉(zhuǎn)為向東傾斜,EW向又轉(zhuǎn)為向東傾斜,變化幅度NS向約1.3角秒/年,EW向約1.4角秒/年,總體趨勢上來說NS向北傾,EW向西傾。石英擺傾斜儀形態(tài)曲線在幾次中強地震和昆侖山8.1級地震、汶川8.0級地震前均出現(xiàn)了異常變化,銀川臺石英擺傾斜儀五日均值線圖如圖1所示,異常變化特征以及對應(yīng)地震情況如表2所示。
地球固體潮汐是人們事先可以預(yù)測和精確測量到的一種自然現(xiàn)象。通過理論研究和長期觀測確認(rèn),地殼內(nèi)應(yīng)變積累到一定程度,會使介質(zhì)物性發(fā)生變化,從而導(dǎo)致該地實測的固體潮振幅發(fā)生變化。在正常情況下,假設(shè)地球為剛體時所計算的當(dāng)?shù)氐睦碚摮闭穹鶠锳,與當(dāng)?shù)貙崪y的固體潮振幅A′之比AA′稱為傾斜固體潮振幅因子γ,它是十分穩(wěn)定的[4]。使用M2波的潮汐因子值作為實用指標(biāo),維尼迪可夫調(diào)和分析方法進行計算,潮汐因子值經(jīng)過三月一平均,一月一滑動處理,如果潮汐因子值偏離正常值的1.5倍中誤差,便認(rèn)為自偏離時刻起資料進入異常變化時段。圖2為銀川臺石英擺傾斜儀M2波潮汐因子曲線,從圖中可以看出,M2波潮汐因子值曲線在幾次中強地震前出現(xiàn)異常變化,異常變化特征見表3。
本文使用寧夏形變觀測資料,對形變監(jiān)測能力范圍內(nèi)形變異常的成功率進行分析統(tǒng)計。成功率η定義為[5]:
按觀測資料異常特征與地震對應(yīng)情況,統(tǒng)計了銀川臺石英擺傾斜儀觀測資料異常成功率,結(jié)果如表4所示。
圖2 銀川臺石英擺傾斜儀M2波潮汐因子曲線(三月一平均,一月一滑動)Fig.2 The quartz pendulum tiltmeter M2 wave tidal factor curve in Yinchuan Seismic station(on average once every 3 months, sliding once a month)
表3 銀川臺石英擺傾斜儀M2波潮汐因子異常變化特征及對應(yīng)地震情況
表4 銀川臺石英擺傾斜儀觀測資料異常成功率
綜合以上分析可知,銀川臺石英擺傾斜儀對監(jiān)測能力范圍內(nèi)中強以上地震有一定的監(jiān)測能力,震前異常主要體現(xiàn)為趨勢異常、破年變異常、M2波潮汐因子低于1.5倍中誤差。依據(jù)本文對觀測資料的分析總結(jié),給出了銀川臺石英擺傾斜儀觀測資料預(yù)報地震的預(yù)測指標(biāo)如下:
(1)銀川臺石英擺傾斜儀觀測資料出現(xiàn)長趨勢異常、破年變異常,可判斷未來2年將會發(fā)生MS5.0以上地震;特別注意在趨勢異常背景下出現(xiàn)的破年變異常,可判斷較長時間發(fā)生強震的可能;
(2)用維尼迪可夫調(diào)和分析方法按月計算M2波潮汐因子值,經(jīng)過三月一平均,一月一滑動處理,潮汐因子值若偏離正常值1.5倍中誤差,便認(rèn)為自偏離時刻起進入異常變化時段,異常變化出現(xiàn)明顯轉(zhuǎn)折,可判斷未來1年左右將會發(fā)生MS5.0以上地震。