跨斷層輔助測(cè)線對(duì)非構(gòu)造因素去除的應(yīng)用探討

路珍 郭泉 李瑞莎 張希 唐紅濤 賈鵬 汪翠枝

摘要：將跨斷層場(chǎng)地的輔助測(cè)線應(yīng)用于跨斷層形變資料的數(shù)據(jù)處理中，引入譜減法對(duì)沿鮮水河斷裂帶的5個(gè)跨斷層場(chǎng)地進(jìn)行分析處理，試圖將輔助測(cè)線與跨斷層測(cè)線都包含的非構(gòu)造因素的影響從跨斷層數(shù)據(jù)中去除，從而突出與地震發(fā)生相關(guān)的斷層活動(dòng)信息。結(jié)果表明：利用譜減法處理過(guò)的跨斷層測(cè)線在多數(shù)場(chǎng)地都表現(xiàn)出了時(shí)頻成分上較明顯的異常，異常出現(xiàn)時(shí)間與周期特點(diǎn)都有較好的一致性。由鮮水河斷裂帶附近地震活動(dòng)性可得異常在時(shí)間上存在與后續(xù)地震活動(dòng)性增加相關(guān)的可能性。

關(guān)鍵詞：輔助測(cè)線；跨斷層測(cè)線；譜減法；非構(gòu)造因素；震前異常

中圖分類號(hào)：P315725文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2017）02-0193-10

0引言

已有的研究認(rèn)為地震發(fā)生的機(jī)理是沿?cái)鄬訋У牟环€(wěn)定粘滑過(guò)程（Brace，Byerlee，1966；Ruina，1983），以往震例總結(jié)和構(gòu)造物理實(shí)驗(yàn)揭示地震發(fā)生之前往往在應(yīng)力應(yīng)變、位移等方面都有跡象可尋（牛安福等，2012；李瑞莎等，2012；薄萬(wàn)舉，華彩虹，2001；李杰等，2011；馬瑾等，2012；卓燕群等，2013；季靈運(yùn)等，2015）?？鐢鄬訙y(cè)量是指利用布設(shè)在斷層帶兩盤的短水準(zhǔn)和短基線測(cè)量來(lái)獲得橫跨斷層帶兩測(cè)點(diǎn)的高差和位移隨時(shí)間的變化，可以較直接、較可靠地監(jiān)測(cè)斷層帶的活動(dòng)，是捕捉地震前斷層異?；顒?dòng)直接有力的方法之一。

近50年的跨斷層形變監(jiān)測(cè)不但積累了豐富的觀測(cè)資料和地震前兆異常，也在我國(guó)地震預(yù)報(bào)實(shí)踐中發(fā)揮了重要的作用，曾多次較為準(zhǔn)確地預(yù)測(cè)地震，如1975年2月4日海城MS73地震（朱鳳鳴，吳戈，1982）、1995年7月22日永登MS58地震（江在森等，1998）、2000年6月6日景泰MS59地震（王雙緒等，2001），還有近期的2013年4月20日蘆山MS70地震（蘇琴等，2014）。由此可見(jiàn)，對(duì)跨斷層形變數(shù)據(jù)中包含的與地震前兆信號(hào)有關(guān)的異常特征進(jìn)行識(shí)別、捕捉和跟蹤具有重要且實(shí)用的地震預(yù)報(bào)意義。

跨斷層形變觀測(cè)中不但包含斷層構(gòu)造運(yùn)動(dòng)的信息，同時(shí)也包含非構(gòu)造的環(huán)境影響等干擾信息，如氣象、地下水位、振動(dòng)、荷載變化、人為干擾等（樓關(guān)壽等，2010）。非構(gòu)造干擾信息與有用的斷層活動(dòng)信息一起被記錄下來(lái)，影響著對(duì)跨斷層數(shù)據(jù)的解釋，為構(gòu)造活動(dòng)的研究和前兆信息的確認(rèn)帶來(lái)了很多不確定因素。地震預(yù)報(bào)研究人員一直致力于對(duì)干擾因素的判定和提出，試圖提取出和孕震直接相關(guān)的構(gòu)造活動(dòng)信息（薄萬(wàn)舉，1997；施順英等，2007；李臘月等，2013a，b；張希等，2012；劉權(quán)威，1996；Moro，Zadro，1998；Matsumoto et al，2003）。然而，已有的對(duì)跨斷層數(shù)據(jù)的處理方法大多是基于跨斷層帶測(cè)線進(jìn)行分析的，而場(chǎng)地的輔助測(cè)線則一直被忽略。本文中的輔助測(cè)線是指為了測(cè)量精度（閉合差在規(guī)定范圍內(nèi)）或者為了對(duì)場(chǎng)地更好地控制而在斷層帶周邊設(shè)立的未跨過(guò)斷層帶的測(cè)線，其往往設(shè)在斷層帶的一盤，且在斷層帶附近。

