基于地震屬性特征的河流相疊置砂巖儲層預測方法

井涌泉 欒東肖 張雨晴 陳 飛 范廷恩 胡光義

(中海油研究總院，北京 100028)

1 引言

地震屬性分析技術的研究和應用貫穿于地震資料儲層解釋(儲層預測[1-5]、巖性解釋[6-8]、流體預測[9-11])的各個環節，目前振幅屬性分析[12-14]、譜分解[15-17]、地震相分析[18-21]以及多地震屬性聚類分析[22，23]等技術在預測小斷層、儲層及其物性參數等方面起著重要作用。地震屬性分析技術通過建立儲層與地震屬性間的相關關系解決儲層預測問題，因此，在特定地質條件下，明確地震屬性的地質意義，優選地震屬性或者地震屬性組合則是實現儲層預測的關鍵問題。

渤海M油田明化鎮組(Nm)下段發育的河流相砂巖是主要的油氣儲集體，單砂層較薄(3～5m)，地震資料調諧厚度約為12m，地震反射通常反映垂向多期砂巖儲層疊合特征。因此，筆者以單井砂體不同垂向疊合模式地震反射特征的統計分析為基礎，優選與河道砂體疊合特征相關的地震屬性，并應用地震屬性聚類分析技術預測實際砂巖儲層，展現出較好的應用前景。

2 砂巖儲層疊合模式地震響應

根據物源供給與可容納空間比值(A/S)，將研究區新近系曲流河相沉積復合砂體疊合模式分為堆疊型、側疊型以及孤立型[24-27]。本區M油田Nm下段河道砂體疊合模式主要為側疊型和孤立型，堆疊型不甚發育[23，24]。根據該劃分模式，通過單井井震標定及聯井小層對比各疊合模式的地震響應特征。

2.1 側疊型地震響應

在基準面上升早期，物源供給與可容納空間比A/S較小，但此時可容納空間增大，河道彎曲度增大，單套砂體規模變大，河道遷移擺動能力相對較強，河道呈大規模沖刷充填，以大規模側向增生為主，形成連片狀河道砂體。單井測井相多為指型與箱型，或者是指型與鐘型組合。基于地質模式劃分，結合等時地層特征，根據相對時間差異，將側疊型進一步劃分為多邊分叉式、多邊合并式以及單邊式三種模式(圖1～圖3)。對三種模式進行地震正演及實際地震響應分析表明，由于各期河道疊置連片，整體上地震響應呈連續的中強反射特征，在河道砂體側疊位置，由于存在局部微小高程差，地震響應呈能量增強、拉伸或波形非對稱畸變等特征，指示相互疊合的兩期河道砂體頂部存在規模差或高程差。

圖1 側疊型多邊分叉式地震響應特征剖面

圖2 側疊型多邊合并式地震響應特征剖面

圖3 側疊型單邊式地震響應特征剖面

2.2 孤立型地震響應

隨著基準面的不斷上升，A/S比值增大，單井測井相呈指型、箱型或正韻律鐘型，反映該沉積時期水流能力較弱，物源供給少，使河道規模變小，砂體變薄，泛濫平原泥巖沉積增加，河道砂體呈孤立型分布，砂體彼此孤立，一般互不連通。實際地震響應及地震正演結果呈透鏡型斷續反射，砂體尖滅部位的同相軸明顯不連續，可見非常明顯的復波分解和單反射的消失，組成一個典型的退積式組合，能量明顯比砂體疊置時小，較難識別(圖4)。

圖4 孤立型地震響應特征剖面

3 河流相疊合砂巖儲層地震屬性組合特征

實際工區連井剖面及地震正演分析表明，河流相砂巖儲層厚度小于地震資料調諧厚度(12m)時，河流相砂巖儲層通常以疊合形式存在，不同疊合砂巖儲層的地震響應均表現為單一波形，但如振幅、波形、頻率等特征互不相同，導致地震屬性的差異。本文在優選敏感地震屬性集合的基礎上，基于河道砂巖疊合模型，系統地分析了地震屬性值與砂巖儲層在不同疊合情況下的統計相關關系，建立了多維敏感地震屬性網狀圖。

3.1 地震屬性建立敏感屬性集合

圖5為三維河道砂體疊合模型，表征了單期河道砂體厚度變化、尖滅、多期河道砂體側向遷移疊合等特征，單期河道砂體最大厚度為7m，寬度為300m。圖6為Line 224河道砂巖疊合模型剖面及其地震響應。由圖可見：在砂巖儲層厚度小于地震資料調諧厚度的情況下，單期砂巖厚度與振幅呈正相關關系；疊置時波形出現拉伸、峰谷不對稱及能量減弱等現象，難以準確分辨砂巖儲層的疊合情況，但相對的波形變化特征較明確。

圖5 三維河道砂體疊合模型 單位：m

圖6 Line 224河道砂巖疊合模型剖面(a)及其地震響應(b)震源為主頻50Hz的Ricker子波

