李 盛 解曉靜 郭昱琴 黃天保
1) 海南省地震局,海南???70203
2) 瓊中基準(zhǔn)地震臺,海南瓊中572900
?
瓊中臺連續(xù)重力觀測記錄地震特征分析*
李盛1)解曉靜1)郭昱琴2)黃天保2)
1) 海南省地震局,海南???70203
2) 瓊中基準(zhǔn)地震臺,海南瓊中572900
摘要利用以瓊中臺為圓心,半徑為3500 km范圍內(nèi)2009—2014年MS≥2.0地震目錄資料,統(tǒng)計(jì)分析瓊中臺連續(xù)重力觀測數(shù)據(jù)曲線與地震的對應(yīng)情況,并對有關(guān)震例在連續(xù)重力原始秒數(shù)據(jù)曲線中產(chǎn)生的地震波進(jìn)行特征分析,總結(jié)規(guī)律,得出瓊中臺連續(xù)重力觀測對震中距3500 km以內(nèi)地震的記錄能力。同時分析了地震事件對重力觀測數(shù)據(jù)、 潮汐和非潮汐參數(shù)的影響,認(rèn)為地震不僅在重力觀測數(shù)據(jù)曲線上產(chǎn)生脈沖或階躍,在非潮汐重力殘差曲線上也會產(chǎn)生脈沖,而對重力M2波潮汐因子產(chǎn)生的影響較小。
關(guān)鍵詞重力觀測; 同震響應(yīng); 震中距; 記錄地震
引言
重力觀測數(shù)據(jù)是地下不同深度、 不同尺度和不同形狀的密度不均勻體所形成的場在地表的綜合反映[1]。臺站重力測量是一種定點(diǎn)觀測,通過對重力加速度的連續(xù)觀測得到重力潮汐資料,通過分析獲得潮汐和非潮汐結(jié)果[2]。經(jīng)過中國地震局“十五”項(xiàng)目和背景場項(xiàng)目的實(shí)施,目前,全國共有75個重力觀測臺站,80套重力儀在運(yùn)行產(chǎn)出數(shù)據(jù),空間分辨率達(dá)400 km,正為我國的地震監(jiān)測預(yù)報(bào)事業(yè)發(fā)揮著重要作用。所謂記錄地震能力,簡單說就是一般情況下,某觀測手段能夠在不同震中距記錄到地震事件的最小震級,或不同震級條件下記錄到地震的最遠(yuǎn)距離[3]。對于連續(xù)重力觀測記錄地震的能力,目前國內(nèi)尚無系統(tǒng)研究和明確結(jié)論。而其他觀測儀器方面,師婭芳等對貴陽臺鉆孔應(yīng)變觀測資料分析認(rèn)為貴陽臺鉆孔應(yīng)變觀測能清晰記錄固體潮變化、 全球中強(qiáng)以上及周邊中小地震同震應(yīng)變變化[4]; 張清秀等對永春井水位觀測資料分析認(rèn)為水位同震效應(yīng)最大震幅總體上與震級成正比,與震中距成反比[5]。
瓊中臺連續(xù)重力觀測項(xiàng)目是中國地震局“十五”重點(diǎn)項(xiàng)目,所用儀器是LCR-PET型重力儀,于2008年3月29日安裝完畢并產(chǎn)出數(shù)據(jù)。
1臺站和儀器概況
瓊中基準(zhǔn)地震臺位于海南省瓊中縣縣城營根鎮(zhèn)東南,處于海南島中部五指山區(qū)木頭山東北側(cè)的小山間盆地內(nèi),屬國家Ⅰ類臺,隸屬于海南省地震局; 受亞熱帶季風(fēng)氣候影響,高溫多雨的氣候特征顯著,年降水量達(dá)2000 mm,年溫度在8~37℃,月溫差為10℃,每年4~10月為雷雨季節(jié),7~10月為臺風(fēng)季節(jié)。臺站周圍10 km內(nèi)無大斷裂通過,地處區(qū)域性基底隆起的花崗巖侵入地區(qū),地殼穩(wěn)定性很好,巖基為海西—印支期花崗巖,巖性致密堅(jiān)硬,該基巖屬混合花崗巖。
瓊中臺連續(xù)重力觀測所用的LCR-PET型重力儀是由美國Micro-g LaCoste,Inc公司生產(chǎn)的一種全自動型重力儀,可用于測量地球各點(diǎn)重力隨時間的變化,也可以檢測周期從秒到若干年的地球物理、 地球動力學(xué)效應(yīng)所導(dǎo)致的微小的重力場變化。儀器主要技術(shù)指標(biāo)如表1所示。
表1 LCR-PET重力儀主要技術(shù)指標(biāo)
2研究資料的選取
2008年3月29日至12月31日,瓊中重力觀測數(shù)據(jù)因儀器故障的原因缺數(shù)較多,在此不作分析。本文研究的瓊中重力觀測數(shù)據(jù)時間區(qū)間是2009年1月至2014年12月。2009年1月以來,瓊中臺重力儀運(yùn)行穩(wěn)定,數(shù)據(jù)可靠,儀器運(yùn)行率達(dá)99.1%,觀測數(shù)據(jù)固體潮清晰,具有明顯的日變形態(tài)(如圖1所示)。
圖1 2009年2月1日瓊中臺重力觀測預(yù)處理分鐘值數(shù)據(jù)曲線
本文所采用的地震目錄主要以中國地震局統(tǒng)一快報(bào)目錄為主要地震目錄來源,對于在研究中涉及到我國境外的地震,則還參考美國地震信息中心地震目錄(從http:∥earthquake.usgs.gov/regional/neic/下載),若二者有不同之處,則以美國地震信息中心地震目錄為準(zhǔn),地震目錄的范圍是2009—2014年以瓊中臺(109.8°E,19.02°N)為圓心,半徑為3500 km范圍內(nèi)所有MS≥2.0地震以及全球MS7.0以上地震,時間為北京時。
3瓊中臺重力觀測記錄地震特征分析
限于篇幅,本特征分析僅對記錄震中距小于等于3500 km的地震進(jìn)行特征分析。下面從地震震中距由近及遠(yuǎn)的角度對記錄地震特征進(jìn)行逐次分析。
2009—2014年該區(qū)域范圍內(nèi)(本區(qū)域范圍都位于我國境內(nèi),完全采用中國地震局統(tǒng)一快報(bào)目錄)共發(fā)生MS2.0~2.9地震5次,MS3.0~3.9地震1次,為2012年11月5日萬寧近海MS3.6地震,未發(fā)生MS4.0級以上地震(震中分布如圖2所示)。瓊中臺重力觀測對上述5次MS2.0~2.9地震未有記錄,但記錄到2012年11月5日萬寧近海MS3.6級地震(該地震信息如表2所示,重力觀測原始秒數(shù)據(jù)曲線如圖3所示),從圖3中可清晰看到地震產(chǎn)生的地震波,且P波和S波清楚,最大振幅時間出現(xiàn)在19:52分,最大振幅為361×10-8m·s-2,因震中距僅129.4 km,地震波到達(dá)時間短,在地震發(fā)生后的30 s內(nèi)地震波即已到達(dá),又因震級較小,同震效應(yīng)持續(xù)時間也比較短,約為3 min。
表2 2012年11月5日萬寧近海MS3.6地震信息
圖2 瓊中臺200 km內(nèi)MS≥2.0地震震中分布圖(2009-01-01~2014-12-31)
圖3 2012年11月5日19:00~20:00時瓊中臺重力觀測原始秒值數(shù)據(jù)曲線
2009—2014年該區(qū)域范圍內(nèi)共發(fā)生MS3.0~3.9地震9次,發(fā)生MS4.0~4.9地震1次,即為2013年2月20日廣西巴馬MS4.5地震。瓊中臺重力觀測未記錄到上述9次MS3.0~3.9地震,但記錄到廣西巴馬MS4.5地震(該地震信息如表3所示,重力觀測原始秒數(shù)據(jù)曲線如圖4所示)。圖4上,因觀測曲線毛刺較多的原因,無法清晰分辨P波到達(dá)時間,地震引起的觀測數(shù)據(jù)最大震幅為67.19×10-8m·s-2,振幅小于萬寧近海MS3.6地震引起的振幅,同震效應(yīng)持續(xù)時間約為5 min。
圖4 2013年2月20日03:00~04:00時瓊中臺重力觀測秒值原始數(shù)據(jù)曲線
表3 2013年2月20日廣西巴馬MS4.5地震信息
2009—2014年該區(qū)域范圍內(nèi)共發(fā)生MS4.0~9.9地震8次,未發(fā)生MS5.0以上地震。上述8次地震中,瓊中臺重力觀測記錄到3次,且該3次地震震級均大于MS4.5(地震信息如表4所示); 其余5次未記錄到,地震震級均介于MS4.0~4.3之間。
表4 瓊中臺重力觀測所記錄到的地震信息(600 km 2009—2014年該區(qū)域范圍內(nèi)共發(fā)生MS4.0~4.9地震155次、MS5.0~5.4地震17次和MS5.5~5.9地震12次。瓊中臺重力觀測對其中的MS4.0~4.9地震均未有記錄; 對MS5.0~5.4地震記錄到12次,5次未記錄到; 對MS5.5~5.9地震均有記錄。由此可以認(rèn)為對于震中距大于1200 km的MS4.9以下地震,瓊中重力觀測無記錄能力。 2009—2014年該區(qū)域范圍內(nèi)共發(fā)生MS5.0~5.4地震67次、MS5.5~5.9地震27次和MS6.0~6.9地震23次。瓊中臺重力觀測對其中的MS5.0~5.4地震均未有記錄;MS5.5~5.9地震中記錄到17次,9次未記錄到,另有1次地震發(fā)生時瓊中臺重力儀無觀測數(shù)據(jù),故無法確定是否可以記錄到該地震;MS6.0~6.9地震均有記錄。由此可以認(rèn)為,對于震中距大于1500 km的MS5.4以下地震,瓊中重力觀測無記錄能力。 2009~2014年該區(qū)域范圍內(nèi)共發(fā)生MS5.5~5.9地震223次,MS6.0~6.4地震110次。瓊中臺重力觀測對于該震中距范圍內(nèi)的MS6.0~6.4地震均有記錄,而對于MS5.5~5.9地震,我們分別統(tǒng)計(jì)了瓊中臺重力觀測對1500 km 以2012年9月7日云南彝良MS5.7和MS5.6雙震(圖5)為例進(jìn)行記錄地震特征分析。該雙震的震中距分別是1115.9 km和1117.3 km,從圖5中可看出,彝良5.7級地震的同震響應(yīng)振幅、 地震引起的脈沖持續(xù)時間明顯大于5.6級地震,說明在震中距相同情況下,地震引起的脈沖振幅、 脈沖持續(xù)時間與震級成正比。從地震發(fā)生到P波初動,耗時約4 min,且地震波結(jié)束后數(shù)據(jù)曲線無階躍,儀器線性漂移速率也未發(fā)生改變。 圖5 2012年9月7日11:00~12:00時瓊中臺重力觀測秒值原始數(shù)據(jù)曲線 表5 瓊中臺重力觀測記錄MS5.5~5.9地震的次數(shù)與未記錄到次數(shù)之比 4地震動對重力儀記錄數(shù)據(jù)的影響 連續(xù)重力觀測的意義在于從重力觀測數(shù)據(jù)中分離出重力潮汐信息和非潮汐信息,對其進(jìn)行研究,發(fā)現(xiàn)可能存在的地震前兆異常, 以應(yīng)用于地震預(yù)測預(yù)報(bào),而不在于記錄地震。地震在重力觀測數(shù)據(jù)曲線上引起的波動嚴(yán)格意義上講是一種干擾,這種干擾是否會影響重力數(shù)據(jù)的應(yīng)用是一個非常關(guān)鍵的問題。重力觀測數(shù)據(jù)曲線記錄地震的表現(xiàn)形式主要有脈沖和同震階躍兩種。對于瓊中臺連續(xù)重力觀測數(shù)據(jù)曲線而言,除2012年4月11日印尼MS8.6地震導(dǎo)致重力儀粘擺(圖6),開擺后觀測數(shù)據(jù)產(chǎn)生階躍外(階躍后的線性漂移速率未發(fā)生改變),其他地震產(chǎn)生的同震響應(yīng)形式均為脈沖。目前為止,重力觀測曲線尚未出現(xiàn)不明脈沖而無地震與之對應(yīng)的情況。 圖6 2012年4月11日至2012年4月13日瓊中臺重力觀測分鐘值預(yù)處理數(shù)據(jù)曲線 地震除了會對重力觀測數(shù)據(jù)產(chǎn)生影響外,對重力潮汐因子和非潮汐重力殘差(所謂非潮汐重力殘差值就是重力觀測數(shù)據(jù)消除理論固體潮和儀器線性漂移后所獲得的數(shù)據(jù))會產(chǎn)生什么樣的影響,是我們在利用連續(xù)重力數(shù)據(jù)進(jìn)行地震分析預(yù)報(bào)中所關(guān)心的。在此,以2014年1月至2014年8月瓊中臺重力預(yù)處理數(shù)據(jù)曲線(圖7上)與經(jīng)計(jì)算后獲得的重力殘差值曲線(圖7中)和重力M2波潮汐因子曲線(圖7下)為例進(jìn)行分析,對于在預(yù)處理曲線上產(chǎn)生明顯同震效應(yīng)的地震在圖上進(jìn)行了標(biāo)注。從圖中可明顯看到,在預(yù)處理曲線上產(chǎn)生明顯同震效應(yīng)的地震在重力殘差值曲線中亦有同震效應(yīng),且在重力殘差值曲線中同震效應(yīng)的表現(xiàn)形式也是脈沖,無階躍情況發(fā)生; 而在M2波曲線上則無明顯的同震效應(yīng),說明地震對重力潮汐參數(shù)的影響較小,對我們利用重力潮汐參數(shù)進(jìn)行地震分析預(yù)報(bào)的影響不大。 圖7 2014年瓊中臺重力預(yù)處理數(shù)據(jù)、 重力殘差值和M2波潮汐因子數(shù)據(jù)曲線 5認(rèn)識與討論 (1) 隨著震中距的增加,瓊中臺重力觀測對地震的記錄能力減弱。其中,對MS<3.0地震無記錄能力; 對震中距超過200 km以上的MS3.0~3.9地震無記錄能力; 對震中距超過800 km以上的MS4.0~4.9地震無記錄能力; 對震中距超過1200 km以上的MS5.0~5.4地震無記錄能力; 對震中距超過3500 km以上MS5.5~5.9地震幾乎無記錄能力。 (2) 瓊中臺連續(xù)重力觀測原始秒數(shù)據(jù)記錄到的地震波具有地震波到達(dá)時間延遲、 對重力觀測曲線產(chǎn)生干擾且干擾的主要表現(xiàn)形式為脈沖、 同震響應(yīng)結(jié)束后觀測曲線的線性漂移不發(fā)生改變等特征。 (3) 瓊中臺重力觀測記錄到的地震事件對重力觀測曲線和非潮汐重力殘差值曲線均會產(chǎn)生脈沖形式的干擾,而對M2波潮汐因子產(chǎn)生的影響較小,無明顯的同震響應(yīng)。 以上認(rèn)識由于地震樣本數(shù)不是很充足、 地震分布不均勻、 時間范圍有限(僅對2009—2014年所發(fā)生的地震進(jìn)行了分析),且未考慮地震震源深度等因素,結(jié)論或有不盡完善之處,今后仍需繼續(xù)研究以對認(rèn)識加以修訂延伸。 (作者電子信箱,李盛: lishyiyuan@163.com) 參 考 文 獻(xiàn) [1] 重力勘探資料解釋手冊編寫組. 重力勘探資料解釋手冊. 北京: 地質(zhì)出版社,1983 [2] 中國地震局監(jiān)測預(yù)報(bào)司. 地形變測量. 北京: 地震出版社,2008: 59-61 [3] 師婭芳,師家升,崔慶谷,等. 貴陽臺鉆孔應(yīng)變觀測記震能力研究. 國際地震動態(tài),2012(9): 20-26 [4] 張清秀,鄭小菁,林木金,等. 永春井水位對遠(yuǎn)大地震的記震能力及同震效應(yīng)特征分析. 華南地震,2011,31(3): 34-40 [5] 宋茉,邱澤華,唐磊. 鉆孔應(yīng)變臺站記震能力研究. 大地測量與地球動力學(xué), 2010,30(增Ⅱ): 32-35 地震科普 The signature of seismic events recorded by the continuous gravimeter at Qiongzhong seismic station Li Sheng1),Xie Xiaojing1),Guo Yuqin2),Huang Tianbao2) 1) Hainan Earthquake Administration,Haikou 570203,China 2) Qiongzhong Seismic Station,Qiongzhong 572900,China AbstractIn this paper,we used the earthquake catalog(MS≥2.0) from 2009 to 2014 within 3500 km away from the Qiongzhong seismic station to analyze statistically the correlation between the continuous gravity observation data recorded by Qiongzhong station and earthquakes. By investigating the characteristic of the seismic waves generated by earthquakes in the original second consecutive gravity data,we get the capabilites of Qiongzhong continuous gravity observation to record the earthquake within 3500 km. At the same time,we analyze the influence of earthquakes to the gravity observation data,tidal and non-tidal data parameters. It is considered the earthquakes not only give rise to a pulse or step on the gravity observation data,but also on a non-tidal gravity residual curve,and cause little influence on the gravityM2 tidal factor. Keywordsgravity observations; the coseismic response; epicentral distance; recorded earthquake 基金項(xiàng)目:中國地震局監(jiān)測、 預(yù)報(bào)、 科研三結(jié)合課題(152102)和中國地震局地震科技星火計(jì)劃(XH15034Y)。 收稿日期:*2015-04-21; 采用日期: 2015-10-10。 中圖分類號:P315.72+6; 文獻(xiàn)標(biāo)識碼:A; doi: 10.3969/j.issn.0235-4975.2015.12.0043.4 對800 km
3.5 對1200 km
3.6 對1500 km