由于輔助測(cè)線通常都在斷層帶附近，因此輔助測(cè)線和跨斷層測(cè)線都包含場(chǎng)地非構(gòu)造影響因素，其主要區(qū)別是跨斷層數(shù)據(jù)包含了斷層帶的活動(dòng)信息，而輔助測(cè)線理論上并沒(méi)有。因此，本文對(duì)跨斷層測(cè)線與輔助測(cè)線數(shù)據(jù)進(jìn)行分析處理，試圖用譜減法將二者共同部分非構(gòu)造因素的影響從中去除，并將所得結(jié)果與以往震例進(jìn)行分析對(duì)比，對(duì)所用分析方法進(jìn)行客觀評(píng)估。

1研究思路及方法

11數(shù)據(jù)分析方法[BT）]

本研究將語(yǔ)音通信中信號(hào)增強(qiáng)技術(shù)中所常用的譜減法引入跨斷層數(shù)據(jù)的分析處理中，其基本算法是將非目標(biāo)信號(hào)作為加性噪聲疊加在目標(biāo)信號(hào)上形成混合信號(hào)，在頻率域中從混合信號(hào)的功率譜中減去非目標(biāo)信號(hào)功率譜，從而得到較為純凈的目標(biāo)信號(hào)。該方法原理比較簡(jiǎn)單，實(shí)驗(yàn)結(jié)果也表明此法在語(yǔ)音信號(hào)中的增強(qiáng)效果很好，可以有效去除背景噪聲（黃斌，曾慶寧，2008；林琴等，2007）。本文所使用的跨斷層數(shù)據(jù)包含斷層活動(dòng)信息和非構(gòu)造影響因素，可看作是帶噪聲的混合信號(hào)。未跨斷層測(cè)線數(shù)據(jù)原理上不包含斷層運(yùn)動(dòng)信息，雖然可能會(huì)有一定的構(gòu)造意義，如可反映斷層一盤的掀斜運(yùn)動(dòng)，而掀斜運(yùn)動(dòng)可認(rèn)為是旋轉(zhuǎn)運(yùn)動(dòng)，其一般較微弱，效應(yīng)并不顯著（Stupazzini et al，2009；Trifunac，2009）。另外，跨斷層場(chǎng)地測(cè)線長(zhǎng)度通常較短，一般只有幾百米，能反映的范圍相當(dāng)有限。因此，基于上述考慮，我們認(rèn)為未跨斷層數(shù)據(jù)主要為非構(gòu)造因素，即可對(duì)應(yīng)為噪聲非目標(biāo)信號(hào)，可用于譜減法的操作。譜減法的原理如下所示（黃斌，曾慶寧，2008；林琴等，2007）。

具體數(shù)據(jù)處理步驟簡(jiǎn)述如下：

（1）對(duì)觀測(cè)資料進(jìn)行連續(xù)性檢驗(yàn)，剔除重復(fù)數(shù)據(jù)，用線性插值法填補(bǔ)缺失的數(shù)據(jù)，并使數(shù)據(jù)周期更加均勻，得到每年12個(gè)周期的連續(xù)性較好的數(shù)據(jù)。

（2）采用小波變換方法對(duì)數(shù)據(jù)進(jìn)行時(shí)頻分析，得到數(shù)據(jù)的時(shí)頻結(jié)果。數(shù)據(jù)處理是在Matlab的環(huán)境里編譯cmor3小波算法進(jìn)行的，具體算法如下（胡昌華等，1999）：

morl小波函數(shù)的定義為：[KH*1]

