李洪麗
(吉林大學(xué) 地球探測科學(xué)與技術(shù)學(xué)院,長春 130026)
微動方法沒有固定的震源,它是利用地球本身存在的由于人類活動和自然現(xiàn)象等帶來的微動信號來研究地殼淺層信息。該方法以其操作簡單、施工方便、探測成本低、對周圍環(huán)境無破壞、適用范圍廣等優(yōu)勢,被越來越多的國內(nèi)、外研究者采用。通過微動的頻譜特性來研究微動在時間和空間上的變化,進而獲得探測區(qū)域地殼淺層速度結(jié)構(gòu)。地殼淺層速度結(jié)構(gòu)能為活動斷裂、強地面運動模擬、地?zé)豳Y源勘查以及地震工程場地評估等地質(zhì)問題提供幫助[1]。高分辨S波速度結(jié)構(gòu)反演還能很好地反映地下介質(zhì)低速層的空間變化,為地球淺部圈層耦合性關(guān)系、強震孕震環(huán)境及區(qū)域動力學(xué)機制等研究提供重要參考[2]。隨著國內(nèi)、外關(guān)于地震學(xué)理論和方法研究的不斷發(fā)展,基于地震背景噪聲數(shù)據(jù)的微動方法在研究地殼淺層速度結(jié)構(gòu)成像、地球內(nèi)部波速結(jié)構(gòu)隨時間變化關(guān)系,以及測定地震水平位置等方面獲得了眾多有意義的研究成果。在地殼淺部數(shù)十米至數(shù)百米深度的波速結(jié)構(gòu)探測中,微動探測方法已經(jīng)在地質(zhì)調(diào)查、土木工程等實際應(yīng)用領(lǐng)域得到了大力推廣,并且取得了良好的應(yīng)用效果。
從微動信號中提取面波頻散曲線的方法主要有,空間自相關(guān)法(SPAC 法)、頻率-波數(shù)法(F-K 法)、H/V譜比法等,其中空間自相關(guān)法最為常用[3-4]。Aki[5]提出了空間自相關(guān)法的概念,研究了在平穩(wěn)隨機波動理論假設(shè)下一維和二維的空間自相關(guān)理論及應(yīng)用實例;何正勤等[6]采用空間自相關(guān)法,研究了地殼淺層1 km~2 km深度的橫波速度結(jié)構(gòu),通過與鉆井結(jié)果比對,驗證了該方法對于確定物性差異性大的地層界面是可靠的;徐佩芬等[7]基于自相關(guān)法提取地震臺站微震信號相速度頻散曲線,并反演得到地殼淺層S波速度結(jié)構(gòu),結(jié)果與人工地震探測結(jié)果較為吻合,驗證了SPAC法估算地殼S波速度結(jié)構(gòu)的可行性;李細兵等[8]利用空間自相關(guān)法獲得一個凹陷盆地基本速度結(jié)構(gòu)和各層深度分布范圍,結(jié)果顯示該區(qū)域的土層結(jié)構(gòu)接近人工地震勘探解釋結(jié)果;張寶龍等[1]利用密集臺陣接收的數(shù)據(jù),基于微動和面波層析成像方法,獲得了在五大連池火山區(qū)尾山火山錐地區(qū)地表至700 m深度的地殼淺層三維S波速度結(jié)構(gòu);謝朋等[9]利用微動探測方法獲取潛江市浩口中學(xué)地下介質(zhì)結(jié)構(gòu)分層,并與鉆孔結(jié)果對比,表明微動探測方法能較好地分辨地層結(jié)構(gòu);盛勇等[10]通過提取頻散曲線和視S波速度分層,描述了地層的速度結(jié)構(gòu)及地?zé)岙惓L卣鳎砻魑涌碧侥苡行У卦诘責(zé)豳Y源勘查中應(yīng)用。上述研究成果表明,微動法近年來在地質(zhì)調(diào)查實際中應(yīng)用已經(jīng)十分普遍,并且該方法也是行之有效的研究方法之一。
采用微動方法(空間自相關(guān)),對吉林市某地區(qū)地?zé)峥碧絽^(qū)域地殼淺層S波速度結(jié)構(gòu)進行研究,研究過程包括野外采集數(shù)據(jù),對采集到的原始地震數(shù)據(jù)進行瑞雷波相速度頻散曲線提取,再對提取得到的頻散曲線進行擬合反演,進而獲得研究區(qū)域地下1 000 m 深度的地殼淺層二維S波速度結(jié)構(gòu)剖面。野外觀測儀器為短周期地震儀(圖1)10臺,為了減小由于噪聲源方位分布不均勻而帶來的對面波頻散曲線可靠性的影響,野外觀測所用10臺短周期地震儀布置位置采用多重圓觀測系統(tǒng)(圖2)。每組如圖2所示布設(shè)10臺短周期地震儀,其中中心點為觀測目標(biāo)點,放一臺短周期地震儀,外圍距中心點的距離分別為50 m、100 m、200 m,共計觀測5個目標(biāo)點,形成一條觀測剖面,每個目標(biāo)點觀測時間不少于15 h,研究深度為地表以下1 000 m以內(nèi)地層。五個觀測目標(biāo)點,分別記作S4-1、S4-2、S4-3、S4-4、S1-8。
圖1 短周期地震儀
圖2 微動單點觀測系統(tǒng)(三重圓形臺陣)示意圖
微動方法主要目的是,從觀測的噪聲信號中提取瑞雷波相速度頻散曲線,并用瑞雷波相速度頻散反演地下S波速度結(jié)構(gòu),所以提取準(zhǔn)確的頻散曲線至關(guān)重要。本研究采用空間自相關(guān)(SPAC)方法,該方法是基于微動信號在時間和空間上平穩(wěn)隨機分布這一假設(shè)。在該假設(shè)條件下,對相距為r的兩個點的微動記錄作空間相關(guān)運算。Aki給出了方位平均后的空間自相關(guān)函,數(shù)可以用零階貝塞爾函數(shù)來表示:
ρ(r,f)=J0(2πfr/c(f))
(1)
式中:r為兩臺站間的距離;f為頻率;c(f)為瑞雷波相速度;ρ(r,f)為方位平均后的自相關(guān)系數(shù);J0為第一類零階貝塞爾函數(shù)[5]。式(1)適用于基階面波能量為主,利用垂向記錄分量進行空間自相關(guān)提取瑞雷面波頻散的情形。
空間自相關(guān)系數(shù)的計算大致分為時域(窄帶濾波法)和頻域兩大類。時域空間自相關(guān)系數(shù)求取的步驟為:①獲取記錄良好的原始記錄,對每段選取的記錄平均分成若干段,分段時要保證可以覆蓋足夠多的頻率信息;②對數(shù)據(jù)進行不同中心頻率的窄帶濾波處理,得到以濾波器中心頻率為自變量的近似單色波記錄;③將單頻記錄的距離為r的臺站對按式(2)計算空間自相關(guān)系數(shù),并將相同臺站對的結(jié)果進行平均;④循環(huán)步驟②和步驟③,計算出一系列頻率值所對應(yīng)臺站對的空間自相關(guān)系數(shù),對所得的空間自相關(guān)系數(shù)按式(3)兩兩取方位平均,得到用于計算頻散的平均自相關(guān)系數(shù)。
(2)
(3)
空間自相關(guān)系數(shù)的頻率域計算公式為式(4)。
(4)
平面上圓形排列臺陣中各臺站的位置一般用極坐標(biāo)形式表示。以參考點為圓心,θ為臺站方位角,S(r,θ,f)為參考點記錄與其他記錄點的互功率譜,S(0,f)和S(r,f)分別為參考點和與其距離r的另一點記錄的自功率譜;Re為取實部,確保空間自相關(guān)系數(shù)為實數(shù)。
頻域計算相對于時域計算的優(yōu)勢是去掉了多次窄帶濾波過程,只利用一次傅里葉變換后完全在頻率域內(nèi)進行計算,因此大大加快了運算速度,所以本研究是在頻率域進行計算的。
頻域空間自相關(guān)系數(shù)求取的步驟:①選取良好的記錄數(shù)據(jù),將數(shù)據(jù)分段;②將分段后的數(shù)據(jù)進行傅里葉變換得到頻率域數(shù)據(jù),在頻率域內(nèi)分別計算對應(yīng)數(shù)據(jù)段的自功率譜和互功率譜,然后將各段數(shù)據(jù)結(jié)果進行平均;③將平均后的自功率譜和交叉譜結(jié)果代入式(4)計算方位平均,便可得到方位平均的空間自相關(guān)系數(shù)。
野外數(shù)據(jù)采集過程中盡量避開有大量人為因素干擾的區(qū)域,經(jīng)過5個觀測點累計50臺短周期地震儀的記錄,完成5個觀測目標(biāo)點的微動信息數(shù)據(jù)采集工作。其中觀測目標(biāo)點S4-1的部分原始采集數(shù)據(jù)波形如圖 3 所示,可以看出各個地震儀記錄波形一致性較好,表明采集到的微動信息可靠。
圖3 S1點的10臺短周期地震儀接收的部分連續(xù)地震背景噪聲數(shù)據(jù)
在計算各觀測目標(biāo)點空間自相關(guān)系數(shù)前,先整理5個觀測目標(biāo)點的各地震儀接收到的原始數(shù)據(jù),再分別進行去均值、去傾斜、去尖脈沖等數(shù)據(jù)預(yù)處理工作。將每個短周期地震儀垂直分量噪聲數(shù)據(jù)按照時間進行分段處理,然后對所有分段數(shù)據(jù)進行傅里葉變換,按照用一臺短周期地震儀的傅里葉譜與另一個臺的傅里葉譜的共軛計算自相關(guān)系數(shù),再對自相關(guān)系數(shù)波形進行平滑處理,得到各觀測目標(biāo)點對應(yīng)的自相關(guān)系數(shù)曲線。數(shù)據(jù)處理過程中將信噪比較低、未能呈現(xiàn)第一類零階貝塞爾曲線形態(tài)的自相關(guān)系數(shù)曲線剔除。最后通過分析自相關(guān)系數(shù)波形的零點和極點,與第一類零階貝塞爾函數(shù)零點和極點進行擬合,得到五個觀測目標(biāo)點的面波相速度頻散曲線(圖4)。
圖4 瑞雷面波相速度頻散曲線
采用CPS330軟件分別對上述5個觀測目標(biāo)點的面波相速度頻散曲線進行線性反演,求取各觀測目標(biāo)點的一維速度模型[11]。再分別對各觀測目標(biāo)點的一維速度模型進行插值處理并合成,最終獲得該研究區(qū)域地表以下1 000 m以內(nèi)的淺表地下介質(zhì)二維S波速度剖面(圖5)。
圖5 研究區(qū)域地殼淺層二維S波速度剖面
S波速度對于地層內(nèi)的流體、溫度變化和介質(zhì)性質(zhì)十分敏感。S波速度結(jié)構(gòu)研究有助于分析研究區(qū)域地下介質(zhì)物理狀態(tài)。地下介質(zhì)巖性是影響地震波速的決定因素,介質(zhì)中存在的裂隙和孔隙度是基本影響因素。另外,地下介質(zhì)孔隙中的不同填充物、介質(zhì)的風(fēng)化和破碎帶以及埋藏深度都會對地層波速產(chǎn)生影響。結(jié)合吉林市周邊地區(qū)地層劃分表,從獲得的1 000 m深度的S波速度結(jié)構(gòu)剖面圖(圖5)中可以分析出:研究區(qū)域內(nèi)S波速度整體由地表至地下隨地層的深度增大呈增大趨勢,研究區(qū)域在深度上可分為三層,第一層速度較低,底界面深度在100 m~150 m左右,為第四紀(jì)和晚新生代(新近紀(jì))沉積地層;第二層S波速度在1 000 m/s左右,底界面深度在550 m左右,推測為早新生代(古近紀(jì))沉積地層;第三層地層速度約1 800 m/s,推測為中生代(侏羅紀(jì))沉積地層。由于不同巖性的地層速度不同,再者由于沉積年代以及沉積環(huán)境的差別,即使是同一巖性的地層速度也會存在一定的差異,所以依據(jù)單一S波速度參數(shù)所分層位可能具有一定的誤差。
另外,由于地震波速在氣體和液體中傳播速度低于在固體(巖石骨架)中的傳播速度,裂隙的存在會使地震波傳播速度減慢。由此可推測出在第三層中顯示的兩個明顯的低速區(qū)可能為含水區(qū)。
利用由10 臺短周期地震儀組成的短周期地震儀臺陣,在吉林某工區(qū)開展了連續(xù)地震背景噪聲觀測,共選取觀測目標(biāo)點5個,組成一個觀測剖面。利用微動方法首次獲得該區(qū)域地表至地下1 000 m深度的淺層S波速度結(jié)構(gòu)剖面。結(jié)果表明:地表至地下1 000 m之間S波速度整體隨深度具有線性增加的趨勢;研究區(qū)整體分為三層:由淺至深依次為第四紀(jì)和新近紀(jì)沉積層,古近紀(jì)沉積層和中生代(侏羅紀(jì))沉積層;在第三層中550 m~950 m區(qū)域存在兩個低速異常區(qū),推測為富含水區(qū)域。該微動觀測結(jié)果將為該區(qū)域今后開展地?zé)崽綔y、野外地質(zhì)調(diào)查等研究工作提供重要參考。