用微振動(dòng)勘探方法探測(cè)隱伏斷裂構(gòu)造

范小平 李細(xì)兵 趙啟光 宋 浩 戴 波 張帥帥

1） 中國(guó)南京210014江蘇省地震局 2） 中國(guó)長(zhǎng)春130061吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院

?

用微振動(dòng)勘探方法探測(cè)隱伏斷裂構(gòu)造

1） 中國(guó)南京210014江蘇省地震局 2） 中國(guó)長(zhǎng)春130061吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院

介紹了運(yùn)用微振動(dòng)勘探方法探測(cè)隱伏斷裂構(gòu)造的方法和思路， 并對(duì)郯廬斷裂帶兩條分支斷裂進(jìn)行了試驗(yàn)性勘探． 通過與人工地震勘探結(jié)果對(duì)比分析， 發(fā)現(xiàn)微振動(dòng)勘探方法對(duì)隱伏斷裂構(gòu)造具有較好的空間定位能力， 其對(duì)巖性差異面形成的斷裂定位誤差與人工地震勘探結(jié)果相當(dāng)； 又因其具有較大的勘探深度， 因而可以有效揭示勘探區(qū)域下方地層及斷裂的空間分布特征及其接觸關(guān)系， 這對(duì)于識(shí)別和確定斷裂構(gòu)造的主破裂面具有重要意義. 但微振動(dòng)勘探方法對(duì)地層或斷裂構(gòu)造細(xì)節(jié)特征的揭示能力較差． 因此， 在實(shí)際地層或斷裂勘探中， 宜采用人工地震勘探方法與微振動(dòng)勘探方法相結(jié)合的方式確定斷裂的深淺構(gòu)造關(guān)系， 可以用微振動(dòng)勘探方法進(jìn)行普查， 在確定斷裂初步位置及空間展布特征后， 再選用人工地震勘探方法進(jìn)行斷裂精確定位．

微振動(dòng)方法 隱伏斷裂構(gòu)造 淺層人工地震

引言

探測(cè)地下地質(zhì)構(gòu)造的手段主要是地震勘探方法， 而對(duì)隱伏斷裂的定位則主要以地震反射波方法為主． 反射波方法不僅具有較高的勘探精度， 而且最終的疊加剖面有助于判定斷層的存在與發(fā)育形態(tài)， 該方法已廣泛應(yīng)用于城市活斷層探測(cè)工作中（劉保金等， 2002）． 對(duì)于斷點(diǎn)埋藏較深或者復(fù)雜環(huán)境下進(jìn)行的斷裂定位， 反射波方法仍面臨不少困境， 如勘探深度有限、 有效波信噪比較低． 反射波方法勘探效果不僅與地震波的激發(fā)與接收條件等因素有關(guān)， 還與地下介質(zhì)物性成分和結(jié)構(gòu)有關(guān)（劉保金等， 2002； 潘紀(jì)順等， 2002）． 在炸藥震源限制使用的情況下， 雖可通過增加可控震源的數(shù)量來獲取較強(qiáng)能量的激發(fā)， 但從經(jīng)濟(jì)效益和施工效率來講， 反射波方法仍不是優(yōu)選的方案． 而得以重新發(fā)展的微振動(dòng)勘探方法（microtremor survey method， 簡(jiǎn)稱為MSM）則可以彌補(bǔ)人工地震勘探方法的欠缺．

微振動(dòng)方法是利用臺(tái)陣觀測(cè)技術(shù)采集地面微振動(dòng)信號(hào)， 采用數(shù)據(jù)處理與分析技術(shù)提取面波（瑞雷波）頻散信息， 再通過瑞雷波反演技術(shù)獲取地下介質(zhì)S波速度結(jié)構(gòu)的地球物理勘探方法（Okada, Suto, 2003; Okada, 2006; 徐佩芬等, 2013）． 微振動(dòng)信號(hào)是一種在時(shí)間域、 空間域都極為不規(guī)則的震動(dòng)現(xiàn)象， 既包含了諸如海浪、 風(fēng)速、 潮汐等自然現(xiàn)象產(chǎn)生的震動(dòng)， 也包含了人類活動(dòng)產(chǎn)生的震動(dòng)； 微振動(dòng)信號(hào)中既包含體波， 也包含面波， 其中面波成分占絕對(duì)優(yōu)勢(shì)． 微振動(dòng)勘探方法則利用面波（瑞雷波）信息來反演地下介質(zhì)結(jié)構(gòu)（徐佩芬等, 2009， 2013）．