ψ（t）=Ce-t[HT4.5]2[HT7"]/2cos5t[JY]（8）[KH*1D]

經(jīng)伸縮和平移后得：[KH*1]

ψa，b（t）=[SX（]1[KF（][JB<1|]a[JB>1|][KF）][SX）]ψ[JB<2（][SX（]t-ba[SX）][JB>2）]（a，b∈R；a≠0）[JY]（9）[KH*2D]

其中：a為伸縮因子；b為平移因子。

對(duì)于任意的函數(shù)f（t）∈L2（R）〔L2（R）表示平方可積的實(shí)數(shù)空間〕的連續(xù)小波變換為：[KH*1]

Wf（a，b）[KG-*2]=[KG-*4][WB]〈f，ψa，b〉=

[DW][KG-*2][JB<1|]a[JB>1|]-1/2∫Rf（t）ψ[JB<2（][SX（]t-ba[SX）][JB>2）][TX-]dt[JY]（10）[KH*1]

（3）基于前述譜減法的原理，在Matlab環(huán)境中編譯算法，進(jìn)行相應(yīng)分析，得到數(shù)據(jù)的譜減法分析結(jié)果。[HJ2mm]

[BT2]12場(chǎng)地介紹

本文選取沿鮮水河斷裂帶且具有跨斷層測(cè)線和輔助測(cè)線的5個(gè)水準(zhǔn)場(chǎng)地作為研究對(duì)象，這些場(chǎng)地測(cè)量精度優(yōu)于10-6，具有數(shù)據(jù)質(zhì)量高、復(fù)測(cè)周期短、資料連續(xù)性好等特點(diǎn)。沿?cái)嗔褞нx取有利于基于斷裂構(gòu)造背景進(jìn)行數(shù)據(jù)分析。所選場(chǎng)地為侏倭、格簍壩子、虛墟、龍燈壩和老乾寧場(chǎng)地，如圖1所示，各場(chǎng)地測(cè)線布設(shè)如圖2所示。

2研究結(jié)果及分析

21原始數(shù)據(jù)分析[BT）]

為了更好地將各個(gè)場(chǎng)地間的數(shù)據(jù)進(jìn)行對(duì)比分析，統(tǒng)一選取1984—2016年的數(shù)據(jù)進(jìn)行分析。圖3列出了各個(gè)場(chǎng)地跨斷層測(cè)線與輔助測(cè)線的原始數(shù)據(jù)曲線及其時(shí)頻分析的結(jié)果。圖中藍(lán)色虛線為沿鮮水河斷裂帶200 km范圍內(nèi)MS≥50地震的發(fā)生時(shí)刻，黑色虛線為沿鮮水河斷裂帶300 km范圍內(nèi)的MS≥60地震的發(fā)生時(shí)刻。

（1）圖3a所示為侏倭水準(zhǔn)測(cè)量的結(jié)果，A←B和A←C測(cè)線都為跨斷層測(cè)線，C-B測(cè)線為輔助測(cè)線。對(duì)于侏倭場(chǎng)地，跨斷層測(cè)線對(duì)斷層的長(zhǎng)期變化趨勢(shì)有一個(gè)較好的反映，都表現(xiàn)出下盤相對(duì)上盤呈上升的趨勢(shì)。而輔助測(cè)線C-B則表現(xiàn)出了較大的不同，只有微弱的遞減趨勢(shì)，且呈鋸齒狀（圖3a-5）。跨斷層測(cè)線與輔助測(cè)線的時(shí)頻特征都表現(xiàn)為周期為100 d以上的能量占主要成分，1995年以后周期為30～50 d以上的能量成分有所增加。由于原始跨斷層曲線呈平滑上升趨勢(shì)，與地震事件的對(duì)應(yīng)并不好（圖3a-1，a-3）。但其時(shí)頻分析結(jié)果對(duì)地震的對(duì)應(yīng)程度有所上升，如白玉MS55、雅江MS60、汶川MS80和魯?shù)镸S65地震都不同程度的對(duì)應(yīng)時(shí)頻圖上能量的增加（圖3a-2，a-4）。

