張會(huì)卿，周宗良，曹國(guó)明，燕 云

(中國(guó)石油大港油田公司勘探開發(fā)研究院 天津300280)

0 引 言

大港 X開發(fā)區(qū)位于北大港構(gòu)造帶東部，構(gòu)造異常破碎，主力層明化鎮(zhèn)組受曲流河河道遷移擺動(dòng)影響，儲(chǔ)層具有單砂體厚度薄、橫向變化快等特點(diǎn)，均屬于 1/4地震波波長(zhǎng)以下的薄砂體[1]，無(wú)法從地震剖面上直接識(shí)別。為了解決上述問題，將“兩寬一高”地震勘探理念引入該區(qū)復(fù)雜油藏薄砂體精細(xì)預(yù)測(cè)中。本次依托高密度地震資料，創(chuàng)新提出基于小入射角道集疊加的最大熵地震頻譜分解預(yù)測(cè)方法，綜合利用小入射角道集疊加處理、薄層調(diào)諧分析、井震多參數(shù)交會(huì)等方法，破解了復(fù)雜地表?xiàng)l件下地震采集頻率低的業(yè)界難題，建立了調(diào)諧頻率、調(diào)諧厚度和砂巖厚度的關(guān)系，實(shí)現(xiàn)了單一薄砂體的準(zhǔn)確預(yù)測(cè)。

1 小入射角疊加處理

1.1 技術(shù)原理

俞壽朋[2]指出，地震勘探分辨率隨炮檢距增大而降低。由于地下介質(zhì)是一組或多組由不同斷層、裂縫等組成的地層，地震波在傳播過程中要發(fā)生不同程度的衰減。總體而言，同一深度下，入射角越小，炮檢距越小，地震波能量和頻率的衰減程度越低。

下面論證地震波入射角與分辨率之間的理論關(guān)系[3]。假設(shè)存在一水平地層，地層頂面、底面埋深分別為H1、H2，地層厚度為h，地層頂、底面地震波入射角分別為θ1、θ2，地層上覆介質(zhì)速度為 V1，地層速度為 V2。由于 H2、H1h，則有：

地層頂界面雙程旅行時(shí) t1和底界面雙程旅行時(shí)t2可分別表達(dá)為：

由于 V1≈V2=V，則地層頂?shù)捉缑骐p程旅行時(shí)平方差為：

因?yàn)镠?h，所以上式可整理為：

根據(jù)Rayleigh準(zhǔn)則，雙程旅行時(shí)差為1/2周期時(shí)達(dá)到地震縱向分辨率極限，則有：

此時(shí)可分辨地層厚度近似為：

根據(jù)公式 9可得，對(duì)于某一時(shí)間深度(t0)處的地層，入射角越小，則可分辨儲(chǔ)層厚度越小，即分辨率越高。因此，選取小入射角道集進(jìn)行疊加處理，可提高地震分辨率。

常規(guī)地震資料由于覆蓋次數(shù)不足，難以確保小入射角道集的有效疊加處理，而高密度地震資料則由于地震采樣點(diǎn)密集、覆蓋次數(shù)高的特點(diǎn)，破解了這一難題。

1.2 應(yīng)用實(shí)例

本次應(yīng)用“兩寬一高”地震數(shù)據(jù)體覆蓋次數(shù)為360次。其中，0～10°入射角覆蓋次數(shù)為 32～37次，0～20°入射角覆蓋次數(shù)為 140次，大于 25°的大角度入射角覆蓋次數(shù)為 78～86次。本次以確保入射角較小和覆蓋次數(shù)較高為原則，選取 0～20°入射角道集進(jìn)行疊加處理。

