郝軍麗 朱成林 張正帥 趙 洋 王合乾,2

1 山東省地震局,濟(jì)南市文化東路20號(hào),250014 2 煙臺(tái)地震監(jiān)測(cè)中心站,山東省煙臺(tái)市南山溝17號(hào),264000

北京時(shí)間2023-02-06 09:17土耳其中南部地區(qū)(37.03°E,37.17°N)發(fā)生7.8級(jí)強(qiáng)震,當(dāng)天18:24再次發(fā)生7.8級(jí)強(qiáng)震(37.20°E,38.02°N),2次地震震源深度均較淺,斷層破裂位于地表以下0～20 km。

煙臺(tái)地震監(jiān)測(cè)中心站(簡(jiǎn)稱煙臺(tái)臺(tái))伸縮儀EW分量在土耳其地震前后出現(xiàn)趨勢(shì)轉(zhuǎn)折變化,本文對(duì)該固體潮畸變進(jìn)行分析,初步探討其與土耳其強(qiáng)震的關(guān)系。

1 臺(tái)站及儀器概況

煙臺(tái)市位于山東省東部的膠東半島區(qū)域,地處華北地震區(qū)東南部,地質(zhì)構(gòu)造復(fù)雜。中國(guó)東部規(guī)模最大的郯廬斷裂帶縱貫半島西部,渤海-蓬萊-威海斷裂帶穿越渤海直達(dá)膠東半島北部,并與郯廬斷裂帶在渤海中部交會(huì),膠東半島以南海域存在南黃海斷陷盆地。半島內(nèi)陸還分布著其他許多規(guī)模不等、方向不一的次級(jí)活動(dòng)斷裂。距臺(tái)站10 km遠(yuǎn)處有2條中等規(guī)模的活動(dòng)斷裂:NE向桃村斷裂和NW向蓬萊-威海斷裂。

煙臺(tái)臺(tái)位于煙臺(tái)市南山公園東側(cè)(121°23′E,37°31′N),地處丘陵緩坡園林區(qū),海拔高度45 m,離海約2 200 m,觀測(cè)環(huán)境良好。臺(tái)基以云母角閃片巖為主,山洞覆蓋植被較好,高程約54 m,洞深分別長(zhǎng)約30.60 m和22.60 m,室溫年變幅約為0.2 ℃。常年地下水位埋深為14～17 m。洞內(nèi)布設(shè)水管傾斜儀、水平擺傾斜儀、垂直擺傾斜儀、SS-Y伸縮儀、鉆孔應(yīng)變觀測(cè)、相對(duì)重力觀測(cè)及各種輔助觀測(cè),共計(jì)21個(gè)分量。

煙臺(tái)臺(tái)SS-Y伸縮儀NS分量基線長(zhǎng)21.21 m,EW分量基線長(zhǎng)28.73 m,儀器布設(shè)方位角分別為14°及104°。2017年以來(lái)儀器運(yùn)行基本正常,數(shù)據(jù)質(zhì)量較好,觀測(cè)精度較高。

2 EW分量同震響應(yīng)特征

2023-02-06 09:17土耳其發(fā)生7.8級(jí)地震,煙臺(tái)臺(tái)伸縮儀EW分量同震形變波初動(dòng)時(shí)刻為09:36,與09:28的6.7級(jí)地震重疊,延續(xù)時(shí)間為236 min,變化形態(tài)為振蕩型,同震波形初動(dòng)性質(zhì)為張性,最大變化幅度為3.37×10-8。18:24土耳其再次發(fā)生7.8級(jí)地震,煙臺(tái)臺(tái)伸縮儀EW分量同震形變波動(dòng)初動(dòng)時(shí)刻為18:35,延續(xù)時(shí)間為78 min,變化形態(tài)為振蕩型,同震波形初動(dòng)性質(zhì)為張性,最大變化幅度為2.61×10-8(圖1)。

圖1 煙臺(tái)臺(tái)伸縮儀EW分量記錄的土耳其地震同震波形Fig.1 Co-seismic waveform of tensometer of EW component at Yantai station

3 強(qiáng)震前后EW分量固體潮畸變

3.1 固體潮畸變特征

3.1.1 地震前后固體潮短期畸變特征

煙臺(tái)臺(tái)伸縮儀EW分量2023-02-04開始出現(xiàn)急速上升趨勢(shì)變化,至02-06 01:00達(dá)到峰值,之后開始下降,最大變化幅度達(dá)到9.58×10-8。2023-02-06 09:17、18:24土耳其分別發(fā)生2次7.8級(jí)地震,與煙臺(tái)臺(tái)震中距分別為7 182 km和7 122 km,震源深度20 km。震后固體潮趨勢(shì)短暫回落后仍處于上升趨勢(shì),但速率較震前變緩。圖2為2023-02-01～10煙臺(tái)臺(tái)伸縮儀觀測(cè)值、EW向理論值、水位分鐘值和氣壓分鐘值曲線。

