耿祥峰(1.山東省地礦局八〇一水文地質工程地質大隊(山東省地礦工程勘察院),山東 濟南 250014;2.山東省地下水環境保護與修復工程技術研究中心,山東 濟南 250014)
國內外探測巖溶塌陷常用的地球物理方法包括:常規電法勘探、地質雷達法、重力測量方法、地震勘探方法等。高密度電法勘探是目前應用效果最好的電阻率法,具有較高的分辨能力。地質雷達具有分辨率高、施工快速方便以及適應能較強等特點,在地質地球物理條件適宜的情況下,成為巖溶塌陷探測的常規手段之一。重力勘探方法最初應用于區域地質構造、金屬礦和非金屬礦勘探,隨著微伽級高精度重力儀的誕生,重力測量精度大幅提高,其分辨能力和對微弱異常的識別能力大大提高,已經成功應用于溶洞、地下空洞調查。瞬變電磁法、電磁波CT、多道瞬態瑞雷波技術以及可控源音頻大地電磁法等,在巖溶塌陷勘察中得到了良好的應用效果。
反射地震勘探技術是依據地下介質在物性差異界面上地震波的動力學以及運動學特征,探測地層或者基巖的埋深及其速度結構,具有分辨率高、探測結果可靠等特點,在巖溶塌陷勘察中具有良好的應用前景。本文選擇淺層地震勘探反射波法對溶塌陷區進行精細探測,取得良好的應用效果。
巖溶發育首先應具有可溶性巖層,包括碳酸鹽類巖石(如白云巖、石灰巖)、硫酸鹽類巖石(如石膏)以及鹵素類巖石(如巖鹽)等,其次還應具備溶解力和足夠流量的水,同時還要有地表水以及地下水的運移通道。巖溶發育區第四系地層、基巖和溶洞、裂隙以及充填黏土、水等介質的溶洞、裂隙之間,其密度和地震波傳播速度具有明顯差異,形成多個波阻抗界面,地震波通過時會形成不同特征的反射波信號。其異常特征如下:
(1)反射波出現半波損失,地震波相位反轉。

式中:ρ1、V1分別為上地層的密度和地震波傳播速度;ρ2、V2分別為下地層的密度和地震波傳播速度;正常地層層序中,下伏巖層的密度和地震波傳播速度均大于上地層,此時反射系數為正值。當地層中存在巖溶或者裂隙發育時,其密度以及地震波傳播速度低于圍巖,反射系數為負值,地震剖面上巖溶塌陷處地震波會出現相位反轉、相位不連續現象。
(2)反射波同相軸錯斷、反射凌亂或者出現空白反射。受重力作用的影響,地層中存在巖溶或者裂隙時,不完整巖層區域下沉,致使反射波同相軸出現錯斷現象,巖溶發育區域地層層序規律遭到破壞,出現雜亂無章的反射波。若存在一定規模的空洞,地震波經過時會出現空白反射區域。
(3)振幅變化特征。地層中巖溶或者裂隙發育程度較低時,巖溶發育區域為黏性土、碎石等充填,地震波經過時,能量衰減迅速,致使反射波能量變弱,表現為反射波振幅變小。地層中存在充填空氣的空洞時,其反射系數可視為-1,地震波經過時,幾乎沒有透射波能量損失,相對于圍巖反射波特征,其反射能量加強,表現為反射波振幅變大。
(4)繞射波。地層錯動造成的棱角點、地層尖滅點、洞穴等,往往都能滿足繞射波的形成條件。繞射波時距曲線的極小值點位于繞射點的正上方,且動校正后繞射波時距曲線仍然是雙曲線。
(5)第四系地層與基巖分界面反射波同相軸特征。受地下水沖刷作用影響,溶塌陷區內第四系地層中黏性土、碎石等會隨地下水運移,第四系地層底界面,反射波往往會出現空白反射或者弱反射現象。
巖溶塌陷區溶洞或者裂隙的規模不大,為了提高分辨率和探測成果的可靠性,地震反射波勘探中一般選擇高頻檢波器,寬頻接收,采用小道間距和小炮間距。為了盡可能減少復雜地形的影響,多采用小排列接收。下面就巖溶塌陷區中地震反射波法勘察實例,討論地震反射波法在巖溶塌陷區的應用效果。
本節給出的是山東省濟南市巖溶塌陷區地震反射波法勘探實例,塌陷區上覆第四系地層,厚度10余米,下伏奧陶系灰巖、泥灰巖,厚度數百米。這種地質條件下開展淺層反射波地震勘探的難點是:第四系地層埋深淺,淺部的有效信號容易被震源附近強烈的干擾波覆蓋,很難區分基巖界面頂板埋深。針對此情形,采用2 m道間距,4 m偏移距的接收方式,16次水平覆蓋觀測系統[1]。
根據野外采集的數據,通過一系列數據處理流程,得到了反射波組清晰的時間剖面圖。巖溶塌陷區地震勘探結果如圖1所示,反射波組中可以劃分4個反射層位,分別標注T1~T4。

