郭春生， 斯 琴， 關(guān)冬曉， 張嘉敏， 朱治國(guó)

(新疆維吾爾自治區(qū)地震局，新疆 烏魯木齊 830011)

形變是地震孕育的動(dòng)力學(xué)背景、是地震預(yù)報(bào)直接可靠的力學(xué)型前兆，震例總結(jié)表明，中強(qiáng)地震發(fā)生前，地傾斜會(huì)存在不同程度的異?，F(xiàn)象[1-3]。但同一觀測(cè)環(huán)境下，不同儀器，同物理量觀測(cè)結(jié)果不一致[4]、觀測(cè)數(shù)據(jù)間缺少整體性的協(xié)調(diào)對(duì)應(yīng)關(guān)系[5]，仍是困惑地震預(yù)報(bào)的一個(gè)問(wèn)題。另一方面，形變臺(tái)站受到各種因素的影響，會(huì)激勵(lì)、調(diào)制出相應(yīng)的周期性變化，并為觀測(cè)儀器所記錄，因此觀測(cè)值序列必然蘊(yùn)含了多種物理數(shù)學(xué)信息[6]。因此，若觀測(cè)數(shù)據(jù)完全或最大限度地反映的是地球物理場(chǎng)信息或同源干擾信息，那么同臺(tái)站不同的地傾斜觀測(cè)儀器對(duì)同一測(cè)項(xiàng)的觀測(cè)結(jié)果在變化形態(tài)、變幅等方面應(yīng)具有較強(qiáng)的一致性；若觀測(cè)數(shù)據(jù)明顯受到觀測(cè)儀器自身性能的影響，以至于記錄到的地球物理場(chǎng)信息嚴(yán)重失真，則該數(shù)據(jù)應(yīng)用于地震預(yù)報(bào)，是沒(méi)有意義的。因此研究多種觀測(cè)手段下觀測(cè)數(shù)據(jù)的一致性問(wèn)題，以及干擾因素及干擾特征，將有助于認(rèn)識(shí)觀測(cè)數(shù)據(jù)所包含信息物理意義的認(rèn)識(shí)，而這些認(rèn)識(shí)將有助于我們判斷提取的觀測(cè)數(shù)據(jù)的異常信息是來(lái)自于地球內(nèi)部的應(yīng)力積累還是外界的某種因素干擾，這對(duì)震情的跟蹤分析具有積極的意義。

同一觀測(cè)環(huán)境下，不同儀器相同物理量觀測(cè)一致性問(wèn)題及形變觀測(cè)影響因素分析，已有專家進(jìn)行探討，例如：王長(zhǎng)江[8]利用形態(tài)對(duì)比法、潮汐因子幅度、潮汐因子精度等方面分析了觀測(cè)不一致的現(xiàn)象；趙慧琴[4]運(yùn)用形態(tài)對(duì)比法、Venedikov調(diào)和分析法討論了易縣地震臺(tái)形變資料的一致性，結(jié)果表明，同一觀測(cè)環(huán)境下，不同儀器相同物理量觀測(cè)不一致，同物理量觀測(cè)的潮汐因子幅度、精度不一致；王梅[9]對(duì)泰安基準(zhǔn)臺(tái)固體潮形變觀測(cè)各分量的曲線形態(tài)、潮汐參數(shù)等進(jìn)行了對(duì)比分析；孫伶俐[10]運(yùn)用聚類及PCA方法對(duì)湖北省7個(gè)臺(tái)站洞體應(yīng)變觀測(cè)記錄進(jìn)行相關(guān)性分析，歸納總結(jié)了潮汐參數(shù)變化特征差異，并針對(duì)觀測(cè)中自然環(huán)境影響進(jìn)行定量分析，自然環(huán)境因素(氣壓、降水、雷雨、溫度)對(duì)湖北地區(qū)洞體應(yīng)變觀測(cè)影響較為顯著。許璐以溫泉體應(yīng)變?yōu)槔?，利用小波分析和相關(guān)性分析等方法，分析了溫泉體應(yīng)變觀測(cè)各周期變化的影響因素及特征，且不同影響因素對(duì)觀測(cè)儀器的影響的周期特征不同[11]。本文中選取庫(kù)爾勒地震臺(tái)3套地傾斜數(shù)據(jù)為研究對(duì)象，首次嘗試將數(shù)據(jù)分解為日波、半日波、月波、年波等頻段，研究不同頻段信息的一致性及氣溫、氣壓的干擾。

