基于疊前OVT域偏移的河道砂體預測方法

周華建

（中國石油大慶油田有限責任公司勘探開發研究院，黑龍江大慶163712）

0 引言

大慶長垣薩爾圖油田為大型河流—三角洲體系的儲層沉積特征［1］，縱向上油層多，砂泥互層頻繁，平面上相變快，剩余油分布呈現總體高度分散、局部相對富集的特征。利用常規疊后地震儲層預測方法只能預測約1/4相對較厚的河道砂體。目前疊前地震OVT域處理解釋技術（炮檢距向量片或炮檢距矢量片）在勘探階段應用較多，在識別巖性及流體中展現了較大的潛力［2-4］，但在開發階段，針對砂泥巖薄互層儲層預測方面處于探索階段，該技術的最大優勢是能以地層的方位各向異性分析為手段預測儲層［5-6］。

從疊前地震OVT域處理工作入手，以地震沉積學為手段，充分挖掘疊前地震資料中河道砂體的地震反射信息，分析不同方位地層切片反映儲層砂體特征的振幅屬性，結合測井數據進行儲層描述，可進一步提高儲層河道砂體展布特征預測的可靠性，以期更好地指導井組間剩余油描述和挖潛。

1 研究區地震資料品質分析

大慶長垣薩爾圖油田采用的是16線6炮正交線束狀觀測系統，采用3 520（220道×16線）道接收、10 m（縱）×10 m（橫）面元、總覆蓋次數為80次（縱向10次，橫向8次）覆蓋。其地震資料主頻約為45 Hz，頻帶寬度為5～150 Hz，從合成記錄與剖面的對比分析，可以得出井震匹配程度很高，波組特征明顯，地層接觸關系清楚，層間信息豐富，信噪比和分辨率均較高，全區能量相對均衡。研究區的地震資料品質能夠滿足疊前OVT域處理和解釋方法的研究要求。

2 地震疊前OVT域儲層預測方法

地震疊前OVT域儲層預測方法采用疊前地震處理和解釋一體化工作模式，通過疊前高分辨率處理、基于炮檢距向量片成像及疊前地震儲層預測方法研究，提高儲層描述精度。

2.1 地震疊前OVT域處理

在常規地震資料處理的基礎上，通過大量的試驗，初步確定了地震疊前OVT域處理的基本步驟、關鍵環節及其質量控制方法，實現了OVT域規則化處理和偏移成像成果。流程中最主要的2個環節是各向異性分析與校正和OVT域數據插值及偏移成像（圖 1）。

圖1 OVT域處理技術流程Fig.1 OVT domain processing technology

2.1.1 各向異性地震響應特征分析與校正

隨著地震技術的發展，大炮檢距和寬方位角采集技術開始被廣泛地采用，在得到更多地質信息的同時，也使得各向異性現象凸顯，以地球介質各向同性的假設為基礎的研究已經不能滿足需要。地震資料處理中，地震各向異性通常表現為速度上的各向異性，即在同一個共成像點道集中，來自不同方向地震波射線的地層速度存在差異，造成動校正后的同相軸難以被真正拉平［7-10］。此時速度不再是標量，它的大小不僅與空間位置相關，還與傳播方向有關。因此，通常在同一點使用相同速度的動校正技術，難以使CMP道集同相軸得到拉平。

與動校正過程中速度過高或過低導致的同相軸隨炮檢距變化而下彎或上翹現象不同，在按炮檢距大小排列的CMP道集中，同相軸呈現參差不齊、隨機跳動現象的特征［圖2（a）］，而且這種情況很難判斷是哪些因素引起的，因為噪音干擾、靜校正不準確以及方位速度差異都有可能產生類似的現象。將動校正后CMP道集按方位角大小進行排列，可知，同相軸隨著方位角起伏變化［圖2（b）］，起伏時差可高達5 ms，其起伏變化具有一定的周期性，當方位角相差180°時，相當于檢波點和激發點位置進行了互換，地震射線路徑未發生改變，反射時間基本一致，即其變化周期大約為180°。

圖2CMP道集顯示Fig.2 CMPgather display

速度方位各向異性校正處理方法可分為3個步驟：①采用基于模型的相關時差計算方法，獲取非正常反射時差；②根據各向異性分析的最大時差大小，剔除異常值；③將得到的時差作用于地震數據，實現各向異性校正。得到的結果是同相軸橫向連續性增強，規律性更明顯（圖3）。

圖3 方位各向異性校正前后對比Fig.3 Comparison before and after azimuthal anisotropy correction