圖1 濟南市巖溶塌陷區地震反射波法勘探結果
可以看出T1~T4反射層位在68~95 m范圍內出現明顯異常變化,表現在T1層位出現空白反射或者反射能量較弱,T2~T4層位出現相位反轉、同相軸錯斷以及反射凌亂等現象,T2反射層位能量減弱,T4反射層位能量增強,推斷解釋此范圍內巖溶發育,并在78 m處布置驗證鉆孔。反演計算78 m處,T1層位埋深約為12.7 m,T2層位埋深約為34.5 m,T3層位埋深約為46.8 m。
驗證孔揭示,第四系地層埋深12.5 m,34.5 m為條帶狀灰巖和礫屑灰巖分界面,46.5 m以下為中風化灰巖。其中17.2~18 m發育巖溶空洞,21~23 m、27~33.5 m廣泛發育蜂窩狀溶孔,部分鉆孔揭露結果如圖2所示。淺層地震勘探反射波法推測結果與鉆孔揭露結果吻合程度較好,淺層地震勘探反射波法在基巖埋深較淺情況下的巖溶勘察中具有良好的應用效果。

圖2 濟南市巖溶塌陷區鉆探揭露結果
本實例的工作區位于山東省濟寧市,工作區上覆第四系地層,厚度數十米,下伏奧陶系灰巖、白云質灰巖以及泥質灰巖,厚度數百米,巖溶發育程度較弱。這種地質條件下進行淺層地震反射波法勘探需要重點考慮橫向分辨能力和縱向分辨能力,以區分規模較小的溶洞和裂隙[2]。
橫向分辨能力取決于菲涅爾帶半徑R的大小,通常水平分辨能力為菲涅爾帶的兩倍,

式中:h為界面埋深;v為該界面以上的地層速度;f為地震子波的主頻。
縱向分辨力指的是在垂直方向上能夠分辨的巖性單元的最小厚度。縱向分辨能力可以用時間域分辨率或者厚度分辨率兩種方式表示。理論上講,當地震子波為零相位子波時,地震勘探剖面中時間分辨率TR和厚度分辨率HR表示為:

式中:v為該界面以上的地層速度;f為地震子波的主頻。縱向分辨力取決于激發地震子波頻率和傳播速度,頻率越高、速度越小,其縱向分辨能力越高。
針對該工作區,采用全排列接收方式,2 m道間距,18次水平覆蓋觀測系統。圖3所示為該巖溶塌陷區地震反射波法時間剖面圖。反射波組中可以劃分3個反射層位,分別標注T1~T3。時間剖面中可以看到,T1、T2在82~123 m范圍內出現明顯異常變化,表現為T1反射層位出現弱反射并相位反轉,T2反射層位出現錯動和相位變化等現象,而T3反射層位連續可追蹤。推斷解釋此范圍內巖溶發育,并在105 m處布置驗證鉆孔。反演計算105 m處,T1層位埋深約為36.0 m,T2層位埋深約為58.3 m,T3層位埋深約為73.1 m。
驗證孔揭露,第四系地層埋深36.1 m,58.2 m為條帶狀灰巖和礫屑泥灰巖分界面,72.5 m以下為中風化灰巖。其中42.7 m和59.1 m發育巖溶空洞,角礫狀泥灰巖芯中廣泛發育溶腔、蜂窩狀溶孔。淺層地震勘探反射波法推測結果與鉆孔揭露結果吻合程度較高。
文章通過淺層地震勘探反射波法在巖溶塌陷勘察中的應用實例,給出了不同地質條件下,淺層地震勘探反射波法在巖溶塌陷勘察中的應用效果,分析認為淺層地震勘探反射波法在巖溶塌陷勘察中具有良好的應用前景。另外巖溶塌陷區規模不一、發育程度不同,尤其是溶洞、溶孔以及裂隙規模較小,淺層地震勘探反射波法很難分辨。巖溶塌陷區地層傾角較大時,也會對淺層地震勘探反射波法勘探效果造成不良影響。