微振動(dòng)勘探方法已廣泛應(yīng)用于盆地結(jié)構(gòu)調(diào)查（何正勤等， 2007）、 工程地質(zhì)勘察（徐佩芬等， 2012）、 場(chǎng)地穩(wěn)定性評(píng)價(jià)（陶夏新等, 2001）等多個(gè)領(lǐng)域， 并且在地?zé)豳Y源調(diào)查（徐佩芬等， 2013）、 煤礦采區(qū)構(gòu)造及采空區(qū)探測(cè)（徐佩芬等， 2013）、 城市地質(zhì)調(diào)查（王振東, 1986； 許建聰, 孫紅月, 2004； 師黎靜等， 2006; 董連成等， 2008； 王偉君等, 2009； 于凱等, 2011）等多個(gè)領(lǐng)域取得了應(yīng)用性成果． 其探測(cè)深度能從近地表到地下2—3 km （徐佩芬等， 2013）．

到目前為止， 利用微振動(dòng)勘探方法探測(cè)城市區(qū)隱伏斷裂構(gòu)造的研究和應(yīng)用并不多見. 為進(jìn)一步探討微振動(dòng)勘探方法探測(cè)城市區(qū)隱伏斷裂構(gòu)造的可能性及可靠性， 進(jìn)一步認(rèn)識(shí)該方法在隱伏斷裂構(gòu)造探測(cè)中的優(yōu)勢(shì)及不足， 本文利用微振動(dòng)勘探方法， 針對(duì)郯廬斷裂帶兩條分支斷裂即王莊—蘇圩斷裂（F1）和大官莊—雙莊斷裂（F2）進(jìn)行試驗(yàn)性勘探， 并將微振動(dòng)勘探結(jié)果與人工地震勘探結(jié)果進(jìn)行對(duì)比研究．

1 微振動(dòng)方法的基本原理

1.1 頻散曲線計(jì)算

微振動(dòng)是一種由體波（P波和S波）和面波（瑞雷波和勒夫波）組成的復(fù)雜振動(dòng)， 并且面波的能量占總能量的70%以上（Toks?z, Lacoss, 1968）． 盡管微振動(dòng)信號(hào)的振幅和形態(tài)隨時(shí)間變化， 但在一定時(shí)空范圍內(nèi)具有統(tǒng)計(jì)穩(wěn)定性， 可用時(shí)間與空間上的平穩(wěn)隨機(jī)過程描述（Aki, 1957）． 微振動(dòng)勘探方法就是以平穩(wěn)隨機(jī)過程為理論依據(jù)， 從微振動(dòng)信號(hào)中提取面波（瑞雷波）頻散曲線， 通過對(duì)頻散曲線反演獲取地下介質(zhì)的橫波速度結(jié)構(gòu)． 頻散曲線的穩(wěn)定性對(duì)探測(cè)結(jié)果具有決定性的影響， 且微振動(dòng)信號(hào)中既包含了長(zhǎng)周期信號(hào)， 也包含了短周期信號(hào)， 為了最大限度提取微振動(dòng)信號(hào)中的有用信息， 同時(shí)保證頻散曲線的穩(wěn)定性， 本文綜合空間自相關(guān)（spatial autocorrelation, 簡(jiǎn)稱為SPAC）法和頻率-波數(shù)法來提取微振動(dòng)觀測(cè)頻散曲線．

Aki（1957）提出的空間自相關(guān)法， 把復(fù)雜的地震波場(chǎng)作為時(shí)間與空間的穩(wěn)態(tài)隨機(jī)過程， 確定地震波的時(shí)間與空間頻譜關(guān)系， 把空間自相關(guān)系數(shù)ρ（r,ω）作為相速度c（ω）和頻率ω的函數(shù)， 即

（1）

式中， J0（x）是x的第一類零階貝塞爾函數(shù)，r為觀測(cè)半徑. 對(duì)不同r處的地震記錄進(jìn)行窄帶濾波， 獲取頻率ω對(duì)應(yīng)的空間自相關(guān)系數(shù)ρ（r,ω）， 再帶入式（1）便可求取不同ω相對(duì)應(yīng)的c（ω）． 經(jīng)多次對(duì)修正空間自相關(guān)（modified spatial autocorrelation, 簡(jiǎn)稱為MSPAC） （Bettigetal, 2001）、 擴(kuò)展空間自相關(guān)（extended spatial autocorrelation, 簡(jiǎn)稱為ESPAC） （Ling, Okada, 1993）及無中心點(diǎn)觀測(cè)排列（centerless circular array, 簡(jiǎn)稱為CAA） （Choetal, 2006）等方法的改進(jìn)， SPAC方法已不再局限于圓形臺(tái)陣幾何排列方式， 而是使用更為靈活的臺(tái)陣布設(shè)來適應(yīng)觀測(cè)環(huán)境． 本文采用修正空間自相關(guān)方法來提取頻散曲線．

頻率-波數(shù)法（f-k）（Capon, 1969）即將微振動(dòng)信號(hào)從時(shí)間-空間域轉(zhuǎn)換到頻率-波數(shù)域進(jìn)行分析， 利用波數(shù)矢量的方向性， 從頻率-波數(shù)譜圖中推算出入射波的方向和瑞雷波相速度的大?。?Capon（1969）提出了高分辨率集束（high resolution beamform, 簡(jiǎn)稱為HSBF）方法， 該方法用傳感器臺(tái)陣的互功率譜密度來估計(jì)f-k譜， 用相干來完善之．f-k譜可寫為

