地震縱波橫波聯合勘探方法尋找近地表第四系內隱伏斷裂方法研究

李顏貴 ， 蔣正中 ， 劉子龍 ， 雍 凡 ， 羅水余

(中國地質科學院 地球物理地球化學勘查研究所，廊坊 065000)

0 前言

隨著人類社會和經濟的飛躍發展，工程建設的規模越來越大，各種自然和人為因素對建筑和工程設施的影響作用也越來越強，其中以活動斷裂為主的地質災害的影響日益成為各國的高度重視。在區域地殼穩定性研究中，活動斷裂成為重要的考慮因素，活動斷裂的活動性主要反映在對地面裂縫、地面沉降以及地面塌陷的控制和誘導作用并可能引起地震的發生，其形成的條件是活動斷裂由于地應力長期的集中，在一定條件下的突然釋放。或者活動斷裂長期緩慢的蠕滑而造成地表產生地裂縫、地面塌陷等。京津地區有大量的活動斷裂和大的構造，所以研究京津地區地殼的穩定性，防范地質災害，就要充分研究京津地區的活動構造，尤其是隱伏在平原區松散層之下的活動斷裂，掌握其分布范圍、活動規律以及影響范圍和影響程度。在活動斷裂的調查方面，前人做了大量工作，徐明才等[1-5]進行了大量縱波尋找活動斷裂方面的研究，取得了顯著效果；在利用橫波尋找活動斷裂方面，白云等[6]中國地震局蘭州地震研究所[7]等進行了初步嘗試，取得了一定的效果；RogerA. Young[8]利用SH波在美國Oklahoma州Norman附近的一個老垃圾場進行橫波調查，在使用錘擊震源取得了較好的勘探效果；使用重錘激發，能量有限，激發能量人為因素影響大，勘探的深度和效果不理想，同時激發的主頻低，阻礙橫波應用[9]，但這局限性為研究橫波勘探活動斷裂提供了很好的研究方向：使用能量大的可控震源，提高激發的能量和頻率；因此本次實驗選擇能量更大的可控震源激發橫波，通過調整震源震動方向及橫波檢波器放置方向來接收SH、SV兩種橫波進行研究，通過使用可控震源，縱橫波勘探深度和資料的信噪比有了明顯提高，對準確尋找地表第四紀未固結沉積層中的隱伏斷裂上斷裂點提供了技術支撐。

1 地質概況與地震地質條件

縱、橫波實驗測線位于北京馬坊環島，與2011年進行的“北京平原區活動斷裂監測專項地質調查項目”中地震勘探項目WT26線位于同一位置，該斷裂為夏墊斷裂的北部，基巖埋深較淺，夏墊斷裂是北京平原區一條重要的活動斷裂，北起平谷馬坊、夏墊、西集、鳳河營到廊坊以西一帶。走向NE30°，傾向SE，傾角為65°～70°，為正斷層。1679年平谷馬坊8級大地震就發生在該斷裂帶上，是一條活動斷裂，在第四系地層內，有不同顆粒的砂層、粘土層及亞粘土層等不同巖性的地層。由此可推測，除該區基巖界面產生的反射信號外，蓋層內反射地震信號將比較豐富。除了粘土和細砂層之間界面的波阻抗差異能引起反射外，砂層中夾的粘土薄層和粘土層中夾的薄砂層都能產生地震響應而被記錄下來，風化層與未風化的基巖地層之間的速度和密度差異較大，這些界面也能夠形成較強的反射波被記錄下來。在該地區分布有四口鉆井，其中測線經過ZK8和ZK9鉆井，已有鉆井資料可以對試驗結果進行驗證，粘土層、礫石層和細沙層分布不均，這恰好為地震解釋研究層位的變化提供依據。

2 方法技術

2.1 縱向分辨率與橫向分辨率

