席超強，王 偉，周文龍

（安徽理工大學(xué)地球與環(huán)境學(xué)院， 安徽淮南232001）

多道面波勘探排列參數(shù)的影響研究
——以雙層凸起地質(zhì)模型為例

席超強，王 偉，周文龍

（安徽理工大學(xué)地球與環(huán)境學(xué)院， 安徽淮南232001）

多道瞬態(tài)面波法（MASW）獲取地層剪切波速度方便快捷，目前已廣泛應(yīng)用在工程勘探領(lǐng)域，特別是在復(fù)雜地形條件下的應(yīng)用逐步增多;針對復(fù)雜地形條件下，數(shù)據(jù)采集至關(guān)重要，其中排列參數(shù)選擇是一個重要因素，以雙層凸起地質(zhì)模型為例，通過對不同偏移距和道間距的排列系統(tǒng)進行高精度有限差分數(shù)值模擬研究;對比分析數(shù)據(jù)頻率譜，討論不同偏移距與道間距對瑞雷波頻散分辨率特征的影響，為現(xiàn)場數(shù)據(jù)采集提供理論指導(dǎo);結(jié)果表明：在一定數(shù)量的通道數(shù)和橫向分辨率允許的情況下，大道間距是優(yōu)選,偏移距盡量大以避免近場效應(yīng)。

偏移距； 道間距； 排列參數(shù); 復(fù)雜地形

1 引言

多道面波分析法（MASW）利用沿地表水平方向傳播的Rayleigh型面波來提取地層的彈性性質(zhì)，目前已經(jīng)廣泛應(yīng)用于場地安全評價[1]、基巖埋深探查[2]、巖溶注漿檢測[3]、反演獲取剪切波速度剖面[4]。多道面波分析法需要分析瑞雷波的水平分量，由于面波從距離震源傳播一定距離才開始趨于水平穩(wěn)定傳播，低頻組分、具有較大的波長的面波需要較長的距離才能充分發(fā)育[5]。

目前對復(fù)雜地形下瑞雷波數(shù)值模擬主要通過使用交錯網(wǎng)格有限差分法。常用的有：利用交錯網(wǎng)格有限差分方法對粘彈性介質(zhì)起伏地表模型中地震傳播進行數(shù)值模擬[6]；模擬分析小凸起小凹陷模型條件下的瑞雷波場傳播特征[7]；采用高階交錯網(wǎng)格有限差分數(shù)值求解二維彈性波方程,然后對常見的層狀介質(zhì)和巖溶介質(zhì)模型分別進行模擬[8]；將笛卡爾坐標系的各向異性介質(zhì)彈性波方程和自由邊界條件變換到曲線坐標系中,采用一種穩(wěn)定的、顯式的二階精度有限差分方法離散介質(zhì)中的彈性波方程,有效地處理自由地表邊界條件[9]；模擬起伏地表條件下的瑞雷波全波場正演，通過坐標變換的方式實現(xiàn)起伏地表自由邊界條件[10]；進行了傾斜和塌陷地形的彈性波場交錯網(wǎng)格有限差分數(shù)值模擬[11]。以上對瑞雷波在數(shù)值模擬上的分析主要集中在完善對復(fù)雜地形下瑞雷波有限差分模擬的實現(xiàn)，并未對排列參數(shù)對瑞雷波頻散特征影響進行分析。

多道面波分析法的關(guān)鍵在于準確提取多模式的頻散曲線。因為頻散曲線與地層速度性質(zhì)密切相關(guān)。對復(fù)雜地形中的MASW排列參數(shù)進行研究為確保在數(shù)據(jù)處理階段能獲得最準確的頻散曲線。本文將分析偏移距、道間距對復(fù)雜地形條件下頻散譜分辨率的影響。

2 有限差分數(shù)值模擬理論

2.1 彈性波動方程理論

假設(shè)地震波傳播的介質(zhì)為各向同性彈性介質(zhì)。根據(jù)彈性力學(xué)三個基本方程（運動微分方程、幾何方程、本構(gòu)方程）進行數(shù)據(jù)變化得到彈性場的位移方程與彈性方程，二維各向同性介質(zhì)一階速度應(yīng)力波動方程為：

（1）

其中Vx、Vz為速度分量,τzx、σxx、σzz為應(yīng)力分量；λ、μ為拉梅常數(shù)。

2.2 交錯網(wǎng)格有限差分數(shù)值解