（2）

1.2 微振動(dòng)剖面形成

傳統(tǒng)的單點(diǎn)微振動(dòng)測(cè)試， 即利用實(shí)測(cè)瑞雷波頻散曲線反演得到臺(tái)陣中心點(diǎn)下方介質(zhì)的一維S波速度結(jié)構(gòu)． 但對(duì)于構(gòu)造探測(cè)而言， 我們往往需要了解剖面上巖性的相對(duì)變化而無需反演S波速度的絕對(duì)值． 所以， 在獲得單點(diǎn)實(shí)測(cè)相速度頻散曲線后， 利用式（3）將其轉(zhuǎn)換為S波視速度（Vx）隨深度d的變化（Vx-d）曲線， 再對(duì)各測(cè)點(diǎn)的Vx-d曲線進(jìn)行橫向插值、 平滑計(jì)算， 最終獲得二維視速度剖面（徐佩芬等, 2009， 2013）.

（3）

式中:vr為瑞雷波速度;ti為周期， 下標(biāo)i為頻點(diǎn)序數(shù)．Vx雖具有速度量綱， 但與具有明確物理含義的相速度Vr與S波速度VS具有質(zhì)的差異.Vx避免了反演過程中的初始模型設(shè)置、 反演結(jié)果選擇判斷等人為因素的影響. 微振動(dòng)剖面能客觀地、 有效地反映地層巖性變化， 是地質(zhì)解釋的主要依據(jù)． 該方法已成功應(yīng)用于實(shí)際工程勘察中（徐佩芬等， 2013）．

2 數(shù)據(jù)觀測(cè)及處理

2.1 數(shù)據(jù)觀測(cè)

微振動(dòng)觀測(cè)臺(tái)陣的半徑稱為觀測(cè)半徑r, 其決定探測(cè)深度H． 通常情況下， 探測(cè)深度H=（3—5）r（Okada, Suto, 2003； 徐佩芬等， 2013）． 為了滿足勘探深度要求以及數(shù)據(jù)處理的需求， 實(shí)際觀測(cè)中采用多個(gè)觀測(cè)臺(tái)陣組成的多重觀測(cè)系統(tǒng)進(jìn)行組合觀測(cè)． 其中三重圓形臺(tái)陣觀測(cè)系統(tǒng)由10個(gè)觀測(cè)點(diǎn)組成， 1個(gè)置于圓心， 其余9個(gè)分別置于內(nèi)接三角形的頂點(diǎn)（圖1）．

圖1 三重圓形微振動(dòng)觀測(cè)系統(tǒng)示意圖

圖2為微振動(dòng)觀測(cè)點(diǎn)、 斷裂構(gòu)造及地震勘探測(cè)線分布圖． 圖中斷裂構(gòu)造引用了“宿遷市城市活斷層探測(cè)與地震危險(xiǎn)性評(píng)價(jià)”項(xiàng)目對(duì)隱伏斷裂構(gòu)造的勘探成果， 其斷裂構(gòu)造的空間位置精度較高． 圖2中黑色三角形代表微振觀測(cè)臺(tái)陣中心觀測(cè)點(diǎn)的位置． 觀測(cè)點(diǎn)近東西向分布， 橫跨郯廬斷裂帶， 采用GPS-Corss進(jìn)行定位， 水平定位誤差控制在5 cm以內(nèi)． 儀器采用GEIWSR-II型無纜自定位地震儀和CDJ-S2C型2 Hz三分量地震檢波器， 采樣率500 Hz， 每個(gè)觀測(cè)點(diǎn)連續(xù)記錄20分鐘．

圖2 微震動(dòng)觀測(cè)點(diǎn)、 斷裂構(gòu)造及地震勘探測(cè)線分布圖

2.2 數(shù)據(jù)處理

對(duì)于時(shí)長(zhǎng)20分鐘的微振動(dòng)觀測(cè)數(shù)據(jù)， 剔除被場(chǎng)地噪聲明顯干擾的數(shù)據(jù)段或者選擇信號(hào)質(zhì)量較高的時(shí)段進(jìn)行MSPAC和f-k計(jì)算， 并提取頻散曲線． 圖3a為典型微振動(dòng)原始記錄． 可以看出， 各觀測(cè)臺(tái)記錄信號(hào)形態(tài)相似， 具有較強(qiáng)的一致性， 信號(hào)能量在0.1—10 Hz之間分布， 但主要集中在1.4—2.2 Hz （圖3b）． 圖4給出了分別采用f-k法和MSPAC法計(jì)算的頻散曲線. 可以看出，f-k法和MSPAC法提取的頻散曲線趨勢(shì)、 形態(tài)相似， 但細(xì)節(jié)略有差異. 這種差異與頻散曲線提取的算法有關(guān)， 也與人為因素有關(guān); 頻散能量基本上呈帶狀分布. 本文首先確定頻散能量團(tuán)的上、 下限， 再取上、 下限值的均值作為該方法所獲取的頻散曲線． 為保證計(jì)算結(jié)果的穩(wěn)定性， 本文采用了f-k法和MSPAC法提取的頻散曲線均值作為最終觀測(cè)頻散曲線． 在獲取各觀測(cè)點(diǎn)頻散曲線后， 利用式（3）將Vr-d曲線轉(zhuǎn)換成Vx-d曲線， 再對(duì)各測(cè)點(diǎn)的Vx-d曲線進(jìn)行橫向插值、 平滑計(jì)算， 最終獲得二維視速度剖面（徐佩芬等, 2009， 2013） ．