提取的地震屬性主要包括振幅類、頻率類和層序統計類等三大類(圖7)。針對全部地震屬性，首先，在每類屬性中分析地震屬性的相關性，優選地震屬性；其次，在優選的地震屬性集合中，分析地震屬性值在河道砂體不同部位的變化，針對不同的砂體疊置部位，形成敏感地震屬性值集合。

振幅類地震屬性(圖8a)主要突出振幅異常，可有效識別巖性及巖性組合特征。這類地震屬性間相關度很高，敏感程度基本一致，只有平均振幅、總振幅與其他地震屬性相關度較低。總振幅為窗口內正、負振幅的和，平均振幅為在此基礎上除以非零振幅的樣點數，而其他地震屬性則是對振幅值進行絕對值或者平方的計算，減小了負振幅的影響，所以相關度較差。通過分析，優選均方根振幅、振幅變化率與平均振幅作為振幅類屬性中三種敏感地震屬性。

頻率類地震屬性(圖8b)在河流相砂體預測中主要反映巖性、調諧效應等引起的子波變化。由于這類地震屬性從不同的角度利用不同的算法刻畫頻率變化，各地震屬性間不具有相關性。其中弧長是綜合頻率和振幅兩種屬性的融合屬性。通過分析，優選有效帶寬、峰值頻率和弧長作為頻率類屬性中的三種敏感地震屬性。

圖7 河道砂體疊合模型地震屬性平面圖

圖8 振幅類(a)、頻率類(b)、層序類(c)地震屬性交會分析

常將層序統計類地震屬性(圖8c)與其他地震屬性組合區分進積或者退積層序、不整合以及預測砂巖儲層的厚度。層序類屬性中(圖8c)的能量半時間、波峰數、波谷數、能量半時間斜率等四種屬性間與振幅類屬性間均不具有相關性，根據對波形結構的敏感程度，優選能量半時間為敏感地震屬性。

綜上所述，最終優選均方根振幅、平均振幅、振幅變化率、有效帶寬、峰值頻率、弧長、能量半時間等7種為敏感地震屬性，建立敏感地震屬性集合。

3.2 不同疊置模式敏感地震屬性集合特征

根據7種地震屬性主要分析了泥巖區、單期河道砂體主體與邊部、側向疊置以及多期疊置部位的地震屬性變化特征(圖9)。

圖9 7種地震屬性與河道疊置對比統計分析

(1)均方根振幅。泥巖區屬性值最小，河道砂巖區邊部至中心，屬性值在相對偏低至偏高的范圍內變化，疊置區域屬性值中等。總體上均方根振幅的屬性值與砂巖厚度和地震振幅變化具有較強的相關性，對厚度變化較敏感，反映疊置河道平面分布特征。

(2)平均振幅。泥巖區屬性值最小，河道砂體區邊部至中心，屬性值相對居中，對厚度變化不敏感。多期砂體疊置區域屬性值最高。總體上屬性變化較復雜，不僅受上、下圍巖反射影響，同時與河道砂體垂向疊合位置和厚度等因素有關，在實際應用中應結合其他地震屬性綜合分析。

(3)振幅變化率。與均方根振幅等振幅類屬性變化規律一致，反映了砂體厚度變化。與其他振幅類地震屬性相比，放大了振幅變化幅度，適合于振幅變化較小的區域。

(4)有效帶寬。與泥巖區屬性值(最大)相比，砂巖區屬性值變化范圍較小，即具有一個較窄的帶寬。依據該地震屬性特點，較窄的帶寬反映了地震反射較簡單或者平滑的反射特征，是均質介質的特征。當砂巖過渡到泥巖時，該屬性值發生明顯變化，說明該屬性結合其他屬性可以較好地反映非均質性。

(5)峰值頻率。泥巖區屬性值最高，砂體孤立時均表現為隨厚度減小屬性值增大的規律，疊置時屬性值較小，同時出現頻陷的情況(零值點)。

(6)弧長。河道孤立時，屬性值隨厚度較小而減小，表現出與振幅類屬性相同的變化特征。泥巖區屬性值較小，多期疊合時屬性值最大，表明砂體多期疊合存在泥巖夾層時屬性值會減小。

(7)能量半時間。泥巖區屬性最小，河道砂巖區邊部至中心，屬性值由低到高變化，疊置區域屬性值中等，比單期砂體屬性值大。實際應用表明，屬性變化較復雜，與上、下圍巖反射、河道砂體垂向疊置位置和厚度等因素有關，在實際工作中應結合其他地震屬性綜合分析。