[JP2][KG2]（2）虛墟場(chǎng)地的跨斷層測(cè)線B←A與C←A（圖3b）表現(xiàn)為斷層下盤相對(duì)上盤呈下降的趨勢(shì)。輔助測(cè)線B-C隨時(shí)間則為緩慢增加的趨勢(shì)，同時(shí)伴隨鋸齒狀的跳躍（圖3b-5）。其跨斷層測(cè)線與輔助測(cè)線時(shí)頻特征主要表現(xiàn)為周期50～100 d以上的能量，后期為周期30～50 d以上的能量。虛墟場(chǎng)地原始跨斷層曲線整體表現(xiàn)為上下跳動(dòng)的特點(diǎn)，對(duì)部分地震有一定的反映，如馬邊MS57、蘆山MS70和魯?shù)镸S65地震前測(cè)線都有一定程度的向上跳躍，對(duì)應(yīng)拉張作用的增加（圖3b-1，b-3）。其時(shí)頻結(jié)果對(duì)白玉MS55和雅江MS60的對(duì)應(yīng)較原始曲線有所提升，表現(xiàn)為能量的增加（圖3b-2，b-4）。[JP]

[KG（0.15mm]（3）格簍壩子場(chǎng)地的跨斷層測(cè)線A←C（圖3c）表現(xiàn)為斷層下盤相對(duì)上盤先下降后上升的變化，其在2002年左右出現(xiàn)約17 mm的臺(tái)階突跳，2012—2013年底一直呈現(xiàn)上升趨勢(shì)，且增加速率逐漸增大，上升幅度約為25 mm。前者對(duì)應(yīng)新龍MS53地震，后者在上升過(guò)程中發(fā)生了蘆山MS70地震，在上升相對(duì)轉(zhuǎn)平一年多的時(shí)間內(nèi)發(fā)生了魯?shù)镸S65和康定MS63地震（圖3c-1）。其輔助測(cè)線先呈減小趨勢(shì)后趨勢(shì)變化不明顯，同樣一直伴隨鋸齒狀的跳躍（圖3c-3）。該跨斷層測(cè)線與輔助測(cè)線時(shí)頻特征大多為周期50～100 d以上的能量，后期為30～50 d的能量。相比而言，其時(shí)頻圖結(jié)果未及原始曲線的映震效果（圖3c-2）。[KG）]

（4）龍燈壩與老乾寧場(chǎng)地的跨斷層測(cè)線與輔助測(cè)線如圖3d和e所示，其跨斷層測(cè)線大多表現(xiàn)為類周期性的突跳波動(dòng)特點(diǎn)，有微弱的上升趨勢(shì)（圖3d-1，圖3e-3）或者趨勢(shì)變化不明顯（圖3d-3，圖3e-1），在部分地震發(fā)生前有較大幅度向下突跳與之對(duì)應(yīng)，如白玉MS55、寧蒗MS62、雅江MS60等（圖3d-1，d-3，3e-1，e-3）。龍燈壩場(chǎng)地輔助測(cè)線B-D表現(xiàn)出微弱的上升趨勢(shì)，并一直伴隨有尖點(diǎn)突跳的異常（圖3d-5）。老乾寧場(chǎng)地的輔助測(cè)線5-1則表現(xiàn)出明顯的上升趨勢(shì)（圖3e-5）。龍燈壩與老乾寧2個(gè)場(chǎng)地的時(shí)頻分析也主要表現(xiàn)為前期周期50～100 d以上的信號(hào)為主要成分，后期主要為周期30～50 d以上的成分，其時(shí)頻分析結(jié)果也在一定程度對(duì)應(yīng)部分地震（圖3d-2，d-4，3e-2，e-4）。

根據(jù)上述分析，發(fā)現(xiàn)鮮水河斷裂帶上的跨斷層數(shù)據(jù)對(duì)于周圍一定范圍內(nèi)的5～6級(jí)地震具有一定程度的映震能力。但由于跨斷層原始測(cè)線大多表現(xiàn)為隨整個(gè)時(shí)間序列的上下波動(dòng)，所包含的異常較多，其中與地震有關(guān)的異常易被淹沒(méi)其中而難以區(qū)分，導(dǎo)致異常和地震的發(fā)生并不是一一對(duì)應(yīng)的，因此從眾多異常中確認(rèn)出明顯的、不同于其他的地震異常具有一定困難。因此，從跨斷層原始數(shù)據(jù)中突出地震異常的工作則顯得比較必要。

[BT2]22譜減法分析結(jié)果