圖3 微振動(dòng)原始記錄（a）及其頻譜（b）

圖4 采用MSPAC法（a）和f-k法（b）計(jì)算得到的頻散能量分布圖； （c） 所提取的觀測(cè)頻散曲線

3 微振動(dòng)勘探結(jié)果

圖5為A3553測(cè)線（圖2）微振動(dòng)勘探視S波速度剖面， 全長(zhǎng)3.173 km． 根據(jù)視S波速度差異性， 在勘探深度范圍內(nèi)， 可識(shí)別出4套地層， 即第四系、 新近系、 白堊系和太古界-元古界地層．A3553剖面的介質(zhì)結(jié)構(gòu)除在縱向上表現(xiàn)出層狀結(jié)構(gòu)外， 在橫向上呈現(xiàn)出強(qiáng)烈的非均勻性． 第四系（黏土為主）底界面除局部有起伏變化外, 整體較為平坦. 新近系（中粗砂為主）底界面起伏明顯， 該地層在A37、A38、A39、A46、A49、A50號(hào)測(cè)點(diǎn)下方界面起伏明顯， 新近系底界面總體由西向東逐漸變淺；A46號(hào)測(cè)點(diǎn)下方存在連續(xù)的視S波速度梯度帶， 推測(cè)該處存在斷點(diǎn)f1-1． 以A46號(hào)測(cè)點(diǎn)為界， 其西側(cè)下伏基巖為上白堊統(tǒng)王氏組砂巖（K2w）， 東側(cè)為太古界—下元古界片麻巖．

圖5 A3553微振動(dòng)勘探視S波速度Vx剖面虛線為視速度等值線， 實(shí)線為主要地層的分界面

圖6為與微振動(dòng)勘探A3553剖面相鄰的人工地震勘探剖面（QL35）（許漢剛等, 2013）, 其位置見圖2所示. 通過微振動(dòng)勘探剖面（圖5）與人工地震勘探剖面（圖6）對(duì)比可以看出， 微振動(dòng)勘探剖面對(duì)主要地層的解譯與人工地震勘探結(jié)果具有較強(qiáng)的一致性， 例如新近系地層頂面、 基巖頂面起伏均由東向西逐漸加深， 且局部變化明顯； 微振動(dòng)勘探剖面A46觀測(cè)點(diǎn)下解譯的斷點(diǎn)f1-1與QL35人工地震勘探剖面上解譯的斷點(diǎn)形態(tài)具有較強(qiáng)的一致性， 斷點(diǎn)位置也較為接近．A3553剖面上斷點(diǎn)f1-1約在2410#樁號(hào)， 而QL35剖面解譯的斷點(diǎn)在2765#樁號(hào), 造成斷點(diǎn)位置不一致的原因， 除兩條勘探測(cè)線位置不重合外， 也與解釋人員對(duì)資料的判讀和識(shí)別有關(guān)． 在A37—A39測(cè)點(diǎn)范圍內(nèi)， 也存在較為明顯的視S波速度變化， 且該處視S波速度異常與QL35人工地震勘探剖面上f5-1斷點(diǎn)（1035#樁號(hào)）位置一致性較好． 但考慮到A37—A39觀測(cè)點(diǎn)分布偏向斷裂一側(cè)， 且未能有效橫跨擬探測(cè)斷裂構(gòu)造（圖2）， 同時(shí)這些觀測(cè)點(diǎn)觀測(cè)時(shí)受到外界干擾（臨S249省道及S49新?lián)P高速）較為明顯， 因此在對(duì)斷裂解譯時(shí)并未對(duì)該處視S波速度異常進(jìn)行解釋．

圖6 QL35測(cè)線人工地震勘探剖面（引自許漢剛等， 2013）