1 臺(tái)站及觀測(cè)系統(tǒng)介紹

庫(kù)爾勒地傾斜觀測(cè)臺(tái)站位于庫(kù)爾勒市北霍拉山南緣的山洞內(nèi)，地理位置為41.82°N,86.19°E,海拔高程1 100 m。附近較大的斷裂主要有北輪臺(tái)—辛格爾斷裂和興地?cái)嗔?。山洞長(zhǎng)約400 m，呈“U”形，觀測(cè)室進(jìn)深200 m，通道呈“L”形，觀測(cè)室上巖石覆蓋厚度近200 m，臺(tái)基巖性為元古界老變質(zhì)巖，觀測(cè)洞室日溫差≤0.1 ℃,年溫差≤0.5 ℃[12]。山洞架有SQ-70D型石英水平擺傾斜儀、MD跨斷層儀、SS-Y型伸縮儀、DSQ型長(zhǎng)水管傾斜儀、GS-15型重力儀、VP寬頻帶垂直擺傾斜儀、斷層氫及輔助觀測(cè)氣象三要素儀(圖1)。

圖1 庫(kù)爾勒地傾斜布設(shè)方位及平面示意圖Fig.1 Orientation and Schematic diagram of Korla crustal inclination

庫(kù)爾勒水平擺2001年起正式觀測(cè)。儀器型號(hào)為SQ-70D，最初為“九五”數(shù)字化觀測(cè)模式，2006年底改造為“十五”模式。2012年6月水平擺更換CCD傳感器，對(duì)資料影響較大。庫(kù)爾勒水平擺固體潮汐清晰、規(guī)則、資料連續(xù)可靠。垂直擺儀器型號(hào)為VP寬頻帶傾斜儀，2014年初開(kāi)始正式觀測(cè)。儀器架設(shè)以來(lái)，工作狀態(tài)穩(wěn)定，資料連續(xù)率、完整率、精度較高。水管儀于2000年11月架設(shè)，2009年2月“十五”數(shù)字化改造，2013年11月比測(cè)(更換不銹鋼缽體、標(biāo)定器、前置放大盒、液體)；2014年12月重新架設(shè)水管儀并換液。

2 數(shù)據(jù)分析方法介紹

2.1 別爾采夫?yàn)V波

別爾采夫?yàn)V波可以把定點(diǎn)形變觀測(cè)數(shù)據(jù)分解成由日月引力作用形成的潮汐部分和儀器的零漂及其他因素引起的低頻趨勢(shì)變化兩部分。其濾波公式為[13]，

(1)

式中，Ym=yn-y-n。

yp為中心時(shí)刻的濾波值，中心時(shí)刻可選在任何一個(gè)時(shí)刻，y0為中心時(shí)刻的觀測(cè)值，yn和y-n分別為以中心時(shí)刻做時(shí)間起算點(diǎn)前n小時(shí)和后n小時(shí)的觀測(cè)值。

2.2 多項(xiàng)式分段曲線擬合

分段曲線擬合方法主要針對(duì)在自然科學(xué)領(lǐng)域中測(cè)量的數(shù)據(jù)而使用的擬合方法，是在工程實(shí)踐與科學(xué)實(shí)驗(yàn)中從一組帶噪聲的實(shí)驗(yàn)觀測(cè)數(shù)據(jù)(xi,yi)i=1,2,3,…n中找出自變量x與因變量y之間隱含的函數(shù)關(guān)系的一種常用方法[14]，但是在數(shù)據(jù)點(diǎn)較多時(shí)，若多項(xiàng)式擬合的階數(shù)低，則擬合精度和效果就不理想，要提高擬合精度和效果就要提高曲線階數(shù)，而階數(shù)太高會(huì)給計(jì)算帶來(lái)麻煩，為了有效的解決該問(wèn)題，一般采用分段曲線擬合,在每段區(qū)間上進(jìn)行局部最小二乘擬合[15]。