本文試圖利用譜減法從跨斷層測(cè)線中扣除非構(gòu)造影響因素，由于譜減法的原理是從帶噪信號(hào)的功率譜中減去整體噪聲功率譜，因此結(jié)果中所去除的非構(gòu)造因素不是如氣象、地下水位、振動(dòng)、荷載變化、人為干擾等某一項(xiàng)干擾因素的影響，而是一個(gè)場(chǎng)地綜合環(huán)境的影響，是上述可能影響因素的綜合結(jié)果。分析結(jié)果見(jiàn)圖4，從圖中可以看出，譜減法的結(jié)果與原始數(shù)據(jù)特征有所不同，顯示對(duì)于侏倭、格簍壩子和老乾寧3個(gè)場(chǎng)地都出現(xiàn)了明顯的異常。具體表現(xiàn)為：

（1）侏倭場(chǎng)地A←B測(cè)線和A←C測(cè)線的譜減法結(jié)果都出現(xiàn)了明顯的能量的增強(qiáng)（圖4a-1，a-2），能量周期約為30～100 d，在50～100 d之間占主要成分，持續(xù)時(shí)間約為2001—2003年。

（2）格簍壩子場(chǎng)地譜減法結(jié)果也表現(xiàn)出了明顯的能量增強(qiáng)（圖4c），分別在40 d和100 d的2個(gè)周期成分處增強(qiáng)較明顯，持續(xù)時(shí)間大致為2000年下半年至2003年中，這個(gè)時(shí)間段與A←C測(cè)線在2002年左右出現(xiàn)的階躍異常在時(shí)間上有較好的一致（圖3c-1）。另外，對(duì)于格簍壩子場(chǎng)地，能量的增強(qiáng)還出現(xiàn)在2012年左右，主要為30～60 d周期的能量成分，持續(xù)時(shí)間大致為從2010年下半年至2012年下半年，這與其原始測(cè)線在2012—2014年的持續(xù)上升有一定的對(duì)應(yīng)。

（3）譜減法得到的老乾寧場(chǎng)地2條跨斷層測(cè)線1←3和5←3的異常周期都在50～100 d左右，持續(xù)時(shí)間大致為2001—2003年（圖4e-1，e-2）。而虛墟和龍燈壩2個(gè)場(chǎng)地的譜減法處理結(jié)果則較難從中指出明顯的、區(qū)別于其他的異?，F(xiàn)象。

總的來(lái)說(shuō)，通過(guò)譜減法得到的沿鮮水河斷裂帶的侏倭、格簍壩子和老乾寧3個(gè)跨斷層場(chǎng)地所出現(xiàn)異常的周期大致都在50～100 d，且異常出現(xiàn)的時(shí)間也都在2002年左右，表現(xiàn)出3個(gè)場(chǎng)地的異常在周期成分和出現(xiàn)時(shí)間上都有很好的一致性。雖然虛墟場(chǎng)地和龍燈壩場(chǎng)地并沒(méi)有出現(xiàn)如上述3個(gè)場(chǎng)地明顯的能量增強(qiáng)的異常，但是在2002年附近仍存在微弱的能量亮斑（圖4b、d）。

上述結(jié)果說(shuō)明，與原始數(shù)據(jù)易將震前異常信號(hào)淹沒(méi)在背景噪聲中不同，利用譜減法處理過(guò)的跨斷層測(cè)線在時(shí)頻成分上表現(xiàn)出了較明顯的異常，且這種異常出現(xiàn)在所研究的多數(shù)場(chǎng)地中，同時(shí)周期和時(shí)間特點(diǎn)都較一致。可見(jiàn)，利用譜減法去除非構(gòu)造因素后的結(jié)果在異常反映上較原始數(shù)據(jù)有一定程度的提高。

[BT2]23譜減法異常分析