圖1為全疊加剖面與小入射角道集疊加剖面對(duì)比實(shí)例。圖1c顯示，lq5-9、lq7-8等井 NmⅢ-4發(fā)育2層薄砂體NmⅢ-4-1、NmⅢ-4-2，這2層砂體從全疊加剖面(圖1a)難以識(shí)別，但在小入射角道集疊加剖面中(圖1b)卻能夠非常清晰地識(shí)別出來，表明通過選取小入射角道集進(jìn)行疊加處理，可以提高地震垂向分辨能力。

圖1 NmⅢ-4小層地震剖面對(duì)比Fig.1 Comparison of NmⅢ-4 formation seismic sections

2 最大熵譜分解

2.1 技術(shù)原理

20世紀(jì)90年代中期，Partyka等[4]首次提出譜分解概念，即利用時(shí)頻分析方法，通過數(shù)學(xué)變換將地震信息從時(shí)間域轉(zhuǎn)換到頻率域，進(jìn)而精確表征地震信號(hào)中每一時(shí)間點(diǎn)的頻率成分[5]，繼而衍生出離散傅里葉變換、最大熵等多種譜分解算法。

最大熵在 1967年由 J.P.伯格提出，用來衡量信息的不確定性；之后 Aydm Akan[6]等將最大熵引入地震信號(hào)分析，用于計(jì)算地震信號(hào)的功率譜。最大熵譜分解原理如下：地震記錄x(t)為一個(gè)平穩(wěn)的離散時(shí)間序列，其熵的表達(dá)式為[7]：

式中：Hm為熵，J/mol·K；t＝1，2，...，N+1，ms；T(N)為托布里茲形式的信號(hào)矩陣，dB。若 r(0)，r(1)，...，r(N)值已知，則可重復(fù)利用遞推算法求取 r(N+1)，r(N+2)，...的值。

地震數(shù)據(jù)道x(t)可用一個(gè) K階自回歸差分方程模擬：

相應(yīng)的地震信號(hào)功率譜P(f)為：

式中：x(t)為地震數(shù)據(jù)道，CDP's/cm；e(t)為預(yù)測(cè)誤差，CDP's/cm；ΔδΚ為預(yù)測(cè)誤差的方差；P(f)為地震信號(hào)功率譜，dB2/Hz；ak為預(yù)測(cè)誤差濾波器的系數(shù)，無(wú)量綱；k 為常數(shù)，k＝1，2，...，K；f 為角頻率，rad/s；Δt 為采樣間隔，ms。

最大熵譜分解確保地震記錄在分析時(shí)窗內(nèi)的自相關(guān)函數(shù)符合已知數(shù)據(jù)要求，而使時(shí)窗外自相關(guān)函數(shù)的不確定性最大。相比傅里葉變換要求地震信號(hào)在(-∞，+∞)之間取值，最大熵減少了因時(shí)窗長(zhǎng)度有限而產(chǎn)生的誤差。最大熵譜分解可生成調(diào)諧體。

調(diào)諧體預(yù)測(cè)薄砂體厚度源于 Widess模型：Widess楔形模型證實(shí)，1/4波長(zhǎng)以下薄層的頂、底反射界面發(fā)生干涉[1]；Kallweit得出薄層等于 1/4波長(zhǎng)時(shí)將產(chǎn)生調(diào)諧效應(yīng)，且調(diào)諧點(diǎn)處振幅值最大[8-9]。黃緒德[10]論證了薄層振幅譜的周期性陷頻特征，其陷頻周期和頻率互為倒數(shù)。邊立恩等[11]論證了薄層調(diào)諧頻率與調(diào)諧厚度關(guān)系，推導(dǎo)出調(diào)諧厚度估算公式：

式中：ΔZ為調(diào)諧厚度，m；v為目的層速度，m/s；f0為調(diào)諧頻率，Hz。

2.2 應(yīng)用實(shí)例