在淺層地震地勘探中，地震記錄所反映的各種地質構造的清晰度取決于地震資料的分辨率。地震勘探的分辨率就是分辨各種地質體和地層細節的能力，它包括縱向分辨率和橫向分辨率；縱向分辨率也叫垂向分辨率或時間分辨率，它是指地震記錄沿垂直方向能夠分辨的最薄地層的厚度，通常的含義是從地震記錄上能夠正確地識別地層頂、底界面的反射波，通過經典的薄層模型實驗得到規律。縱向分辨率主要與地震波的波長直接相關，即與地震波的速度和頻率有關，橫波的速度低，在該實驗區內橫波的速度是縱波的0.2倍～0.4倍，通過分析實測頻譜曲線可知，SH波頻率相對較低，SV波頻率和縱波頻率接近，這與詹正彬[10]所敘述的SH、SV頻率一致；SV橫波速度低，主頻和縱波頻率接近，這樣可以提高縱向分辨率，對于區分薄層地層，特別是未固結的地層有較好的分辨能力。橫向分辨率也叫水平分辨率或空間分辨率，它是指地震記錄沿水平方向能夠分辨的最小地質體的寬度。地震波是以球面波的形式向前傳播，當遇到波阻抗界面就產生反射，地震波有個延續時間，根據波的疊加原理，當初至波波前面的時差為1/4個周期內將產生相長干涉，其雙程旅行時間比垂直入射到界面的反射被的雙程旅行時間晚半個周期以上，因而產生相長干涉。我們把產生反射波相長干涉的有效界面稱為第一菲涅爾帶，理論證明可知其半徑：

式中：r第一菲涅爾帶半徑；v為地震傳播速度；h為界面深度；f為地震波主頻。由于橫波的速度較低，特別是在為固結的軟土地層，縱波的速度比橫波的速度高幾倍，在主頻變化不大的情況下，橫波的橫向分辨率要比縱波高，從圖3速度譜就可以清楚看到。

2.2 數據采集技術

合理的數據采集技術是提高淺層高分辨率地震勘探的基礎，橫波的激發常采用扣板法激發，能量較小，且適用條件較差；本次實驗采用美國IVI公司制造的MiniVib系列T2500先進可控震源，包含了從縱波輸出到橫波輸出的一套組件，既能激發常規縱波，又能激發橫波，獨特的“minivib”震源橫波輸出技術可以讓操作者隨意調節震源器與接收排列之間的角度，寬頻帶5 Hz～550 Hz的高頻輸出，體積小、重量輕，是理想的淺層地震勘探震源。縱波采用96道、零偏移距、3 m道距，具體參數見表1。橫波通過調整震源震動方向及橫波檢波器放置方向來激發和接收SH、SV兩種橫波，激發SH波和P波野外采集一樣，只需要調整震源的震動方向與檢波器的方向一致，即90°垂直測線方向，激發方向沿Y方向作用的水平力源只能激發SH波，布設橫波檢波器水平芯體正方向指向Y正方向，檢波器在炮點兩端布置一樣，采集的單炮記錄存在震源的兩邊有極性相反的情況，道編輯改變極性。每種橫波單獨激發采集，這樣避免了調整角度的問題；在利用可控震源進行數據采集前，進行野外數據采集試驗：掃描類型試驗、掃描頻率試驗和掃描時間試驗等，采集參數見表1；另外，還進行了垂直疊加次數的試驗：為壓制干擾，除利用可控震源的掃描疊加壓制干擾外，還利用垂直疊加技術(每個物理點垂直疊加2 次～3次)進一步壓制噪聲，提高地震記錄的信噪比，通過增加震動次數和掃描時長，保證在同一地點多次激發能量均衡穩定、振幅譜和自相關穩定；由于縱波勘探深度深，線性升頻掃描的信號在向下傳播的過程中，縱波的高頻部分損失較大，降低了淺層的分辨率；通過分析橫波單炮記錄的振幅、能量可知，在中淺部，地層對能量的衰減是呈現指數型衰減，設計采集衰減系數倒數的指數型掃描信號，彌補了高頻信號能量，壓制了掃描信號的自相關旁瓣幅度；采用地面力回饋信號相關，保證接收到的信號是地層真實的反映；由于激發橫波震源的震動方向在水平方向，震源的能量輸出要相對小些，通過實驗可知，當增加能量輸出在60% 時，8邊形鋸齒與地面耦合發生位移，單炮記錄存在震源干擾，調整掃頻方式為線性掃描和震源輸出能量為40%，震源干擾減弱或消失；另外，由于橫波激發和接收的方向特性匹配原則[10]，激發SH波沿Y方向作用的水平力源，這樣SH波的能量最強，SV波和P波最弱為零；同理，x方向作用的水平力主要產生SV波，遠偏移距，入射角度變大，會產生能量強的轉換PV波，處理較困難；通過實驗及計算，60道180 m偏移距采集，能減弱轉換波在大偏移距的影響；本次實驗為揭露的斷層深度在100 m以內，為了驗證斷層從深至淺的連續性，通過擴展排列實驗和干擾波調查，采用“最佳時窗”原則，使用96道接收，非對稱中間激發， 3 m小道間距，12 m炮距，水平覆蓋次數12次。