[KG（0.15mm]由于所分析的跨斷層場(chǎng)地都是沿鮮水河斷裂帶選取的，因此在分析異常與地震發(fā)生關(guān)系的過(guò)程中，首先考慮鮮水河斷裂帶的地震活動(dòng)情況。圖5匯總了1983年以來(lái)沿鮮水河斷裂帶200 km范圍內(nèi)的MS50～60地震（灰實(shí)線）和300 km范圍內(nèi)的MS≥60地震（紅實(shí)線）。從圖中可以看出，200 km范圍內(nèi)的MS50～60地震有2個(gè)密集發(fā)生的時(shí)段，分別為1986—1996年和2007—2015年。1997—2007年有一個(gè)10年左右的相對(duì)小的平靜，期間只發(fā)生了2次5～6級(jí)地震。300 km范圍內(nèi)的MS≥60地震，集中出現(xiàn)于1989年，發(fā)生了巴塘MS66和小金北MS65等地震，隨后有一個(gè)近20年的MS≥60地震平靜期，期間只發(fā)生了1998年寧蒗MS62和2001年雅江MS60地震，在直到2008年這一平靜被打破，相繼發(fā)生了汶川MS80、玉樹(shù)MS71、蘆山MS70、魯?shù)镸S65、康定MS63等地震?？梢?jiàn)，無(wú)論是地震頻次還是地震量級(jí)，2008年以后都有一個(gè)大幅的增加。本文利用譜減法對(duì)跨斷層數(shù)據(jù)分析所得的異常現(xiàn)象出現(xiàn)在200 km范圍內(nèi)的5～6級(jí)地震的小平靜期內(nèi)，同時(shí)也處于2008年打破鮮水河斷裂帶300 km范圍內(nèi)近20年MS≥60地震平靜之前（圖5中深灰色陰影所示）。所得異常在時(shí)間上存在其與后續(xù)地震活動(dòng)增加、能量增強(qiáng)有關(guān)的可能性。因此，利用譜減法扣除非構(gòu)造因素后，鮮水河斷裂帶跨斷層測(cè)線時(shí)頻信息中表現(xiàn)出明顯且一致的異常頻率信息可能與斷層自身構(gòu)造作用的孕震信息有關(guān)，從而引起后續(xù)地震活動(dòng)性的增強(qiáng)。

3討論與結(jié)論

[KG（0.15mm]目前對(duì)跨斷層數(shù)據(jù)提取和提煉突出其構(gòu)造活動(dòng)特點(diǎn)較常見(jiàn)的方法有灰色關(guān)聯(lián)度（張希等，2012）和主成分分析法（李臘月等，2013a，b；方穎等，2015）。如張希等（2012）基于灰色關(guān)聯(lián)度分析方法提煉了反映鮮水河斷裂構(gòu)造變形特征方式的新綜合指標(biāo)，研究了其與斷裂附近強(qiáng)震孕育發(fā)生的可能關(guān)系，得出經(jīng)過(guò)分析處理可能觀測(cè)到昆侖山口西81級(jí)地震的影響、汶川80級(jí)地震前應(yīng)力場(chǎng)強(qiáng)化或擾動(dòng)影響、玉樹(shù)71級(jí)地震后短期影響。方穎等（2015）利用主成分分析方法對(duì)鮮水河斷裂西北段的跨斷層形變資料進(jìn)行了分析，得到第二主成分的時(shí)間系數(shù)在2001年底出現(xiàn)了大尺度的轉(zhuǎn)折，認(rèn)為可能是昆侖山81級(jí)大地震的影響。另外，其同樣提出在2008年之前，鮮水河斷裂帶滑動(dòng)出現(xiàn)了低頻信號(hào)的增強(qiáng)，可能與后續(xù)地震的發(fā)生有很強(qiáng)的相關(guān)性（方穎等，2015）。本文利用譜減法對(duì)鮮水河斷裂帶進(jìn)行研究，得出其在2002年左右表現(xiàn)出明顯的異常現(xiàn)象，這可能與2008年后續(xù)大地震的孕育有關(guān)，同時(shí)也不排除與昆侖山口西81級(jí)遠(yuǎn)震間的關(guān)系。可見(jiàn)，本文結(jié)果與前人相關(guān)工作具有相當(dāng)程度的一致性。[KG）]

前人已有方法與本文方法都是試圖用某種方法在原始數(shù)據(jù)的基礎(chǔ)上去除相關(guān)干擾信息，最大限度提煉出與斷層構(gòu)造活動(dòng)特征直接相關(guān)的信息。但前人工作主要都是利用跨斷層測(cè)線數(shù)據(jù)進(jìn)行分析的，而本研究引入了未跨斷層的數(shù)據(jù)，對(duì)觀測(cè)數(shù)據(jù)進(jìn)行了充分應(yīng)用。