圖7為A0121測(cè)線（圖2）微振動(dòng)勘探視S波速度剖面． 圖中細(xì)虛線表示視S波速度等值線， 剖面全長(zhǎng)2.778 km． 根據(jù)視S波速度的差異性及視速度梯度變化， 在勘探深度范圍內(nèi)， 可識(shí)別出4套地層， 即第四系、 新近系、 上白堊統(tǒng)和下白堊統(tǒng)地層． 第四系（黏土為主）底界面除局部有起伏變化外， 整體較為平坦； 新近系（中粗砂為主）底界面起伏明顯， 埋深總體由西向東逐漸加深， 受斷裂作用明顯； 上白堊統(tǒng)王氏組砂巖（K2w）和下白堊統(tǒng)青山組凝灰?guī)r（K1q）頂界面埋深起伏明顯， 可能與斷裂構(gòu)造發(fā)育有關(guān)． 在A8、A9、A10、A14、A15測(cè)點(diǎn)下面， 均存在顯著的視S波速度梯度帶， 特別是A14和A15測(cè)點(diǎn)下， 有連續(xù)多個(gè)地層（K--Q）均為視S波速度梯度帶， 推測(cè)該處可能存在斷點(diǎn)f2； 在A9測(cè)點(diǎn)下面， 從下白堊統(tǒng)（K2w）到第四系（Q）也為視S波速度梯度帶， 但下白堊統(tǒng)（K1q）地層連續(xù)， 推測(cè)該處存在斷點(diǎn)f2-1． 從視S波速度梯度帶可以看出， 斷點(diǎn)f2-1切割深度及規(guī)模均小于斷點(diǎn)f2．

圖7 A0121微振動(dòng)勘探視S波速度剖面

圖8 測(cè)線人工地震勘探剖面QL15（引自許漢剛等， 2013）

圖8為與微振動(dòng)勘探A0121剖面位置相鄰的人工地震勘探剖面QL15． 通過微振動(dòng)勘探剖面（圖7）與人工地震勘探剖面（圖8）對(duì)比可以看出，A0121剖面對(duì)主要地層及斷裂構(gòu)造的解譯與QL15剖面解譯結(jié)果具有較好的一致性， 例如上白堊統(tǒng)地層頂面埋深由西向東逐漸加深、 新近系地層頂面起伏不明顯的特征等．A0121剖面A9觀測(cè)點(diǎn)下面解譯的斷點(diǎn)f2-1與QL15剖面1784#樁號(hào)下方解譯的斷點(diǎn)f2-1位置一致性也較好；A0121剖面A14和A15測(cè)點(diǎn)下面的解譯斷點(diǎn)f2與QL15剖面2488#樁號(hào)下方基巖頂面反射波組不連續(xù)的位置也具有較強(qiáng)的一致性．

4 討論與結(jié)論

微振動(dòng)勘探A0121剖面、A3553剖面與人工地震勘探QL15剖面、QL35剖面位置相鄰或相近， 兩種不同的勘探方法雖對(duì)主要地層及斷裂的解譯具有較好的一致性， 但對(duì)地層或斷裂細(xì)節(jié)特征的反映則有所差異．

人工地震勘探剖面對(duì)斷裂空間位置具有精確的定位能力， 如對(duì)斷點(diǎn)的定位可以精確到樁號(hào)， 空間誤差控制在幾個(gè)共深度點(diǎn)（CDP）范圍內(nèi)， 也能較好地揭示各地層之間的接觸關(guān)系和斷裂的發(fā)育形態(tài)等， 如QL35剖面在斷裂F1發(fā)育處顯示上白堊統(tǒng)王氏組砂巖與太古界─下元古界片麻巖呈不整合接觸關(guān)系， 并表現(xiàn)出擠壓逆沖的特征. 同時(shí)人工地震勘探對(duì)地層埋藏深度的確定也有較高的精度， 可以獲取各層位的絕對(duì)埋深． 受激發(fā)條件所限， 人工地震勘探對(duì)基巖內(nèi)部層位的確定能力明顯較弱， 因而對(duì)斷裂構(gòu)造較深層位發(fā)育特征的確定受到一定的限制．

微振動(dòng)勘探方法最顯著的特征是勘探深度大， 可以有效揭示基巖內(nèi)部層位的發(fā)育． 例如A0121剖面、A3553剖面的有效勘探深度可達(dá)900 m左右， 且給出勘探區(qū)域4套地層的發(fā)育特征、 接觸關(guān)系、 埋深及相對(duì)厚度等信息．A0121剖面顯示下白堊統(tǒng)青山組（K1q） 地層頂面埋深約600 m左右， 上白堊統(tǒng)王氏組（K2w）地層厚度約400 m左右；A3553剖面顯示上白堊統(tǒng)王氏組（K2w）地層厚度約400 m左右， 其頂面埋深由西向東逐漸變淺， 并在F1斷裂附近呈尖滅形態(tài)， 超覆于太古界─下元古界片麻巖之上． 上述結(jié)果與前人在該區(qū)所獲取的勘探結(jié)果相一致（謝瑞征等， 1994）．