表1 縱橫波采集參數表

2.3 數據處理

地震數據處理主要包括：解編、定義觀測系統、編輯、振幅補償、疊前去噪、速度分析、動校正、疊加、疊后去噪、偏移等。為達到理想的處理效果，在正式處理數據前，進行處理參數分析和處理效果對比試驗，選用有效的處理模塊編排流程，選擇了最佳的處理參數。本次資料處理工作也借鑒了以往深層地震探測資料處理解釋的經驗，并做了一些改進；通過試處理，發現有效的疊前去噪技術、速度分析技術是此次數據處理流程中的關鍵步驟，整個縱波數據處理流程如圖1所示。橫波的數據處理有其本身的特點，傳播的方式和速度不同，但數據處理的基本方法是相同的，處理流程如圖1所示，去噪、振幅恢復采用些非常規的方法處理，處理的難點主要：①橫波的速度分析：橫波速度低，速度在150 m/s～300 m/s，本次處理軟件采用的是帕拉代姆帕拉代姆公司Focus5.4處理軟件，拾取速度和速度掃描最低速度在320 m/s，不能使用該模塊；②橫波的靜校正：雖然地表起伏不大，但是橫波的靜校正求取與常規的不同，特別是橫波的折射波，其上面的P波、初至波等干擾波比橫波折射波速度高，如圖3(b)速度分析上面速度較高的就是P波反射波和折射波，LOVE波干擾大偏移距的SH反射波信號及折射波，折射波識別和提取困難，靜校正較難，如圖3(b)所示單炮記錄200 ms以上各種干擾波(反射波、折射波、震源線性干擾波、面波和初至波)；③橫波振幅恢復和能量補償：在潛水面處有很強的波阻抗界面，其上伏地層因其覆蓋次數低疊加效果不好，在其下屏蔽了較深層的反射波，能量很弱，看不到反射層，如圖2(a)SH單炮記錄所示；④橫波數據弱信號的提取，SH的反射波在大偏移距信號受Love波干擾，由于其速度接近與SH波，常規的FK域去噪差，而在小偏移距內沒有受到Love波干擾，但有線性干擾，使反射波難以識別；⑤反褶積，對可控震源激發的地震信號進行反褶積有兩種途徑，直接對震動記錄進行，另外是對相關記錄進行，現在的處理基本上是最小相位的反褶積，而反射波的基本上是混合相位的，存在淺層應用反褶積效果反而變差的情況。

根據以上難點，采用了適合本工區的處理方法和處理流程：①對速度分析和常速掃描，對現有觀測系統進行改進，通過理論計算，提升速度10倍，也滿足地震勘探原理，這樣處理模塊就能正常使用，經過對比，效果較好；②采用先去強噪聲后，在進行地表一致性振幅恢復，利用速度抬升后，拾取的速度進行速度和頻率域二維濾波、相關濾波、傾角FK濾波和疊前多域去噪(干擾波在不同的域上的表現形式不同，分別在炮域、共檢波點域、共偏移距域去噪)，如圖3(b)所示，去噪效果明顯，速度譜能量更集中，拾取更準確，反射波更清晰，特別是深層的反射信號；

圖1 縱橫波數據處理流程圖Fig．1 P and S wave processing flow chart

圖2 SH波多域去噪、壓制面波、線性干擾和震源干擾對比剖面Fig．2 SH pave gather is compared with multi-domain suppress interface wave, linearity noise and source noise after suppressing noise(a) SH單炮記錄；(b)CDP域去噪；(c)炮域串聯反褶積；(d)多域組合LIFT去噪

圖3 速度抬升前后壓噪及速度分析對比剖面Fig．3 Gather is compared with suppressing noise and velocity analysis after velocity uplifts(a) 未速度抬升；(b)速度抬升后；(c) NMO