雖然我們利用譜減法對(duì)跨斷層數(shù)據(jù)進(jìn)行處理得到了較明顯一致且與地震相關(guān)性較高的異常現(xiàn)象，但目前對(duì)譜減法在跨斷層數(shù)據(jù)處理中還存在一些弊端。由于譜減法得到的時(shí)頻信息在進(jìn)行逆變換時(shí)會(huì)造成相位信息的缺失（黃斌，曾慶寧，2008；林琴等，2007），同時(shí)又受到跨斷層資料采樣率的限制，因此目前還較難將譜減法所得到的結(jié)果還原為時(shí)域上的信息。與已有方法相比，時(shí)域信息的缺失會(huì)限制對(duì)斷層活動(dòng)信息的獲取和對(duì)異常信息的解釋，使譜減法的應(yīng)用目前只停留在更好地突出異常，而限制了對(duì)異常的動(dòng)力學(xué)來(lái)源進(jìn)行判斷。因此，將譜減法的結(jié)果轉(zhuǎn)化至?xí)r域，獲得斷層帶隨時(shí)間變化的運(yùn)動(dòng)信息是以后研究的努力方向。

本文將以往都忽略的輔助測(cè)線引入到了跨斷層數(shù)據(jù)處理中，結(jié)合輔助測(cè)線并利用譜減法將跨斷層測(cè)線數(shù)據(jù)中的非構(gòu)造因素去除，從而突出與未來(lái)地震發(fā)生相關(guān)的構(gòu)造異常信息。研究結(jié)果顯示多數(shù)臺(tái)站、相同時(shí)間段都出現(xiàn)了類似的異?，F(xiàn)象，其表明譜減法對(duì)于去除跨斷層非構(gòu)造因素在頻域上是可行的，為更進(jìn)一步工作打下了一個(gè)基礎(chǔ)。盡管目前對(duì)異常的解釋受到了限制，譜減法所得到的異常結(jié)果對(duì)于強(qiáng)震跟蹤與地震預(yù)測(cè)預(yù)報(bào)還是有積極意義的。[KH1*2D]

參考文獻(xiàn)：

薄萬(wàn)舉，華彩虹2001地形變強(qiáng)震前兆指標(biāo)探討[J].地震，21（1）：25-32

薄萬(wàn)舉1997信息流合成方法在震情分析研究中的應(yīng)用[J].西北地震學(xué)報(bào)，19（2）：42-48

方穎，張晶，江在森，等2015用跨斷層形變資料分析鮮水河斷裂西北段的運(yùn)動(dòng)特征[J].地球物理學(xué)報(bào)，58（5）：1645-1653

胡昌華，張軍波，夏軍，等1999基于MATLAB的系統(tǒng)分析與設(shè)計(jì)——小波分析[M].西安：先電子科技大學(xué)出版社

黃斌，曾慶寧2008基于譜減法語(yǔ)音增強(qiáng)研究[J].大眾科技，（10）：22-24，26

季靈運(yùn)，劉立煒，郝明2015利用InSAR技術(shù)研究滇西南鎮(zhèn)康—永德地區(qū)現(xiàn)今地殼形變特征[J].地震研究，38（1）：84-89

江在森，祝意青，王慶良，等1998永登58級(jí)地震孕育發(fā)生過(guò)程中的斷層形變與重力場(chǎng)動(dòng)態(tài)圖象特征[J].地震學(xué)報(bào)，20（3）：41-48

李杰，殷海濤，吳晨，等2011沂沭斷裂帶跨斷層形變異常及其解釋[J].大地測(cè)量與地球動(dòng)力學(xué)，31（4）：18-21

李臘月，邢成起，武安緒，等2013a用主成分分析方法綜合分析首都圈跨斷層形變資料[J].中國(guó)地震，29（2）：246-255

李臘月，邢成起，武安緒，等2013b張北62級(jí)地震前后斷層活動(dòng)時(shí)空演化特征分析[J].大地測(cè)量與地球動(dòng)力學(xué)，33（1）：29-34

李瑞莎，張希，唐紅濤，等2012玉門、民樂(lè)等地震前跨斷層形變異常及與定點(diǎn)前兆觀測(cè)的相互印證[J].西北地震學(xué)報(bào)，34（4）：388-392

林琴，張道信，吳小培2007一種基于改進(jìn)譜減法的語(yǔ)音去噪新方法[J].計(jì)算機(jī)技術(shù)與發(fā)展，17（7）：63-66

劉權(quán)威1996地殼形變干擾因素有效性的判別[J].測(cè)繪學(xué)報(bào)，25（2）：16-21，23-24

樓關(guān)壽，周偉，金鵬，等2010跨斷層形變觀測(cè)干擾因素的調(diào)查[J].大地測(cè)量與地球動(dòng)力學(xué)，30（增刊2）：68-74