得益于較深的勘探深度， 微振動(dòng)勘探剖面可以較清楚地揭示斷裂構(gòu)造在淺部地層的發(fā)育特征及其形態(tài)． 例如斷點(diǎn)f2-1在微振動(dòng)勘探剖面和人工地震勘探剖面均有顯示， 但從人工地震勘探剖面上并不能識(shí)別出斷點(diǎn)f2-1下方下白堊統(tǒng)青山組（K1q）地層的連續(xù)性， 只顯示出上白堊統(tǒng)王氏組（K2w）地層頂面的不連續(xù)性. 由A0121剖面可以看出，f2斷點(diǎn)不僅錯(cuò)斷了第四系、 新近系及上白堊統(tǒng)地層， 也錯(cuò)斷了下白堊統(tǒng)地層． 就斷點(diǎn)規(guī)模而言，f2斷點(diǎn)的規(guī)模明顯大f2-1斷點(diǎn)， 推測(cè)f2斷點(diǎn)為主斷面的可能性較大．A3553剖面和對(duì)斷點(diǎn)f1-1的識(shí)別與QL35剖面解譯結(jié)果具有較強(qiáng)的一致性. 其原因之一就是斷點(diǎn)f1-1兩側(cè)巖性存在顯著差異， 微振動(dòng)勘探方法對(duì)斷裂構(gòu)造解譯主要依據(jù)視S波速度梯度， 因而對(duì)巖性差異面的識(shí)別具有較好的分辨率． 微振動(dòng)勘探方法的空間分辨率與觀測(cè)臺(tái)陣半徑有關(guān)， 而觀測(cè)半徑又與勘探深度相關(guān)， 因此在確定了勘探目的層和微振動(dòng)勘探觀測(cè)系統(tǒng)后， 微振動(dòng)勘探的空間分辨率就基本確定了． 就空間分辨率而言， 淺層人工地震勘探法適合對(duì)斷裂構(gòu)造的詳查， 而微振動(dòng)勘探方法較適合對(duì)斷裂構(gòu)造普查．

綜上所述， 微振動(dòng)勘探方法對(duì)城市區(qū)隱伏斷裂具有較好的空間定位能力， 特別是對(duì)斷裂作用下的巖性差異面定位與人工地震勘探定位精度相當(dāng)． 得益于其較深的勘探深度， 微振動(dòng)勘探可以有效揭示勘探區(qū)域下覆地層及斷裂的空間展布特征， 對(duì)識(shí)別和判斷主斷裂面有重要意義， 但其對(duì)地層或斷裂細(xì)節(jié)特征的揭示能力較差． 因此， 在實(shí)際地層或斷裂勘探中， 可選用人工地震勘探方法與微振動(dòng)勘探方法相結(jié)合的方法， 二者取長(zhǎng)補(bǔ)短， 以獲取斷裂構(gòu)造的深淺層發(fā)育特征．

微振動(dòng)成像方法是以視S波速度來反映地層或斷裂特征， 視S波速度可以有效突出介質(zhì)結(jié)構(gòu)的差異性以及減少反演過程中的人為因素干擾， 但在對(duì)視S波速度剖面進(jìn)行量化解釋時(shí)， 則需要參考其它資料．

審稿專家為本文的完善提出了合理化建議， 浙江省地震工程研究所和江蘇省地震工程研究院的同仁參與微振動(dòng)野外觀測(cè)， 吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院提供了觀測(cè)儀器， 作者在此一并表示衷心的感謝．

董連成, 陶夏新, 師黎靜, 李廣影. 2008. 利用地脈動(dòng)臺(tái)陣記錄反演場(chǎng)地淺層VS結(jié)構(gòu)[J]. 巖土力學(xué), 29（2）: 553--556.

Dong L C, Tao X X, Shi L J, Li G Y. 2008. Inversion of shallow shear wave velocity profile of engineering site from recordings associated with microtremors array[J].RockandSoilMechanics, 29（2）: 553--556 （in Chinese）.

何正勤, 丁志峰, 賈輝, 葉太蘭. 2007. 用微動(dòng)中的面波信息探測(cè)地殼淺部的速度結(jié)構(gòu)[J]. 地球物理學(xué)報(bào), 50（2）: 492--498.

He Z Q, Ding Z F, Jia H, Ye T L. 2007. To determine the velocity structure of shallow crust with surface wave information in microtremors[J].ChineseJournalofGeophysics, 50（2）: 492--498 （in Chinese）.

劉保金, 張先康, 方盛明, 趙成斌, 段永紅, 朱金芳, 黃昭, 黃宗林, 王善雄, 鄭德剛. 2002. 城市活斷層探測(cè)的高分辨率淺層地震數(shù)據(jù)采集技術(shù)[J]. 地震地質(zhì), 24（4）: 524--532.

Liu B J, Zhang X K, Fang S M, Zhao C B, Duan Y H, Zhu J F, Huang Z, Huang Z L, Wang S X, Zheng D G. 2002. Acquisition technique of high-resolution shallow seismic data for surveying of urban active faults[J].SeismologyandGeology, 24（4）: 524--532 （in Chinese）.

潘紀(jì)順, 劉保金, 朱金芳, 張先康, 方盛明, 王夫運(yùn), 段永紅, 徐朝繁. 2002. 城市活斷層高分辨率地震勘探震源對(duì)比試驗(yàn)研究[J]. 地震地質(zhì), 24（4）: 533--541.