采用交錯網(wǎng)格法（時間M=2，空間N=12）對公式（1）進行運算。取x=iΔx、z=jΔz、t=kΔt，其中i，j，k分別表示空間x方向、空間z方向與時間t的網(wǎng)格點數(shù)，Δx、Δz、Δt分別為空間x方向、空間z方向的空間步長與時間t的時間步長。

在運用交錯網(wǎng)格差分法求解一階彈性波動時，Vx、Vz值于半節(jié)點處運算，τxz、σxx、σzz值需在網(wǎng)格節(jié)點處計算，公式（1）第1個方程計算格式如下：

同理，可得公式（1）其余4個方程的2階時間12階空間差分近似方程：

2.3 波動方程的吸收邊界條件

對地震波場進行數(shù)值模擬時，定義的人工邊界會產(chǎn)生干擾波如邊界反射，影響數(shù)據(jù)計算效果，因此采用完全匹配層（PML）吸收邊界條件來消除自定義人工邊界產(chǎn)生的干擾波對波場模擬運算的影響。

3 模型試驗

使用雙層凸起地形地質(zhì)模型實現(xiàn)對復(fù)雜地形的數(shù)值模擬。模型尺寸：100m×100m，網(wǎng)格點數(shù)為1 000×1 000，模型上部地層為橫波速度200m/s的低速層，模型下部地層為橫波速度400m/s的高速層。地質(zhì)體模型參數(shù)見表1。

表1 雙層凸起地質(zhì)模型參數(shù)表

圖1 雙層凸起地質(zhì)模型示意圖

在不同偏移距對比試驗中，炮點分別于20m，14m，10m處進行激發(fā)，使用的震源為主頻20Hz的Ricker子波，相應(yīng)偏移距2m，8m，12m。共使用26個檢波器，各檢波器通道間距2m，檢波器位置x=22～72m。數(shù)據(jù)采樣時間600ms，采樣間隔0.2ms，采樣點數(shù)3k，如圖2所示。

圖2 不同偏移距對比試驗排列示意圖

在道間距對比試驗中炮點于14m處激發(fā)，偏移距8m，共進行三組模擬試驗：道間距為0.5m時，檢波器范圍為x=22～34.5m；道間距為1m時檢波器分布范圍為x=22～47m，道間距為2m時，檢波器分布范圍為x=22～72m，數(shù)據(jù)采樣時間600ms，采樣間隔0.2ms，采樣點數(shù)3k，三組模擬試驗排列如下圖3所示。

圖3 不同道間距對比試驗排列示意圖

4 頻譜特征分析

本文使用GeogigaSurfaceWavePlus8.0 模塊進行面波數(shù)據(jù)處理。載入有限差分運算所得的地震記錄，通過視窗可對面波進行干擾波切除與面波視速度測量。圖4為偏移距為8m，道間距為2m的排列面波單炮記錄。從圖中能夠看出“掃帚”狀面波信號，信號沒有受到人工邊界的反射影響，模擬數(shù)據(jù)質(zhì)量較高。

圖4 模擬運算面波單炮記錄

SurfaceWavePlus模塊提供三種頻散計算方法：相移法、F-K法、自相關(guān)法（SPAC），本文使用F-K法進行頻譜運算，速度分析范圍0～500m/s，頻率分析碼范圍0～80Hz。進行地震道平衡與頻譜平衡處理以突出能量顯示效果。聚焦系數(shù)為5。聚焦系數(shù)范圍為1～5，系數(shù)越大頻散收斂程度越高。統(tǒng)一頻散計算參數(shù)，對比不同排列參數(shù)對面波頻散的影響。

4.1 不同偏移距與頻散的影響

圖5（a）為偏移距2m時頻散譜，25～31Hz、50～64Hz范圍內(nèi)頻率不連續(xù)有跳點出現(xiàn)。這是由于上下兩組地層速度之間存在差異，面波能量向高階發(fā)育的現(xiàn)象。7～25Hz頻率范圍為面波基階能量；31～50Hz頻率范圍為第一高階能量；64～80頻率范圍為第二高階能量范圍。圖5（b）為偏移距8m時頻散譜，在25Hz、63Hz處頻散能量產(chǎn)生跳點。相對于圖5（a）能量譜完整性較好，散點較少。7～25Hz頻率范圍為面波基階能量、25～63Hz頻率范圍第一高階能量，63～80為第二高階能量。圖5（c）為偏移距12m時頻散譜，頻譜圖能量整體收斂程度低于圖5（a）和圖5（b）中能量譜。各階能量團相連，分離識別難度大。7～26Hz為面波基階能量，26～58Hz為第一高階能量，58～60Hz為第二高階能量。