馬瑾，SHERMAN S I，郭彥雙2012地震前亞失穩(wěn)應(yīng)力狀態(tài)的識(shí)別——以5°拐折斷層變形溫度場(chǎng)演化的實(shí)驗(yàn)為例[J].中國(guó)科學(xué)：地球科學(xué)，42（5）：633-645

牛安福，張凌空，閆偉，等2012汶川地震近震源區(qū)地形變短期前兆現(xiàn)象的解析[J].地震，32（2）：52-63

施順英，張燕，吳云，等2007基于跨斷層形變異常預(yù)測(cè)云南地震的試驗(yàn)[J].大地測(cè)量與地球動(dòng)力學(xué)，27（5）：82-87

蘇琴，楊永林，鄭兵，等20144·20蘆山70級(jí)地震預(yù)測(cè)思路及過(guò)程回顧[J].地震地質(zhì)，36（4）：1077-1093

王雙緒，江在森，陳文勝，等2001景泰59級(jí)地震的斷層形變異常及中短期預(yù)報(bào)[J].地震學(xué)報(bào)，23（2）：151-158

張希，唐紅濤，李瑞莎，等2012基于灰色關(guān)聯(lián)度的鮮水河斷裂活動(dòng)特性與大震關(guān)系研究[J].地震研究，35（4）：500-505

朱鳳鳴，吳戈1982一九七五年海城地震[M].北京：地震出版社

卓燕群，郭彥雙，汲云濤，等2013平直走滑斷層亞失穩(wěn)狀態(tài)的位移協(xié)同化特征——基于數(shù)字圖像相關(guān)方法的實(shí)驗(yàn)研究[J].中國(guó)科學(xué)：地球科學(xué)，43（10）：1643-1650

BRACE W，BYERLEE J1966Stick-slip as a mechanism for earthquakes[J].Science，153（3739）：990-992

MORO G D，ZADRO M1998Subsurface deformations induced by rainfall and atmospheric pressure：tilt/strain measurements in the NE-Italy seismic area[J].Earth and Planetary Science Letters，164（1-2）：193-203

MATSUMOTO N，KITAGAWA G，ROELOFFS E A2003Hydrological response to earthquakes in the Haibara well，central Japan – IGroundwater level changes revealed using state space decomposition of atmospheric pressure，rainfall and tidal responses[J].Geophysical Journal International，155（3）：885-898

RUINA A1983Slip instability and state variable friction laws[J].J geophys Res，88（B12）：359-310

STUPAZZINI M，DE LA PUENTE J，SMERZINI C，et al2009Study of Rotational Ground Motion in the Near-Field Region[J].Bulletin of the Seismological Society of America，99（2B）：1271-1286

TRIFUNAC M D2009Review：Rotations in Structural Response[J].Bulletin of the Seismological Society of America，99（2B）：968-979[ZK）][HJ][FL）][LM]

[STHZ][WT4HZ][JZ]A Tentative Discussion on the Application of Aided Measuring[JZ]Line to Remove the Non-tectonic Factors from[JZ]Cross-fault Deformation Data

[WT5B1][STBZ][JZ]LU Zhen1，GUO Quan2，LI Ruisha1，ZHANG Xi1，TANG Hongtao1，JIA Peng1，WANG Cuizhi1

[WT5"B1X][JZ]（1 The Second Monitoring and Application Center，CEA，Xian 710054，Shaanxi，China）[JZ]

（2 The Institute of Crustal Dynamics，CEA，Beijing 100085，China）

[WT5HZ][JZ]Abstract[WTB1]

This paper applied the aided measuring line of cross-fault site to the cross-fault deformation data processing In order to remove the non-tectonic influence from cross-fault deformation data and highlight the information of fault activities which is related to the earthquake，we used spectral subtraction to analyze the cross-fault deformation data which are obtained along Xianshuihe fault zone and contain aided measuring line By the processing of spectral subtraction，the results show that the cross-fault deformation data from most sites are all have obvious anomalies on time-frequency analysis And the occurrence and the period of those anomalies are all consistent Through analyzing the seismic activity along the Xianshuihe fault zone，the occurrence time of anomalies support the possibility that the anomalies are related to the increased seismic activities

[WTHZ]Keywords：[WTB1]aided measuring line；cross-fault measuring line；spectral subtraction method；non-tectonic factor； anomaly before earthquake