基于模型地震屬性分析，通過對比實際工區井點地震屬性特征，得到不同砂體疊置模式與敏感地震屬性相對關系統計表(表1)。不同砂體疊置模式的不同單一地震屬性存在相近特征，但是7種敏感地震屬性的組合特征各不相同，可以較好地區分砂巖疊置模式。根據表1，得到河道砂體疊置地震屬性特征圖(圖10)，可以較直觀地認識不同砂體疊置模式地震屬性值域組合范圍。即在小于地震資料分辨率的情況下，采用地震屬性聚類分析等技術手段預測曲流河相砂體疊置模式、識別儲層疊置、泥巖夾層分布等情況，對油田開發階段井間儲層連通情況進行細化識別，指導油田開發生產。

表1 不同砂體疊置模式與敏感地震屬性相對關系統計

圖10 河道砂體疊合地震屬性特征

4 應用效果分析

渤海M油田Nm下段油藏主要為曲流河相沉積，由于沉積不穩定，頻繁的改道和遷移使砂體縱向互相疊合、橫向變化快，不同沉積微相砂體連通性差異很大。對測井砂、泥巖速度分析表明，砂巖的速度、密度略小于泥巖，差值約為250m/s、0.2g/cm3，在“泥包砂”背景下波阻抗差異較明顯，造成Nm下段地震剖面形成相對能量較強的地震反射。分析地震響應特征(圖11)表明，由于主力油層NmⅠ-3和NmⅡ-3復合河道側向遷移、疊加形成的砂體在平面上基本疊合連片，地震響應表現為強振幅、連續性較好，并可見到一部分復合河道遷移、疊合的響應特征。在兩套強反射之間表現為能量較弱、強弱變化明顯、橫向連續性較差、波形變化復雜的特征，反映了孤立砂體橫向尖滅的地震響應。

為了研究不同砂巖儲層疊置引起的地震屬性變化，在解釋NmⅡ-3小層砂體時，在精細井震標定的基礎上，基于砂體頂、底波峰/波谷解釋方案約束，進一步解釋波形整體形態特征(圖12)。通過提取均方根振幅、振幅變化率、平均振幅、有效帶寬、弧長、峰值頻率、能量半時間等7種屬性，將模型試驗確定的4類地震預測砂巖儲層疊合模式(圖10)作為7種敏感地震屬性聚類分析的4類目標預測結果，得到 NmⅡ-3小層與聚類分析平面圖(圖13)；結合測井資料驗證，在均方根振幅反映河道砂巖儲層整體分布特征下，聚類分析結果更為精細準確地預測了砂巖儲層內部疊合分布特征。

圖11 過井A1-B2-B3-B4-B5地震剖面

圖12 NmII-3小層地震屬性提取窗口

(1)單砂體(圖13a)。結合測井認識，即NmⅡ-3砂體與NmⅡ-4砂體之間隔層較厚，NmⅡ-3砂體受NmⅡ-4砂體的地震調諧效應影響較小，構成砂體的主體部位，地震響應為強振幅、波形對稱的連續反射。

(2)二層疊置(圖13b)。井上砂體疊置統計表明，地震響應為中等強度振幅、波形對稱，反映了單層較薄砂體或者砂體疊置，由于內部夾層很薄，互相干涉使振幅減弱，主要分布在主體砂體的周圍。

(3)三層疊置(圖13c)。地震響應為中等到低振幅、波峰拉伸、能量減弱、頻率偏高。井震標定統計表明，反映了多期砂巖儲層疊合特征，泥巖夾層很薄，主要分布在A、B平臺井區主體砂體之間。

(4)多層疊置或者泥巖區(圖13d)。振幅很弱，頻率最高，主要反映了泥巖區的分布情況。

圖13 NmII-3均方根振幅(左)與聚類分析(右)平面圖

5 結束語

根據物源供給與可容納空間比值，可將研究區新近系曲流河相沉積復合砂體疊合模式劃分為堆疊型、側疊型以及孤立型。基于該劃分模式，通過井震地震響應特征分析，進而基于地震模型分析河道砂體疊置部位的敏感地震屬性特征，形成了包括均方根振幅、振幅變化率、平均振幅、有效帶寬、弧長、峰值頻率、能量半時間等7種屬性的敏感地震屬性，并將模型試驗確定的4類砂巖儲層疊置模式定義為地震屬性聚類分析的4類預測目標結果。通過渤海實際靶區油田河流相沉積 NmⅡ-3小層的應用驗證，表明該方法可較有效地預測砂巖儲層內部疊置模式的分布特征。

需要指出的是，在海上河流相油田開發階段，在解剖儲量單元內部的儲層結構時，儲層研究尺度一般小于地震資料分辨率尺度。因此，地震響應與儲層疊置模式不具有一一對應關系，但是通過研究河道砂體疊置特征以及一定尺度內(小于調諧厚度)的砂體疊置模式，仍然能夠建立較準確的地震響應與砂巖儲層之間的關系。同時，儲層預測時不再強調準確分辨單層薄砂巖的能力，而是通過波形與地震屬性集合的差異，對薄砂巖儲層的不同疊合特征進行區分識別，進而實現砂巖儲層的精細解釋。