地球物理中波動方程有限差分偏移模型應用及實驗對比

2014-03-26 03:27:24黃智海周艷玲

吉林地質 2014年2期

黃智海，周艷玲，趙峰，趙杰

吉林省地質工程勘察院，吉林長春 130001

1 模型介紹

通過淺層地震反射法獲得反射地震剖面，成功地揭示淺層地下結構，地震偏移基于波動方程，通過對反射層位置的校正以及對繞射能量的收斂來消除反射記錄中的失真。現在，偏移是地震數據處理過程的中心環節。不同偏移方法受到介質速度變化的不同限制，有些偏移方法要求介質是常速度或僅有縱向速度的變化，如Stolt偏移及相移法；有些偏移方法受到地層傾角的限制，如15°和45°有限差分偏移。為了對比傾角對方法的影響，建立具有不同傾角的層狀彎曲地層。根據方法對速度的要求，速度模型將分為常速度和縱向連續變化兩種。模型的反射層模型，速度模型的建立，Kirchhoff正演及速度分析過程由開源軟件Madagascar完成。模型及參數如下：

圖1反射層模型包含從0°到45°變化的傾角。正演時，道數為101道，炮檢距0.1 km，時間采樣總數1 301個，時間間隔0.004 s。深度延拓步長等于時間間隔。對于速度模型一，其速度為常量1.5 kft/s；速度模型二的速度是縱向線性變化的，初值為1.5 kft/s，變化梯度為0.25。實驗中加入了隨機噪聲。為了下文中方便進行處理解釋，從下到上地層分別標注為第一層到第四層。

與圖1進行對比可以看出，圖2中反射層的位置發生了偏移，地層傾角也發生了變化，且反射波在三、四層中蝴蝶結干涉嚴重。圖中第四層端點處的反射是因為吸收邊界的問題，本文有意留下，但它不在偏移處理的范圍之內。

圖3是連續速度模型是為了檢驗有限差分法及相移法的偏移效果，相移法要求速度沒有橫向變化，所以模型只有縱向上的速度變化。

從圖4中可看出除第一層，另外三層反射波位置都有不同程度的偏移；三，四層反射波出現蝴蝶結干涉。

圖1 反射層模型Fig.1 The re fl ective layermodel

圖2 常速度模型下合成的自激自收記錄Fig.2 Transmitting and receiving record composed in the constant velocitymodel

圖3 連續速度模型Fig.3 Continuous velocitymodel

圖4 連續速度模型下合成的自激自收記錄Fig.4 Transmitting and receiving record composed in the continuous velocitymodel

2 處理結果

應用15°有限差分偏移見圖5，常速度模型。偏移后一層沒有發生變化；二層的位置得到較好的校正，彎曲的程度與反射層模型中一樣；三，四層的位置也有不同程度的校正；蝴蝶結干涉有一定的分解，較緩的第三層中的干涉分解較好，第四層中仍有明顯的干涉。

應用45°有限差分偏移見圖6，常速度模型。

與15°有限差分偏移結果有明顯的不同，偏移的效果更徹底。同樣，第一層沒有變化；第二層反射波校正到真實地層位置，反射波很清晰；第三層和第四層反射波校正后位置發生了很大的改變，波形與原反射層模型很相近；蝴蝶結干涉受到不同程度的分解。

應用stolt偏移見圖7，常速度模型。

各層的反射波與有限差分偏移結果相比明顯要清晰一些，層的位置校正的也很好；干涉也得到較好的分解。

對比15°差分偏移和stolt偏移結果見圖8，為了對比不同偏移方法的偏移效果，將偏移后的剖面和原剖面進行疊加，以便清楚地進行對比。

左圖應用的是15°偏移，右圖應用的是Stolt偏移。可以看出兩個偏移方法的偏移結果具有很大差別。單從第二層來看，即對于較小的地層傾角，15°有限差分偏移和stolt偏移在常速度介質中的偏移效果是一樣的；但從三、四層的效果對比中可以看出，stolt偏移對反射層傾角的校正能力要好于15°有限差分偏移，且校正后的反射層也更加清晰；兩種方法都使干涉得到了分解，但對于三、四層這樣大傾角的地層，Stolt偏移的分解能力要好很多。

圖5 15°原剖面Fig.5 The 15°original pro fi le

圖6 15°有限差分偏移剖面Fig.6 The 15° fi nite-difference excursion pro fi le

圖7 45°原剖面Fig.7 The 45°original pro fi le

圖8 45°有限差分偏移剖面Fig.8 The 45° fi nite-difference excursion pro fi le

對比45°有限差分偏移和Stolt偏移結果見圖9。

對比兩圖，第二層的偏移效果是一樣的；對比三層和四層，45°有限差分偏移對反射層位置的校正程度要稍微大一些，但是反射層的清晰程度較Stolt偏移要差很多；從對蝴蝶結干涉分解的程度來看，45°有限差分偏移的效果要好一些，尤其是對第四層的干涉，干涉基本上已能夠完全分解，但是結合圖，可以推測，當反射層傾角繼續增大時，45°有限差分對干涉分解的能力將沒有stolt偏移好。

應用45°有限差分偏移見圖10，連續速度模型。

偏移后的剖面，第一，二層反射層更加清晰，第二層的位置也到了校正，彎曲變得圓滑；對三層和四層反射層的位置有明顯的校正；第三、四層的干涉有很好的分解。

應用相移法，連續速度模型見圖11。

相移法對傾斜反射層的位置校正的很好，干涉也得到了很好的分解。但是由于偏移速度采用的是均方根速度，使其偏移結果不夠徹底。

對比45°有限差分偏移和相移法結果見圖12

由于偏移速度用的是均方根速度，所以對偏移的效果有一定的影響，但是從圖中可以看出，有限差分法和相移法在縱向變速度介質中的偏移效果都很好。