③橫波的頻率和速度與Love波頻率和速度接近，采用頻率域和速度域去除效果不好，但在小偏移距內，SH反射波是標準的雙曲線型，通過拉東變換和速度拉伸切除的方法，去除大偏移距內的噪聲，在小偏移距內，利用線性去噪、傾角FK濾波去除震源線性干擾，對比圖2(d)和3(b)，中間的線性干擾能量很弱，小偏移距內的Love波去除的較干凈，如圖3(c)所示，600 ms大偏移距內的Love波經NMO后，動校拉伸切除和疊加后可以去除大部分Love波；④采用預測反褶積和零相位反褶積組合串聯的方法，可以看出效果明顯，如圖2(c)所示，可以看出，采用地表一致性振幅恢復前后對比圖，300 ms～600 ms的反射波較清晰；疊后的數據采用了譜白化，從而抬升了去噪和疊加后對信號的主頻降低，從而提高了信號的分辨能力，處理過程中沒有加入混波參數，混波雖然能增強振幅，但容易丟失有用的地震振幅信息；⑤橫波在震源兩側，存在極性反轉，采用道編輯極性反轉，或采用負偏移距整體靜校正，校正量為正偏移距的半個周期，消除極性相反的影響。

3 應用效果與分析

如圖6所示，根據該剖面的波組特征，解釋了3組反射波T0、T1、T0-1。結合鉆孔柱狀對比圖(圖5)，可以推斷出T0為第四系底界面產生的反射波,反射波的深度在180 m，T1為白云巖頂界面反射波，深度在240 m左右；T0-1第四系內的粘土與粗沙層產生的強波阻抗界面產生的反射波深度在80 m左右。這些反射波組連續性較好,可連續追蹤對比。CDP108、CDP195和 CDP280處解釋有3條斷層F1-F3,其中F1和F2條斷層都使T1反射波發生了錯斷，而F3沒有，但上面的反射波組有相應的不連續，不能準確地判斷是否延伸到了第四系內部，F2使T0-1反射波發生了錯斷，但是具體錯斷到哪個深度和層位，縱波難以確定，這需要勘探精度更高的橫波勘探進行解釋，這3條斷層相距較近，推測為同一組斷層。F2斷層的可靠性較高，與夏墊斷裂的中深部構造對應較好，夏墊斷裂為正斷裂，從走向、傾向SE和傾角，結合縱波的勘探結果確定的布置鉆井的位置；在斷層附近，反射波的信噪比較低，連續性較差，結合圖4和ZK8～ZK9可知，地層在244 m處白云巖層見棕黃色含礫粘土，222 m處見次級斷裂，且見斷層擦痕，第四系內部，地層變化較劇烈，對于第四系內部的分層縱波分層已難以解決，T0-1反射層部分不連續可能是地層變化劇烈造成的，如果第四系內部小的錯段，不會引起斷層附近的地層發生明顯破碎，大體上會填充斷層泥，而且需要結合橫波的資料。

圖4 第四系等厚線圖Fig．4 Quaternary system isopach map

圖5 鉆孔地層剖面柱狀對比圖Fig．5 Drilling stratigraphic correlation diagram

如圖7和圖8所示，橫波的勘探深度在150 m以內，反射波受斷層破碎的影響較小，其信噪比相對較高。根據地震剖面上的波組特征，結合已掌握的地質資料，認為F1～F3斷層為本次實驗重點探測的夏墊斷裂。根據測區第四系等厚線圖(圖4)和鉆孔柱狀對比圖(圖5)，在圖5的粘土層和礫石層接觸面，有強反射，且和深度剖面對應較好，SH波和SV波的疊加剖面的分辨率有較大差異，是不同類型的波反映的介質特征的差異所致，推斷SV波在斷層附近發生橢圓偏振[10]，斷層附近反生轉換。在已知鉆井驗證、這幾種波的地震特征對比后，聯合探測解釋，主要是相位和振幅，特別在斷層F2和斷層F3處，結合鉆孔揭露層位對比，地層的標定較強的波阻抗差異的層位，橫波在第四系內部分三個層位，從淺到深T0-3、T0-2和T0-1，斷層F2在SH剖面上顯示到達了40 m，SV剖面是60 m，不同類型的橫波反映的介質特性有差別，SV波在F2斷層在60 m處振幅和速度變化不明顯，但SH波明顯，地下地層結構復雜，具體產生的原因，需要探討，在F3斷層附近兩種橫波都較明顯，速度和振幅都發生變化，深度和位置對應交好(圖3)，設礫石層速度等于導出的均方根速度(245 m/s)，SV波主頻為70 Hz，那么粘土層1/4波長厚度就為0.88 m，SH主頻為50 Hz，粘土層1/4波長厚度就為1.23 m，即鉆孔內土層的變化在3 m內都可以在橫波的剖面上反映出來。ZK8和ZK9在第四系內部粘土層分布不均勻，變化較大，橫波剖面的強反射同向軸與粘土層位置、深度變化一致，說明使用橫波勘查第四系土層的變化層位準確、深度可靠。橫波地震剖面推斷F1-F3斷裂，SH波深度剖面F3斷裂已經延伸到了20 m左右層位，因上層反射波層位信噪比較低，已難判斷是否延伸更淺層位；SV波深度剖面在20 m左右層位，其上反射波層位信噪比相對較高，推斷斷層已延伸到該層位，結合SH波資料，綜合二種橫波剖面，推斷F3斷層延伸到了20 m，在T0-3層位之上；結合地質資料，其地層應該是Q3p的晚期，這與地質、電法資料和鉆探揭露的一致。