圖2 煙臺(tái)臺(tái)伸縮儀觀測(cè)值、EW向理論值、水位分鐘值和氣壓分鐘值曲線Fig.2 Observation value, theoretical value of EW component of tensometer and minute value curve of water level and air pressure at Yantai station

3.1.2 煙臺(tái)及山東形變觀測(cè)儀器固體潮曲線長(zhǎng)短期趨勢(shì)變化

煙臺(tái)臺(tái)伸縮儀固體潮EW分量出現(xiàn)短期趨勢(shì)畸變的同時(shí),其他觀測(cè)手段未出現(xiàn)同步變化。圖3為煙臺(tái)臺(tái)伸縮儀及水管儀長(zhǎng)期觀測(cè)曲線?？梢钥闯?2022-10開始伸縮儀EW、NS分量及水管儀NS分量出現(xiàn)向上趨勢(shì)變化,水管儀EW分量2022-10出現(xiàn)向下趨勢(shì)變化。表1為趨勢(shì)變化出現(xiàn)前后干擾因素統(tǒng)計(jì),主要有連續(xù)降雨和調(diào)零校準(zhǔn)等因素。其中,短期連續(xù)降雨僅能影響短期趨勢(shì)變化,而調(diào)零校準(zhǔn)有可能改變觀測(cè)曲線趨勢(shì)變化。基于這些干擾因素的存在,需對(duì)2022-10以來(lái)煙臺(tái)臺(tái)形變觀測(cè)的趨勢(shì)轉(zhuǎn)折原因作進(jìn)一步分析。

表1 2022-09～12干擾因素統(tǒng)計(jì)

圖3 煙臺(tái)臺(tái)伸縮儀和水管儀觀測(cè)曲線 Fig.3 Observation curve of tensometer and water-tube tiltmeter at Yantai station

山東地區(qū)有13個(gè)形變臺(tái)站,包含37套形變儀器,共177個(gè)測(cè)項(xiàng),其中擁有洞體應(yīng)變觀測(cè)手段的有4個(gè)臺(tái)站(泰安、郯城、煙臺(tái)、嘉祥)。梳理山東形變觀測(cè)固體潮長(zhǎng)短期趨勢(shì)變化,發(fā)現(xiàn)除煙臺(tái)臺(tái)伸縮儀NS分量、水管儀兩分量(圖3)及榮成鉆孔傾斜儀EW分量(圖4)外,未發(fā)現(xiàn)有同時(shí)段長(zhǎng)短期趨勢(shì)變化。

圖4 榮成臺(tái)鉆孔傾斜儀觀測(cè)值、水位分鐘值和氣壓分鐘值曲線Fig.4 Observation value of drilling inclinometer and minute value curve of water level and air pressure at Rongcheng station

3.1.3 榮成臺(tái)鉆孔傾斜儀畸變

山東榮成地震臺(tái)鉆孔傾斜儀EW分量自2023-01-08～02-06土耳其大震之前固體潮曲線出現(xiàn)上升趨勢(shì)變化,震后02-07趨勢(shì)恢復(fù)正常(圖4)。榮成臺(tái)與煙臺(tái)臺(tái)同屬于膠東半島,且都與土耳其震中處于同一緯度。對(duì)榮成臺(tái)進(jìn)行現(xiàn)場(chǎng)核實(shí)排查,排除了固體潮畸變期間儀器調(diào)試、臺(tái)站周邊建筑施工以及降雨等環(huán)境變化和人為干擾等影響,但由于榮成鉆孔傾斜儀自身穩(wěn)定性較差,難以判定其EW分量畸變?cè)颉?/p>

3.2 儀器穩(wěn)定性

從M2波潮汐因子相對(duì)中誤差、儀器年零漂、數(shù)據(jù)完整率等方面對(duì)煙臺(tái)臺(tái)伸縮儀EW分量進(jìn)行儀器穩(wěn)定性分析(表2、3)。

表2 伸縮儀EW向M2波潮汐因子相對(duì)中誤差

由表2可以看出,伸縮儀EW分量M2波潮汐因子相對(duì)中誤差<0.02,精度較高。表3中2020～2022年伸縮儀EW向年零漂均小于10-6,符合伸縮儀穩(wěn)定性設(shè)計(jì)指標(biāo)。對(duì)伸縮儀運(yùn)行情況進(jìn)行統(tǒng)計(jì)后發(fā)現(xiàn),2020～2022年缺記率分別為0.241%、0.227 5%、0.000 2%,運(yùn)行率、完整率均在99.5%以上。從圖1伸縮儀EW向固體潮分鐘值看出,地震前后固體潮曲線光滑、清晰。

