基于振幅和頻率的砂體厚度校正新方法及其應用
——以渤海南部X油田為例

薛明星，劉傳奇，高 磊，陳文雄，劉 歡

(中海石油(中國)有限公司天津分公司渤海石油研究院，天津塘沽 300452)

基于振幅和頻率的砂體厚度校正新方法及其應用
——以渤海南部X油田為例

薛明星，劉傳奇，高 磊，陳文雄，劉 歡

(中海石油(中國)有限公司天津分公司渤海石油研究院，天津塘沽 300452)

傳統的砂體厚度圖校正方法基于砂體頂底反射視時間厚度，然而該屬性僅在砂體厚度大于1/4波長時與真厚度近似具有線性關系，當砂體厚度小于1/4波長時將產生較大的校正誤差?；谛钅Ｐ驼?，使用反射振幅與中心頻率構建了一種在3/7波長范圍內與砂體真厚度具有線性關系的新屬性 ，給出了一種全新的砂體厚度校正方法。新方法思路簡潔，易于實現，適用于生產中大多數橫向較均勻砂體的厚度校正問題。應用實例表明，新方法校正結果精確合理，顯著優于傳統方法。

砂體厚度校正；視厚度；反射振幅；中心頻率

在河流相油田開發過程中，沉積砂體厚度的計算精度直接影響油田儲量的確定以及開發井井位的優化。傳統的砂體厚度校正方法基于視厚度進行計算，假定其與真厚度具有線性關系，視厚度能夠反映真厚度的變化趨勢，然而眾所周知，該關系只在砂體厚度大于1/4波長的條件下近似成立，當砂體厚度小于1/4波長時，視厚度基本保持不變，此時不能反映真厚度的變化趨勢[1-5]。

油田實際開發過程中涉及大量整體或局部為薄層的含油砂體，對于此類情況傳統方法必然產生較大誤差，在無井鉆遇處得到的砂體厚度可能與實際厚度相去甚遠。合理的校正方法要求校正的屬性應該與砂體真厚度在較寬的厚度范圍內具有簡單的線性關系。研究表明，可以利用振幅與頻率構造具有這一特征的屬性。本文利用雷克子波合成記錄，給出了聯合利用砂體反射振幅與中心頻率進行厚度校正的新思路[6-9]。

1 基本原理[7-13]

為了探討與真厚度具有線性關系的新屬性，采用雷克子波合成地震記錄，考慮如圖1a所示頂底反射系數等幅異號楔狀砂體模型，利用褶積關系得到合成地震記錄(圖1b，波速2 300 m/s,子波主頻40 Hz)，由圖1可提取楔狀模型振幅與時間厚度曲線(圖2)，振幅在14.4 m(1/4波長)處發生調諧。傳統砂體厚度校正方法利用各井點實鉆厚度值對砂體視厚度進行校正得到最終校正結果，然而圖2表明，當砂體厚度小于1/4波長時，視厚度不能反映真厚度的變化趨勢。中心頻率與時域提取的振幅屬性關系如圖3。

圖1 楔狀模型(a)及其合成地震記錄(b)

由圖3可知，楔狀模型反射振幅曲線隨厚度的增加具有一個極大值點，即1/4波長調諧點，而頻率曲線具有一個極小值點。當厚度小于1/4波長時，振幅與厚度近似具有線性關系，頻率變化則較為緩慢；當厚度大于1/4波長且小于頻率極值點時，振幅下降，同時頻率值也快速降低。因此，可嘗試利用頻率曲線的特點對振幅曲線進行改造，使新曲線在比1/4波長更寬的范圍內與厚度具有線性關系。最直接的改造形式為采用A/F與厚度的關系曲線(A為振幅，F為頻率)，曲線形態如圖4所示，由該圖可知，曲線極值已經向右移動，對應厚度超過1/4波長，但移動量較小且極值點左側曲線線性特征不理想，為此可將變換式修正為A/Fn，加強頻率對曲線的調整作用。圖5為n分別取3，3.5，4時變換式A/Fn的形態。經過反復嘗試，最終確定當n取3.5時，新曲線極值點左側部分與厚度的線性關系最理想，且極值點對應厚度值約為3/7波長(圖5實例中極值點厚度為25 m，3/7波長為24.6 m)。由于單砂體儲層研究以波阻抗反演或90°相移方法為基礎，而一般這兩種方法所能刻畫砂體的厚度大多小于3/7波長，因此采用變換式表達的新屬性對砂體厚度進行校正，能夠滿足大多數生產實際的需要。新的厚度校正方法需要提取砂體頂面(或底面)反射振幅及中心頻率計算新屬性 ，再利用各井點實鉆砂體厚度值對該屬性進行校正，即獲得厚度最終校正結果。