通過比較不同偏移距的頻散譜，有如下認識：圖5（a）三種不同階的能量分離清晰，但有不連續(xù)散點產(chǎn)生，因為排列距震源距離較近，受近場效應(yīng)影響，高階能量比較發(fā)育，各階結(jié)合處受影響形成不連續(xù)散點現(xiàn)象。圖5（b）三種不同階能量分離清晰，在各階能量交界處有散點產(chǎn)生，但是散點主要趨向高階能量，說明在偏移距為8m時，基階波已經(jīng)充分發(fā)育，受高階影響不大，通過拾取三階頻散曲線，可以進行后續(xù)聯(lián)合反演。圖5（c）三種能量未能清晰分離，各階結(jié)合處有所相連，連接處的散點趨于低階能量，說明波經(jīng)過一定距離傳播后，高階波能量減弱，基階波能量增強，基階波能量趨于主導(dǎo)地位。

隨著偏移距的增加，近場效應(yīng)影響減弱，基階波能量會完全占據(jù)主導(dǎo)地位，可以使用基階波單一反演算法進行處理。在實際施工過程中，大偏移距導(dǎo)致震源如錘擊能量在長傳播過程中的更大衰減，以及周圍環(huán)境如施工機械噪音對數(shù)據(jù)的影響更加明顯。通過選擇適中的偏移距可獲取質(zhì)量較高的頻散譜，可以獲取準確的頻散曲線。

（a） 偏移距2m（b） 偏移距8m（c） 偏移距12m

圖5 不同偏移距頻散對比示意圖

4.2 不同道間距對頻散的影響

通過上組試驗的頻譜對比結(jié)果，固定炮點位置，確定偏移距為8m，改變排列的道間距來比較不同道間距對頻散譜的影響。

圖6（a）為道間距為0.5m時頻散譜，通道排列長度為12.5m，頻率譜有多處跳點，7～20Hz范圍為面波基階能量，36～55Hz頻率范圍為第一高階面波能量，60～80Hz為第二高階面波能量。圖6（b）為道間距為1m時頻散譜，7～24Hz頻率范圍為面波基階能量，24～58Hz頻率范圍為第一高階面波能量，58～80Hz頻率范圍為第二高階面波能量。圖6（c）為道間距為2m的頻散譜，7～25Hz頻率范圍為面波基階能量、25～63Hz頻率范圍第一高階面波能量，63～80為第二高階面波能量。

通過比較不同道間距的頻散譜，有如下認識：圖6（a）到圖6（c）的三組道間距從小到大，能量收斂性也逐漸增強。圖6（a）頻散質(zhì)量較差，10Hz以下低頻部分分辨率較低原因在于排列長度太短，根據(jù)采樣定理短排列無法采集低頻大波長的面波組分，20～30Hz這一組分基階能量受第一高階能量影響嚴重，高階能量雖有部分缺失但從形態(tài)特征均能判斷所在模態(tài)。圖6（b）頻散質(zhì)量一般，收斂性優(yōu)于圖6（a），三種不同階能量分離清晰，基階與第一高階、第一高階與第二高階結(jié)合處有能量散點分布，受低階、高階共同影響，速度介于理論低階與理論高階之間。圖6（c）頻散譜能量收斂性最佳，各階能量分離清晰，各階交界處有散點產(chǎn)生，但是散點主要趨向高階能量，各階能量團兩端都受相鄰階能量團的影響。隨道間距的增大，受相鄰階能量團的影響也在逐漸減弱。

隨著道間距增加，排列長度增加，能夠有效分離基階、高階能量，可拾取受相鄰階能量影響較小、趨于理論值的頻散曲線。多道面波分析技術(shù)獲取的通道排列范圍內(nèi)的平均速度。大道間距、長排列導(dǎo)致橫向分辨率降低。通過選擇合適的道間距可獲取質(zhì)量較高的頻散譜，可以獲取準確的頻散曲線。

（a） 道間距0.5m（b） 道間距1m（c） 道間距2m

圖6 不同道間距頻散對比示意圖

5 結(jié)論

本文以雙層凸起介質(zhì)模型模擬復(fù)雜地形，通過對不同偏移距、不同道間距的排列形式進行復(fù)雜地形條件下的數(shù)值模擬。使用F-K法對數(shù)據(jù)進行頻散運算，對頻散譜比較分析得到如下結(jié)論與認識：