Pan J S, Liu B J, Zhu J F, Zhang X K, Fang S M, Wang F Y, Duan Y H, Xu Z F. 2002. Comparative experiment on seismic sources in high-resolution seismic exploration for urban active faults[J].SeismologyandGeology, 24（4）: 533--541 （in Chinese）.

師黎靜, 陶夏新, 趙紀(jì)生. 2006. 地脈動(dòng)臺(tái)陣方法的有效性分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 25（8）: 1683--1690.

Shi L J, Tao X X, Zhao J S. 2006. Validation of microtremors array method[J].ChineseJournalofRockMechanicsandEngineering, 25（8）: 1683--1690 （in Chinese）.

陶夏新, 劉曾武, 郭明珠, 師黎靜, 董連成. 2001. 工程場(chǎng)地條件評(píng)定中的地脈動(dòng)研究[J]. 地震工程與工程振動(dòng), 21（4）: 18--23.

Tao X X, Liu Z W, Guo M Z, Shi L J, Dong L C. 2001. A review of microtremor study in engineering site rating[J].EarthquakeEngineeringandEngineeringVibration, 21（4）: 18--23 （in Chinese）.

王偉君, 劉瀾波, 陳棋福, 張杰. 2009. 應(yīng)用微動(dòng)H/V譜比法和臺(tái)陣技術(shù)探測(cè)場(chǎng)地響應(yīng)和淺層速度結(jié)構(gòu)[J]. 地球物理學(xué)報(bào), 52（6）: 1515--1525 （in Chinese）.

Wang W J, Liu L B, Chen Q F, Zhang J. 2009. Applications of microtremorH/Vspectral ratio and array techniques in assessing the site effect and near surface velocity structure[J].ChineseJournalofGeophysics, 52（6）: 1515--1525 （in Chinese）.

王振東. 1986. 微動(dòng)的空間自相關(guān)法及其實(shí)用技術(shù)[J]. 物探與化探, 10（2）: 123--133.

Wang Z D. 1986. The micromotional spatial autocorrelation method and its practical technique[J].Geophysical&GeochemicalExploration, 10（2）: 123--133 （in Chinese）.

謝瑞征, 唐錦鐵, 顧瑾平, 朱永正, 季幼庭, 朱書俊, 黃偉生, 梁中華. 1994. 宿遷市地震動(dòng)小區(qū)劃工作報(bào)告[R]. 南京: 江蘇省地震局: 78.

Xie R Z， Tang J T, Gu J P, Zhu Y Z, Ji Y T, Zhu S J, Huang W S, Liang Z H. 1994.ReportontheWorkofSeismicMicrozonationofSuqian[R]. Nanjing: Earthquake Administration of Jiangsu Porvince: 78 （in Chinese）.

許漢剛, 顧勤平, 趙啟光. 2012. 宿遷市隱伏斷層的控制性人工地震探測(cè)野外工作報(bào)告（JSE-2012-SQ-Z01）[R]. 南京： 江蘇省地震工程研究院: 185.

Xu H G, Gu Q P, Zhao Q G. 2012.ReportontheShallowSeismicExplorationforDetectingtheBuriedFaultofSuqianCity（JSE-2012-SQ-Z01）[R]. Nanjing: Earthquake Engineering Institute of Jiangsu Province: 185 （in Chinese）.

許建聰, 孫紅月. 2004. 地脈動(dòng)作用下的巖土動(dòng)力響應(yīng)研究[J]. 振動(dòng)工程學(xué)報(bào), 17（2）: 147--150.

Xu J C, Sun H Y. 2004. Study on dynamical response of ground soil to micro-tremor[J].JournalofVibrationEngineering, 17（2）: 147--150 （in Chinese）.

徐佩芬, 李傳金, 凌甦群, 張胤彬, 侯超, 孫勇軍. 2009. 利用微動(dòng)勘察方法探測(cè)煤礦陷落柱[J]. 地球物理學(xué)報(bào), 52（7）: 1923--1930.

Xu P F, Li C J, Ling S Q, Zhang Y B, Hou C, Sun Y J. 2009. Mapping collapsed columns in coal mines utilizing microtremor survey methods[J].ChineseJournalofGeophysics, 52（7）: 1923--1930 （in Chinese）.

徐佩芬, 侍文, 凌蘇群, 郭慧麗, 李志華. 2012. 二維微動(dòng)剖面探測(cè)“孤石”: 以深圳地鐵7號(hào)線為例[J]. 地球物理學(xué)報(bào), 55（6）: 2120--2128.

Xu P F, Shi W, Ling S Q, Guo H L, Li Z H. 2012. Mapping spherically weathered “Boulders” using 2D microtremor profiling method: A case study along subway line 7 in Shenzhen[J].ChineseJournalofGeophysics, 55（6）: 2120--2128 （in Chinese）.