3 數(shù)據(jù)分析

本文中選取庫(kù)爾勒地震臺(tái)地傾斜及氣象資料2016～2020年小時(shí)值數(shù)據(jù)為研究對(duì)象。 2016年以來(lái)，水管儀、水平擺NS分量趨勢(shì)N傾，垂直擺S傾， 2018年后水管儀N傾速率增加，水平擺N傾速率減緩；水管儀EW分量E傾，水平擺和垂直擺W傾，庫(kù)爾勒各類地傾斜資料趨勢(shì)性變化存在矛盾。從年變形態(tài)來(lái)說(shuō)(圖2)，水管儀、水平擺和垂直擺均呈現(xiàn)“夏低冬高”的形態(tài)，水管和水平擺的年變形態(tài)較垂直擺年變形態(tài)清晰。氣溫、氣壓年變形態(tài)清晰且趨勢(shì)變化不明顯。

圖2 庫(kù)爾勒地震臺(tái)地傾斜、氣壓、氣溫日均值曲線圖(a) 水管儀NS (b) 水管儀EW (c) 水平擺NS (d) 水平擺EW(e) 垂直擺NS (f) 垂直擺EW (g) 氣溫 (h) 氣壓Fig.2 Curves of ground tilt, air pressure and air temperature at Korla Seismic Station

3.1 庫(kù)爾勒地傾斜日波、半日波一致性及干擾分析

首先利用別爾采夫?yàn)V波將庫(kù)爾勒地傾斜各分量、氣溫、氣壓分解出日波、半日波(圖3)。

圖3 庫(kù)爾勒地傾斜、氣溫、氣壓的日波、半日波曲線圖(a) 水管儀NS (b) 水管儀EW (c) 水平擺NS (d) 水平擺EW(e) 垂直擺NS (f) 垂直擺EW (g) 氣溫 (h) 氣壓Fig.3 Curves of daily waves, half-day waves of ground tilt, temperature and pressure of korla

計(jì)算庫(kù)爾勒地傾斜各分量日波、半日波與氣溫、氣壓日波、半日波(潮汐頻段)及理論固體潮的相關(guān)系數(shù)，結(jié)果見(jiàn)表1。不失一般性，繪制水平擺與垂直擺NS分量，水管儀與垂直擺EW向散點(diǎn)圖(圖4)。

表1 庫(kù)爾勒地傾斜與氣溫、氣壓、固體潮相關(guān)系數(shù)(日波、半日波頻段)

由表1及圖4可知，庫(kù)爾勒地傾斜日波、半日波頻段觀測(cè)結(jié)果一致性較好，且均與理論固體潮具有較高的線性相關(guān)性，而與氣溫、氣壓的相關(guān)性較差。

圖4 庫(kù)爾勒地傾斜不同分量間的散點(diǎn)圖(a) 水平擺與垂直擺NS分量 (b) 水管儀與垂直擺EW向Fig.4 Scatter plot between different components of Korla ground tilt

3.2 庫(kù)爾勒地傾斜月波、年波頻段的一致性及干擾分析

利用一般多項(xiàng)式分段曲線擬合法[16]將剩余的低頻成分濾出周期在1～2個(gè)月內(nèi)的月波(圖5)、年波(圖6)。計(jì)算2016年1月至2020年4月庫(kù)爾勒地傾斜各分量、氣溫、氣壓月波之間的相關(guān)系數(shù)，結(jié)果見(jiàn)表2。