圖2 楔狀模型振幅、時間與厚度關系

圖3 楔狀模型振幅、中心頻率與厚度關系

圖4 楔狀模型振幅、中心頻率與厚度關系

圖5 n取不同值時A/Fn極大值點左側曲線形態

2 應用實例

以渤海南油X油田F25-1575砂體為例，該含油砂體所在層段1/4波長為14 m，砂體先后被F25、F21兩口井鉆遇(圖6)。兩口井處砂體時間厚度均為12 ms，F25井實鉆厚度9.2 m，按照傳統厚度校正方法預測，F21井處砂體厚度也應為9.2 m，但實鉆表明F21井處砂體厚度僅為5.6 m?？紤]振幅頻率信息，F25井及F21井處的振幅分別為-30 446 dB與-21 262 dB,中心頻率分別為33 Hz與34 Hz,按照新方法計算A/F3.5，分別為0.147 5與0.092 8，由此預測F21井處砂體厚度為5.8 m，與實鉆結果十分接近。根據實鉆信息采用傳統方法對砂體厚度進行校正如圖7a，將其與砂體剖面圖(圖6)對比發現，剖面上砂體最左端高部位(區域A)時間厚度達20毫秒，但校正結果厚度僅為8～9 m，顯然不合理，這是受F21井處較薄的實鉆厚度影響的結果。利用振幅頻率信息采用新方法對砂體厚度進行校正，如圖7b，此時砂體左端高部位校正結果可達12～14 m，該厚度與砂體地震反射特征更協調，計算結果準確合理。

圖6 F25-1575砂體過井90度相移資料剖面

圖7 F25-1575砂體90度相移振幅屬性與厚度等值線疊合

3 結論

為了克服基于頂底反射視厚度進行砂體厚度校正的傳統方法存在的缺陷，本文通過反射振幅和中心頻率構造了一種在3/7波長范圍內與砂體真厚度具有線性關系的新屬性 ，并給出了一種全新的砂體厚度校正方法，該方法簡潔、易于實現。應用實例表明，當薄層存在時，新方法校正結果更加精確，顯著優于傳統方法，值得推廣應用。當然，新屬性的構建基于均勻的楔狀體模型，要求砂體橫向較均勻、物性變化小，這是該方法使用過程中應該注意的地方。對于橫向物性變化較大的砂體厚度的校正問題，還需要開展進一步的研究工作。

[1] 孫魯平，鄭曉東，首皓，等.薄層地震峰值頻率與厚度關系研究[J].石油地球物理勘探，2010，45(2)：254-259.

[2] 王云專，郭雪豹，邢小林，等.薄層峰值頻率特征分析[J].地球物理學進展，2013，28(5)：2515-2523.

[3] 張海燕，李慶忠.幾種常用解析子波的特性分析[J].石油地球物理勘探，2007，42(6)：651-657.

[4] 李雪英，李東慶，白詩緣.薄層研究方法綜述[J].地球物理學進展，2014，29(5)：2197-2203.

[5] 黃真萍，王曉華，王云專.薄層地震屬性參數分析和厚度預測[J].石油物探，1997，36(3)：28-38.

[6] 汪恩華，賀振華，李慶忠.薄儲層厚度計算新方法探索[J].物探化探計算技術，2001，23(1)：22-25.

[7] 寧松華，汪勇.基于振幅屬性分析法的楔狀砂體厚度預測研究[J].石油天然氣學報，2010，32(5)：90-94.

[8] 云美厚，丁偉.地震子波頻率淺析[J].石油物探，2005，44(6)：578-581.

[9] 張兆輝，李勝軍，張忠義，等.目標砂體厚度精細預測技術在歧南地區的應用[J].地球物理學進展，2014，29(5)：2314-2319.

[10] 袁子龍，楊冰，王建國.薄層、薄互層地震反射時間域與頻率域正演模擬研究及應用[J].石油物探，1996，35(3)：14-20.

[11] 鄒華勝，寧書年，楊峰，等.支持向量機在儲層厚度預測和計算中的應用[J].地球物理學進展，2007，22(5)：1597-1601.

[12] 王云專，王曉華，黃真萍.時域頻域多參數薄層厚度預測[J].大慶石油學院學報，1996，20(3)：16-20.

[13] 李雪英，陳樹民，王建民，等。薄層時頻特征的正演模擬[J].地球物理學報，2012，55(10)：3410-3419.

編輯：趙川喜

1673-8217(2017)03-0070-03

2016-10-19

薛明星，工程師，1983年出生，2008年畢業于成都理工大學物探專業，現從事油田開發地震研究工作。

TE112

A