利用最大熵法制作調(diào)諧體并得到調(diào)諧頻率，調(diào)諧體時(shí)窗依據(jù)砂巖時(shí)間厚度確定為 30ms。東西向調(diào)諧體剖面中，A1井調(diào)諧頻率約為 25Hz，向兩翼調(diào)諧頻率急劇升高，反映調(diào)諧厚度減薄(圖2)；南北向調(diào)諧體剖面中，調(diào)諧頻率始終穩(wěn)定在 25～30Hz，反映調(diào)諧厚度基本不變(圖3)。調(diào)諧體剖面反映巖性邊界更加直觀，且與單頻體切片對(duì)應(yīng)性很好。

圖2 A1井-A1-8井最大熵調(diào)諧體剖面Fig.2 Maximum entropy tuning cube profile of wells A1-A1-8

圖3 A1-3井-A1-7井最大熵調(diào)諧體剖面Fig.3 Maximum entropy tuning cube profile of wells A1-3-A1-7

馬躍華等[12]利用調(diào)諧頻率計(jì)算調(diào)諧厚度，單砂體厚度誤差為15%～83%。誤差源于以下因素[13]：①任何算法均存在誤差，造成調(diào)諧頻率的求取不準(zhǔn)確；②由于單砂體厚度的空變性，難以用一個(gè)確定的調(diào)諧體時(shí)窗確保只有目的層地震信號(hào)參與計(jì)算；③地震信號(hào)中低頻噪聲與有效的低頻信息相夾雜，影響低頻端厚度估算；④受測(cè)不準(zhǔn)原理、測(cè)井-地震觀測(cè)頻率不同、儲(chǔ)層物性參數(shù)不同等因素的影響[14-15]，速度參數(shù)難以準(zhǔn)確獲取。因此，調(diào)諧厚度不能直接表征單砂體厚度，而應(yīng)理解為相對(duì)厚度。在地層速度較穩(wěn)定的地區(qū)，可通過統(tǒng)計(jì)完鉆井資料，建立調(diào)諧厚度與單砂體厚度相對(duì)關(guān)系，從而求取井間單砂體厚度。

本次統(tǒng)計(jì)了 15口完鉆井的儲(chǔ)層參數(shù)，利用聲波時(shí)差獲取單井層速度，代入調(diào)諧厚度估算公式，計(jì)算得到調(diào)諧厚度。通過建立多井調(diào)諧厚度和單砂體厚度關(guān)系，證實(shí)二者呈線性正相關(guān)(圖4)，關(guān)系式為：

式中：Z為單砂體厚度，m；ΔZ為調(diào)諧厚度，m。

3 應(yīng)用效果

應(yīng)用公式(15)完成單砂體厚度預(yù)測(cè)(圖5)，在此基礎(chǔ)上優(yōu)選構(gòu)造有利區(qū)部署2口新井A1-1-1和A1-1-2，新井成功鉆遇 NmⅢ-8-2儲(chǔ)層，且單砂體厚度與預(yù)測(cè)厚度吻合率達(dá) 90%以上。鉆井結(jié)果與預(yù)測(cè)結(jié)果吻合，證明基于最大熵譜分解方法的單砂體定量預(yù)測(cè)方法在本區(qū)適用。

圖4 A1井區(qū)單砂體厚度與調(diào)諧厚度關(guān)系Fig.4 Relation between thickness of single sand body and tuning thickness in well A1 area

圖5 A1井區(qū)NmⅢ-8-2砂巖厚度等值圖Fig.5 Contour map of NmⅢ-8-2 sandstone thickness in well A1 area

4 結(jié) 論

①地震縱向分辨率與地震波入射角近似呈正比，基于“兩寬一高”的小入射角道集疊加處理能夠有效提高地震資料垂向分辨能力。

②應(yīng)用最大熵譜分解可得到薄層調(diào)諧體，利用調(diào)諧頻率可計(jì)算調(diào)諧厚度。

③單一薄砂體厚度和調(diào)諧厚度呈線性正相關(guān)關(guān)系，據(jù)此可實(shí)現(xiàn)井間單一薄砂體厚度預(yù)測(cè)。