表3 伸縮儀EW向年零漂

綜上,煙臺(tái)臺(tái)SS-Y伸縮儀EW分量運(yùn)行情況良好,觀測(cè)質(zhì)量較高。

3.3 場(chǎng)地環(huán)境及氣象因素

圖5(a)為煙臺(tái)臺(tái)伸縮儀EW分量和洞室溫度日均值長(zhǎng)期曲線。可以看出,伸縮儀所在洞室溫度常年穩(wěn)定在14 ℃左右,室溫年變化約為0.2 ℃,2023-02-01～10趨勢(shì)變化前后洞溫日變化小于0.01 ℃,2項(xiàng)指標(biāo)均符合儀器安裝洞室溫度年變小于0.5 ℃、日變小于0.03 ℃的要求。

圖5 伸縮儀EW分量與洞溫、降雨、氣壓觀測(cè)曲線Fig.5 Observation curve of EW component oftensometer, hole temperature, rainfall and air pressure

圖5(b)為煙臺(tái)臺(tái)伸縮儀EW分量與輔助觀測(cè)降雨及氣壓的整點(diǎn)值曲線。EW分量趨勢(shì)變化會(huì)受到連續(xù)降雨及氣壓變化影響,如2022-09-12臺(tái)風(fēng)“梅花”導(dǎo)致2022-09-14～16持續(xù)降雨,降雨累積量為218 mm,造成形變觀測(cè)出現(xiàn)滯后1～2 d的趨勢(shì)變化;2022-10-02、04降雨后再次出現(xiàn)滯后2 d的趨勢(shì)變化,降雨累積量為111 mm。2023-01以來(lái)無(wú)臺(tái)風(fēng)事件,1月降雨累積量?jī)H17 mm,且2023-02-01～10伸縮儀趨勢(shì)變化前后未有降雨或大風(fēng),在固體潮出現(xiàn)畸變的同一時(shí)段內(nèi)輔助觀測(cè)氣壓曲線未出現(xiàn)同步變化。

經(jīng)過現(xiàn)場(chǎng)核實(shí),畸變出現(xiàn)前后臺(tái)站周圍環(huán)境無(wú)明顯短期變化,觀測(cè)環(huán)境基本一致,且未曾對(duì)儀器進(jìn)行過調(diào)試校準(zhǔn)等人為操作。

綜上,可以排除溫度、短期氣象因素及人為干擾對(duì)此次伸縮儀EW分量的干擾。

3.4 地震活動(dòng)

2023-02全球共發(fā)生12次6級(jí)以上地震,3次7級(jí)以上地震。7級(jí)以上地震分別為02-06土耳其2次7.8級(jí)和02-23塔吉克斯坦7.2級(jí)地震。從時(shí)間上看,02-06土耳其雙震型震群發(fā)生在煙臺(tái)臺(tái)伸縮儀EW分量趨勢(shì)變化之后2 d,土耳其地震屬于雙震型巨震,而塔吉克斯坦7.2級(jí)地震震級(jí)相對(duì)較小,其他地震均小于6.3級(jí),且畸變時(shí)間對(duì)應(yīng)上不如土耳其地震。

2023-02全國(guó)范圍內(nèi)發(fā)生5次4級(jí)以上地震,4次發(fā)生在新疆、青海等西部地區(qū),且震級(jí)較小,與煙臺(tái)臺(tái)趨勢(shì)變化的幅度不對(duì)應(yīng)。山東內(nèi)陸未發(fā)生3級(jí)以上地震,可以排除近震影響。

4 討 論

土耳其7.8級(jí)雙震震中位于地中海-喜馬拉雅地震帶和安納托利亞板塊、阿拉伯板塊及非洲板塊縫合帶,與煙臺(tái)臺(tái)同在37°N附近。煙臺(tái)臺(tái)位于膠東半島、歐亞板塊東緣,受太平洋板塊俯沖作用影響。此次土耳其地震,煙臺(tái)臺(tái)除伸縮儀EW分量出現(xiàn)震前畸變外,其他測(cè)項(xiàng)儀器均未出現(xiàn)短期同步變化。經(jīng)現(xiàn)場(chǎng)排查,煙臺(tái)臺(tái)附近場(chǎng)地環(huán)境無(wú)變化,固體潮變化清晰可見,儀器穩(wěn)定性較好,未見支持此次畸變?cè)虻母蓴_源。