徐佩芬, 李世豪, 杜建國(guó), 凌蘇群, 郭慧麗, 田寶卿. 2013. 微動(dòng)探測(cè): 地層分層和隱伏斷裂構(gòu)造探測(cè)的新方法[J]. 巖石學(xué)報(bào), 29（5）: 1841--1845.

Xu P F, Li S H, Du J G, Ling S Q, Guo H L, Tian B Q. 2013. Microtremor survey method: A new geophysical method for dividing strata and detecting the buried fault structures[J].ActaPetrologicaSinica, 29（5）: 1841--1845 （in Chinese）.

于凱, 車愛蘭, 馮少孔, 衫山長(zhǎng)志. 2011. 長(zhǎng)周期地脈動(dòng)觀測(cè)在深部地層調(diào)查中的應(yīng)用[J]. 上海交通大學(xué)學(xué)報(bào), 45（5）: 701--705.

Yu K, Che A L, Feng S K, Siguyama T. 2011. Application in deep stratum of soil investigation using long-period microtremor[J].JournalofShanghaiJiaotongUniversity, 45（5）: 701--705 （in Chinese）.

Aki K. 1957. Space and time spectra of stationary stochastic waves with special reference to microtremors[J].BullEarthqResInstTokyoUniv, 35: 415--456.

Bettig B, Bard P Y, Scherbaum F, Riepl J, Cotton F, Cornou C, Hatzfeld D. 2001. Analysis of dense array noise measurements using the modified spatial auto-correlation method （SPAC）: Application to the Grenoble area[J].BGeofisTeorAppl, 42（3/4）: 281--304.

Capon J. 1969. High-resolution frequency-wavenumber spectrum analysis[J].ProcIEEE, 57（8）: 1408--1418.

Cho I, Tada T, Shinozaki Y. 2006. Centerless circular array method: Inferring phase velocities of Rayleigh waves in broad wavelength ranges using microtremor records[J].JGeophysRes, 111: B09315. doi:10.1029/2005JB004235.

Ling S, Okada H. 1993. An extended use of the spatial autocorrelation method for the estimation of structure using microtremors[C]∥Procofthe89thSEGJCoference. Nagoya: Society of Exploration Geophysicists of Japan: 44--48 （in Japanese）.

Okada H, Suto K. 2003.TheMicrotremorSurveyMethod:GeophysicalMonographSeries[M]. Oklahoma: Society of Exploration Geophysicists: 150.

Okada H. 2006. Theory of efficient array observations of microtremors with special reference to the SPAC method[J].ExplorGeophys, 37（1）: 73--85.

Toks?z M N, Lacoss R T. 1968. Microseisms: Mode structure and sources[J].Science, 159（3817）: 872--873.

Detecting buried fault structures by microtremor survey method

1）EarthquakeAdministrationofJiangsuProvince,Nanjing210014,China2）CollegeofInstrumentation&ElectricalEngineering,JilinUniversity,Changchun130061,China

The methods and schemes of detecting buried fault structures using microtremor survey method are introduced, by which two branch faults of the Tanlu fault zone are explored. Comparison with the shallow seismic exploration results indicates that the spatial locating ability of microtremor survey method is better to detect the buried fault structure, the location error of fault caused by lithological differences is close to that from the results of the shallow seismic exploration. For its deep exploration ability, the microtremor survey method can reveal the spatial distribution characteristics as well as contact relationship between the strata and the faults beneath the microtremor survey arrays, which is important to identify and confirm the main rupture surface of faults. However, the microtremor survey method is difficult to reveal or describe the detail features of faults. To explore the strata or faults, it is the best to combine the shallow seismic exploration method with the microtremor survey method to study the relationship between deep and shallow faults. The microtremor survey method can be used to detect fault roughly, and then the shallow seismic method can be used to determine accurate position of faults after knowing the general location and spatial characteristic of faults.

mircrotremor survey method; buried fault structure; shallow seismic exploration

10.11939/jass.2015.01.012.

江蘇省社會(huì)發(fā)展項(xiàng)目（BE2010753）、 蘇州市“活斷層探測(cè)與地震危險(xiǎn)性評(píng)價(jià)”項(xiàng)目（JSE-2010-SZ）和宿遷市“活斷層探測(cè)與地震危險(xiǎn)性評(píng)價(jià)”項(xiàng)目（JSE-2012-SQ）共同資助.

2014-04-22收到初稿， 2014-08-04決定采用修改稿.

e-mail: nj_fxp@163.com

10.11939/jass.2015.01.012

P315.2

A

范小平, 李細(xì)兵, 趙啟光, 宋浩, 戴波, 張帥帥. 2015. 用微振動(dòng)勘探方法探測(cè)隱伏斷裂構(gòu)造. 地震學(xué)報(bào), 37（1）: 134--143.

Fan X P, Li X B, Zhao Q G, Song H, Dai B, Zhang S S. 2015. Detecting buried fault structures by microtremor survey method.ActaSeismologicaSinica, 37（1）: 134--143. doi:10.11939/jass.2015.01.012.