4 初步結論

通過縱、橫波試驗研究，結合鉆孔層位資料，基本上控制了夏墊斷裂在第四系內部斷層的走向、具體延伸位置及其活動特性，對斷層在第四系內部發育情況有了新的認識：橫波勘探發現了第四系內部發育一定數量的次生隱伏斷裂，并且斷裂延伸到了更新世晚期，斷層深度較淺，橫波勘探對于第四系內部的分層有較好的效果。

圖6 3 m縱波深度地震剖面Fig．6 3 m P wave depth seismic section

圖7 3 m SV波深度地震剖面Fig．7 3 m SV wave depth seismic section

圖8 3 m SH波深度地震剖面Fig．8 3 m SH wave depth seismic section

鑒于橫波資料的特點，橫波的去噪和速度分析是橫波數據處理的難點，通過對橫波資料處理、解釋及與鉆井資料對比結果，尋找近地表第四系內隱伏斷裂方法橫波勘探效果明顯，但還有許多問題值得進一步研究：淺層縱波及橫波深度偏移后基本沒有變化，這與以往的處理有很大差別，推斷可能是速度變化較小。通過開展縱、橫波淺層地震勘探尋找隱伏斷裂方法研究也取得了一些經驗：①從數據采集入手，采用高性能的勘探設備，適宜的震源、最佳的觀測系統，利用合理的抗干擾手段，獲得最佳地震記錄；②淺層地震數據處理以提高分辨率為主要目的，淺層地震數據在達到要求分辨率的基礎上，減少修飾手段，如小的異常經過混波后，異常消失或減弱；在合適的時窗下進行速度、頻率域二維濾波、相關濾波及疊前多域去噪，對提高分辨率有較好的效果；處理過程中對數據要進行時頻分析，去噪以適合為主；速度抬升后的速度分析、常速掃描提高波速拾取準確度，疊加速度準確；③橫波勘探在淺層活斷層調查中具有明顯優勢，縱向分辨率和橫向分辨率都較高，能分辨薄的地層且定位準確。在條件允許的情況下進行多波勘探，能夠彌補傳統縱波勘探丟失的近地表沉積層信息，聯合其他地球物理方法，如高密度電法，間接研究斷裂附近的土層的破碎程度，清晰掌握地下結構由淺至深的空間展布形態有一定效果。

致謝

特別致謝北京地調院王志輝工程師提供的部分地質和鉆井資料。

參考文獻：

[1] 高景華, 徐明才.地震方法確定活動斷裂上斷點的影響因素分析[J].地震地質, 2008, 30(2):338-347．

[2] 高景華, 徐明才.利用地震剖面研究夏墊斷裂西南段的活動性[J].地震地質, 2008, 30(2):497-504．

[3] 徐明才, 高景華.地震方法在斷層活動性評價中的應用研究[J].地震學報, 1998, 20(3):279-284．

[4] 王慶海, 徐明才.地震反射技術在工程地質調查中的應用[J].工程勘察, 1987, 15(2):60-63．

[5] 高景華, 徐明才.高分辨率城市活斷層地震資料處理技術[J].物探化探計算技術, 2005, 27(3):209-216．

[6] 白云, 鐘以章, 盧造勛.地震S波反射法在城市活動斷層探測中的應用與效果研究[J].東北地震研究, 2001, 17(2):69-74．

[7] 盧育霞, 王振明.蘭州盆地活斷層的高分辨率淺層地震勘探技術[J].西北地震學報, 2008, 33 (4):344-353．

[8] ROGER A YOUNG.Near-surface, SH-wave surveys in unconsolidated, alluvial sediments[J].The Leading EDGE, 2001, 20(9):936-948．

[9] 施逸忠，王式.淺層地震橫波反射法最佳應用條件研究[J].河海大學學報, 1995, 23(6):16-21.

[10] 詹正彬，姚姚.多波及橫波地震勘探[M].北京：地質出版社，1994.