在巖石破壞實(shí)驗(yàn)中發(fā)現(xiàn),臨近主大破裂時(shí)會(huì)出現(xiàn)應(yīng)變擾動(dòng),說明在長(zhǎng)趨勢(shì)形變變化背景上還有突升、突降或類似于波的應(yīng)變擾動(dòng)[1]。地殼是一個(gè)開放的動(dòng)力系統(tǒng),地震是在此系統(tǒng)中孕育發(fā)生的一種局部破裂現(xiàn)象。與有限尺度巖石破裂過程顯著不同,地殼巖石圈近乎為半無(wú)限開放空間,無(wú)明確的應(yīng)力邊界。因此,地震前兆現(xiàn)象未必僅分布在孕震震源區(qū)附近,地震前兆現(xiàn)象,尤其是形變前兆分布的空間范圍往往具有廣泛性、復(fù)雜性。國(guó)內(nèi)外震例表明,地震前在近震源區(qū)往往觀測(cè)不到地形變異?，F(xiàn)象,僅存在極少明確的觀測(cè)異常;而在遠(yuǎn)離震源區(qū)域,往往可觀測(cè)到相對(duì)較多的異常。連續(xù)形變觀測(cè)等的短臨突變性信息可能不是直接來(lái)自于震源,而是源于殼幔邊界或軟流圈內(nèi)長(zhǎng)周期應(yīng)變波的遷移,這種遷移對(duì)地震的發(fā)生具有一定的觸發(fā)或促進(jìn)作用。2023-02-06土耳其7.8級(jí)雙震至今,我國(guó)大陸地區(qū)小震活動(dòng)有一波顯著增強(qiáng)現(xiàn)象,而部分臺(tái)站的某些高精度形變觀測(cè)手段抓取到了相應(yīng)的變形信息[2-5]。

王梅等[6]分析發(fā)現(xiàn),2003年及2004年日本7級(jí)以上地震前,泰安地震臺(tái)伸縮儀EW分量出現(xiàn)固體潮畸變,應(yīng)變量分別達(dá)到3×10-8和6.3×10-8,且其他測(cè)項(xiàng)無(wú)同步變化的情況。該畸變形態(tài)、持續(xù)時(shí)間與前驅(qū)波不同,且變化量級(jí)達(dá)到10-8,比前驅(qū)波應(yīng)變量級(jí)10-9稍大,波形特征不顯著,更類似于靜地震[7-8]的遠(yuǎn)震前兆異常。煙臺(tái)臺(tái)伸縮儀EW分量固體潮變化形態(tài)與典型地震變化及異常形態(tài)相似,且出現(xiàn)在雙強(qiáng)震前,畸變應(yīng)變量為9.58×10-8,出現(xiàn)時(shí)間為震前2 d,并延續(xù)到震后一段時(shí)間,這些特征與泰安臺(tái)伸縮儀畸變類似。

綜上認(rèn)為,此次煙臺(tái)臺(tái)伸縮儀EW分量畸變不排除為土耳其遠(yuǎn)震的前兆異常變化,但在機(jī)理上仍需進(jìn)行深入研究。

5 結(jié) 語(yǔ)

1)經(jīng)核實(shí),本次固體潮畸變期間無(wú)人為干擾、未調(diào)試儀器、場(chǎng)地環(huán)境未發(fā)生變化,儀器數(shù)據(jù)真實(shí)可靠。對(duì)于土耳其7.8級(jí)雙震,煙臺(tái)臺(tái)洞體應(yīng)變資料表現(xiàn)出很好的同震映震性能。

2)煙臺(tái)臺(tái)伸縮儀EW分量穩(wěn)定性指標(biāo)正常,畸變出現(xiàn)前后固體潮曲線清晰、光滑,畸變形態(tài)特征符合地殼應(yīng)力變化特征,認(rèn)為該畸變反映了地殼應(yīng)力的短時(shí)變化。

3)煙臺(tái)臺(tái)伸縮儀NS分量無(wú)短期同步畸變變化,同洞室其他測(cè)項(xiàng)也未見同步變化,可能與儀器自身響應(yīng)特性、地下介質(zhì)特性及儀器系統(tǒng)誤差等因素有關(guān)。

4)此次煙臺(tái)臺(tái)伸縮儀EW分量畸變不排除為土耳其遠(yuǎn)震的前兆異常,但其機(jī)理還需作進(jìn)一步探討。土耳其地震后,煙臺(tái)臺(tái)伸縮儀EW分量固體潮趨勢(shì)短暫回落后繼續(xù)維持上升趨勢(shì),該趨勢(shì)是大震后的能量釋放還是儀器自身的漂移現(xiàn)象,也待進(jìn)一步觀察。

致謝：感謝陳時(shí)軍研究員的幫助!