隨偏移距增加，近場效應(yīng)影響減弱，基階波能量在波場中逐漸占據(jù)主導(dǎo)地位，影響高階波分離，選擇合適的偏移距，拾取各階分離較好的頻散曲線。

隨道間距增加、排列長度增大，頻散譜能量趨于收斂，各階能量相互影響減弱，便于拾取精度較高的頻散曲線，但橫向分辨率所有降低。實際勘探過程中需要謹慎選擇合適道間距。

本文僅針對雙層凸起介質(zhì)模型進行復(fù)雜地形下排列參數(shù)分析對比試驗，模型較為單一，后續(xù)須對其他復(fù)雜地形進行數(shù)值模擬并分析頻譜特征，從理論上對復(fù)雜地形面波勘探進行指導(dǎo)，以獲取高質(zhì)量的面波數(shù)據(jù)記錄。

[1] 崔建文,喬森,樊耀新.瞬態(tài)面波勘探技術(shù)在工程地質(zhì)中的應(yīng)用[J].巖土工程學(xué)報,1996,18（3）:35-40.

[2] 楊建強,錢榮毅,鄧桂林,等.瑞雷面波折射層析及井間CT勘探在基巖地基勘察中的應(yīng)用研究[J].工程地球物理學(xué)報,2012,9（2）:231-237.

[3] 林承灝,張平松,郭立全,等.利用多道面波勘探技術(shù)調(diào)查路基巖溶地質(zhì)條件[J].工程地球物理學(xué)報,2012,9（5）:595-599.

[4] 胡家富,楊潤海.利用Rayleigh波反演淺土層的剪切波速度結(jié)構(gòu)[J].地球物理學(xué)報,1999,42（3）:393-400.

[5] 楊成林.瑞雷波法勘探原理及其應(yīng)用[J].物探與化探,1989,13（6）:465-468.

[6] 董良國.復(fù)雜地表條件下地震波傳播數(shù)值模擬[J].勘探地球物理進展,2005,28（3）:187-194.

[7] 熊章強,南坤,唐圣松.起伏地表二維瑞雷面波正演模擬及波場分析[J].物探化探計算技術(shù),2010,32（5）:470-475.

[8] 鄒延延,徐義賢,沙椿,等.彈性介質(zhì)中垂直地震剖面有限差分正演[J].地球物理學(xué)進展,2011,26（6）:1959-1967.

[9] 蘭海強,劉佳, 白志明.VTI介質(zhì)起伏地表地震波場模擬[J].地球物理學(xué)報,2011,54（8）:2072-2084.

[10]江凡.起伏地表條件下坐標變換的瑞雷波場正演模擬[J].福建工程學(xué)院學(xué)報,2013,11（1）:14-18.

[11]琚長輝,王偉,王祥春,等.基于起伏地形的瑞利面波有限差分數(shù)值模擬[J].東華理工大學(xué)學(xué)報（自然科學(xué)版）,2015,38（3）:319-324.

Study on the Influence of Arrangement Parameter of Multi-channel Surface Wave Exploration in Complex Terrain——Taking Double-layer Raised Geologic Model as an Example

XI Chaoqiang,WANG Wei,ZHOU Wenlong

（SchoolofEarthandEnvironment,AnhuiUniversityofScienceandTechnology,HuainanAnhui232001）

Multi-channel transient surface wave method （MASW） is convenient and rapid in acquiring strata shear wave velocity, and has been widely used in engineering exploration, especially in complex terrain. In the case of complex terrain, data acquisition is very important. The selection of arrangement parameter is an important factor. Taking double-layer raised geologic model as an example, this paper studied the high precision finite difference numerical value through the arrangement system of different offset and track pitch. The influence of different offset and track pitch on the resolution characteristics of Rayleigh wave dispersion rate is discussed by comparing the data frequency spectrum, which provides theoretical guidance for field data acquisition. The results show that the aisle spacing is preferable when a certain number of channels and lateral resolution are allowed; the offset should be as large as possible to avoid the near-field effect.

offset distance; track pitch; acquisition parameters; complex terrain

2016-12-07

席超強（1993-）,男，安徽亳州人，碩士，主要研究方向為工程地球物理勘探。

P631

B

1671-4733（2017）01-0022-05

10.3969/j.issn.1671-4733.2017.01.008