圖5 庫(kù)爾勒地傾斜、氣溫、氣壓月波(2018-01-01～2018-02-28)(a) 水管儀NS (b) 水管儀EW (c) 水平擺NS (d) 水平擺EW(e) 垂直擺NS (f) 垂直擺EW (g) 氣溫 (h)氣壓Fig.5 Moon wave of crustal inclination, temperature and pressure of korla

圖6 庫(kù)爾勒地傾斜、氣溫、氣壓年波(a) 水管儀NS (b) 水管儀EW (c) 水平擺NS (d) 水平擺EW(e) 垂直擺NS (f) 垂直擺EW (g)氣溫 (h) 氣壓Fig.6 Year wave of ground tilt, temperature and pressure of korla

由圖5、表2可知，垂直擺EW分量與氣壓月波具有一定的相關(guān)性，其他相應(yīng)觀測(cè)基本不一致。統(tǒng)計(jì)庫(kù)爾勒地傾斜年波滯后氣溫相位，并計(jì)算它們之間的相關(guān)系數(shù)，結(jié)果見(jiàn)表3。

表2 庫(kù)爾勒地傾斜各分量月波之間及與氣溫、氣壓月波相關(guān)系數(shù)

表3 庫(kù)爾勒地傾斜滯后溫度相位及與氣溫、氣壓相關(guān)系數(shù)(年頻段)

由圖6及表3，庫(kù)爾勒地傾斜的年波雖具有一定差異性，經(jīng)相位平移之后，均與氣溫具有較高的線性相關(guān)性。

4 結(jié)論與討論

庫(kù)爾勒地傾斜在日波、半日波頻段一致性較好，均與固體潮具有較高的線性相關(guān)性，這也充分說(shuō)明庫(kù)爾勒地傾斜固體潮記錄清晰可靠，日波、半日波頻段能較好地反映出固體潮的變化。

楊紹富[17]認(rèn)為氣壓對(duì)庫(kù)爾勒水管儀的影響既有對(duì)測(cè)點(diǎn)載荷的影響，也有通過(guò)觀測(cè)系統(tǒng)產(chǎn)生的影響，氣體流動(dòng)會(huì)在水管儀兩端的缽體產(chǎn)生氣壓差，從而影響到液面高度，對(duì)觀測(cè)數(shù)據(jù)產(chǎn)生影響，氣壓的影響始終存在。而本文中的研究結(jié)果表明庫(kù)爾勒地傾斜月波頻段的一致性較差，水管儀雖然與氣壓具有一定的相關(guān)性，但相關(guān)程度較低，僅垂直擺EW分量與氣壓月波具有一定的相關(guān)性。

復(fù)雜地形條件下地傾斜的年變化幅度受到溫度的影響更為顯著，山坡所受影響比近山腳的平地大1個(gè)數(shù)量級(jí)[18]，在假定半無(wú)限空間介質(zhì)模型下，地表溫度作周期變化時(shí)，地下深度y處溫度與地表溫度波動(dòng)的相位差為[19]

(2)

其中，φ為相位差，ω為周年變化角頻率，k為熱擴(kuò)散系數(shù)。上式表明，溫度波動(dòng)相位差與深度成正比。由圖1，垂直擺距洞口較近，水平擺與水管儀距洞口較遠(yuǎn)；根據(jù)山體特征，水平擺、水管儀上層覆蓋較厚，垂直擺上層覆蓋相對(duì)較薄。因此垂直擺相位滯后天數(shù)較水平擺與水管儀相位滯后天數(shù)少。庫(kù)爾勒地傾斜年波頻段經(jīng)相位平移之后，均與氣溫具有較高的線性相關(guān)性，這說(shuō)明庫(kù)爾勒地傾斜年周期變化可能受溫度的影響較大。

綜合上述分析可知，在庫(kù)爾勒地傾斜的數(shù)據(jù)跟蹤分析過(guò)程中，若提取到日波、半日波頻段的異常信息需要考慮剔除固體潮汐的影響，若提取到年變幅度及相位的異常變化時(shí)，需要考慮剔除溫度變化的影響。