通過時差校正，不同OVT單元之間保持一致，層位追蹤和斷層解釋更加精準，能提高斷裂解釋和儲層預測的精度。區域內分布不同的斷裂系統，就會有不同的走向和傾向，同一個斷裂系統不同的斷層也會具有不同的走向和傾向。常規處理的地震數據只能接收垂直于采集方向的斷裂系統的信息，而疊前OVT域數據能夠從不同方位充分接收到各種走向斷裂系統的地震信息，從而能使各種走向的斷裂系統都得到較好的成像結果。

圖4OVT道集劃分Fig.4 OVT gather division

圖5 十字排列中的OVT劃分及成像圖Fig.5 OVT division and imaging map in cross arrangement

2.1.2 OVT數據插值及偏移成像

OVT道集相當于含有相對固定方位角的炮檢距道集，稱為炮檢距向量片。OVT數據集是十字排列道集的自然延伸，是十字排列道集內的1個數據子集［11-12］。圖 4（a）為十字排列與 OVT 關系示意圖，描述了十字排列中OVT劃分方法，圖中箭頭方向指向OVT道集中心方位，在十字排列中按炮線距和檢波線距等距離劃分得到許多小矩形，1個矩形就是1個OVT炮檢距向量片。

十字排列可由正交觀測系統抽取出來［圖4（b）］，即把來自同一炮線和同一檢波線的所有地震道集合起來，因此十字排列的個數與炮線和檢波線交點的數目一樣多，每個OVT都是沿炮線有限范圍內的炮點和沿檢波線有限范圍內的檢波點構成，這2個范圍把OVT的取值限制在1個小區域內，OVT可以限制炮檢距和方位角的范圍。

相對于常規處理地震道集，OVT道集兼具炮檢距和方位角2類信息，在1個OVT道集中，各地震道的炮檢距和方位角大致相同，無論是在近炮檢距、中炮檢距還是遠炮檢距范圍，OVT道集的能量一致性均較好。根據大慶長垣薩爾圖油田研究區地震采集觀測系統設計的參數，包括炮線距、檢波線間距以及80次覆蓋次數等野外施工參數，在十字交叉排列道集基礎上，將地震記錄劃分為80個有效OVT單元（圖5），每個單元對應1個完整的相當于1個單次覆蓋的地震數據體，各自含有較小范圍的炮檢距和方位角信息，為后續疊前解釋奠定了基礎。

2.2 常規地震與OVT域處理資料對比分析

常規地震處理技術采用多次疊加的方法，多次疊加的平均效應消除了地震數據中隨不同的方位角和炮檢距變化的儲層各項異性信息，無法把這些信息用于地震解釋和儲層預測，而疊前OVT域處理技術采用單炮次或部分疊加方法，充分保留了反映儲層方位各項異性的方位角和炮間距信息。如圖6所示，與常規地震處理剖面及頻譜對比，OVT域地震剖面上同相軸數量增多（圖中藍色圈所示），地震資料頻帶更寬，主頻更高，地震分辨率更高，反映儲層不同方位的有效信息加強［13-14］。

從成像效果對比分析來看（圖7），常規地震屬性信噪比高，地層切片反映儲層較清晰，但與OVT域地層切片相比，有些儲層信息沒有預測出來（圖中綠圈所示）。OVT域地震資料與常規地震資料相比，優勢在于能夠從不同方位反映儲層有效信息，最大限度地保留地震信息的有效性，提高儲層預測的精度。

圖6 常規與OVT域地震處理剖面及頻譜對比Fig.6 Comparison of seismic processing section and spectrum by conventional seismic and prestack seismic in OVT domain

圖7 OVT域與常規地震成像效果對比Fig.7 Comparison of effect between OVT domain and conventional seismic imaging

2.3 基于地震沉積學的儲層預測

得到多個不同方位、炮間距的地震數據體后，采用地震沉積學中的地震屬性切片技術進行儲層預測分析，其主要工作流程是：①基于地震構造解釋結果建立小層級等時地層格架，實現小層級地層切片的制作；②通過地層切片與開發地質解釋成果數據計算相關性，實現響應小層地層切片預測結果優選，并結合研究區沉積環境特征及測井解釋結果，分析識別地震屬性切片預測的地質體［15-19］。

2.3.1 等時地層格架建立

利用地層切片技術預測儲層的基礎是建立等時地層框架［20-22］。研究區構造平坦，目的層段無缺失或地層尖滅，地震資料顯示油層組級地層同相軸穩定、連續，油層組與小層級地勢形態基本一致，井震標定相關性達85%以上。通過地震合成記錄標定和地震同相軸追蹤，建立相對較厚時間域的油層組級等時地層格架，然后以油層組SⅡ，SⅢ等時地層格架為約束，在油層組內部按照井點分層的等比例原則建立沉積單元級地層格架（圖8）。所得分層與實際井點分層深度誤差在±1 m內的井達90%以上，精度得到了保證，為后續相應沉積單元的儲層預測精度提供了保障。

2.3.2 地層切片預測結果優選

研究區工區面積20 km2，目的層為分流平原亞相沉積，復合河道砂體和點壩砂體發育［1］。利用上述方法建立全區地層格架，在多個含有不同方位、炮檢距的地震數據體中的一定時窗內提取相應層位的地震樣點值，實現多個地層切片的制作［23-26］。不同方位角體和常規地震體預測同一小層（基于振幅屬性），預測儲層整體趨勢基本一致，但方位角體整體連續性更好，砂體邊界更清晰；不同方位角體預測結果存在細節上的差異。如圖9所示，選取4個區域應用地貌學特征進行詳細對比分析，其中區域1為復合河道砂體發育帶，OVT域地層切片與常規地層切片［圖9（a）］相比，識別出了單一河道砂體［圖9（b）中綠色箭頭所指］；區域2為窄小河道砂體發育帶，OVT域地層切片預測的河道砂體更加連續［圖9（c）中綠色箭頭所指］；區域3和4為點壩砂體發育帶，OVT域地層切片預測的點壩砂體邊界更清晰［圖9（c），（d）中綠色弧線所示］。

圖8 沉積單元級等時地層格架Fig.8 Isochronic stratigraphic framework at sedimentary unit level

利用定量方法對研究區A小層不同方位角體地層切片與常規地層切片預測儲層砂體精度進行對比分析。采用如下Pearson線性相關理論公式，計算井點砂巖厚度值與井旁道地震振幅值的相關性［27-28］。測井解釋儲層砂體發育區與地震振幅屬性預測砂巖發育區域匹配程度越高（圖10），所得相關系數R越大，證明地震預測儲層精度越高。

圖10 井點砂巖厚度與地震振幅關系Fig.10 Relationship between well-point sandstone thickness and seismic amplitude

式中：R表示相關系數；x為井點地震振幅值；y為井點砂巖厚度，m；N為井數，口。

以圖9中對應的4個區域為計算目標區，選取區域中大于2 m的井點砂巖厚度值與井旁道地震振幅值做相關性計算（表1），可知，OVT域地層切片預測的儲層河道砂體精度高于常規地層切片預測的結果。

表1 砂巖厚度值與振幅值相關性計算結果Table 1 Relative calculation results of sandstone thickness and amplitude

3 井震結合河道砂體描述實例

以大慶長垣薩爾圖油田A小層局部發育單一河道砂體區域為例，如圖11所示，圖中儲層預測結果顏色由藍→白→紅變化表示巖性由泥巖向砂巖逐漸過渡，OVT域地層切片［圖11（a）］預測砂巖厚度2 m以上的井點與地震預測結果符合率超過70%，比常規地層切片［圖11（b）］預測的結果符合率提高了15%，而且預測出枝狀的河道砂體，反映的砂體更加連續清晰。

圖11 A小層局部地層切片及井點砂巖解釋厚度Fig.11 Local stratal slices and interpretation thickness of well-point sandstone inAsublayer

依據圖11（a）預測的儲層砂體展布趨勢，結合預測條帶上及周邊井點相應層位砂體的解釋厚度、測井曲線形態等對比確認［29-32］，刻畫出沉積微相圖［圖 12（a）］，按照同樣的方法根據圖 11（b）刻畫出圖 12（b）。以注入井H井為例，按圖12（b）的刻畫結果，該井分布在席狀砂體中，統計該區A小層120口分布在席狀砂體內的水井日均吸水量7 m3，而實際上H井在A小層的日均吸水量為15 m3，兩者存在明顯的差異，開發動態數據統計結果與對該井的地質認識存在矛盾，而按圖12（a）的刻畫結果，H井分布在河道砂體中，統計研究區A小層96口分布在河道砂體內的水井日均吸水量為14 m3，和H井的實際吸水量非常接近，說明H井在A小層應分布在河道砂體內。由此可知，以基于疊前OVT域地層切片預測的儲層砂體展布趨勢為依據的相圖精度更高。

圖12 A小層局部井震結合的相Fig.12 Phase of local well-seismic combination inAsublayer

4 結論

（1）疊前地震OVT域處理解釋技術，與常規地震資料相比，能從不同方位角體的反射記錄中反映出地層的方位各向異性特征，能夠提高儲層地質體預測的準確性，是長垣薩爾圖油田密井網區儲層預測的一種有效方法。

（2）依據OVT域疊前地震屬性儲層預測成果，井震結合對儲層河道砂體進行精細刻畫，增強了儲層河道砂體識別的確定性，對開發調整有很好的指導意義。