王斌，尹路，陳永波，滕團余，李得滋

隨著準噶爾盆地瑪湖凹陷斜坡區(qū)在瑪北和瑪西地區(qū)下三疊統(tǒng)百口泉組（T1b）油氣勘探的突破，針對環(huán)瑪湖凹陷斜坡區(qū)百口泉組砂礫巖儲集層的研究成為熱點[1-2]。已鉆井的資料證實，研究區(qū)內(nèi)目的層埋深大，特別是瑪東地區(qū)，埋深在4 000～5 000 m；有利的沉積相帶為扇三角洲前緣亞相，砂層縱向疊置嚴重，“甜點”儲集層厚度較小，一般為10～15 m[3]。僅利用常規(guī)地震資料進行疊前彈性參數(shù)反演，難以區(qū)分儲集層和非儲集層。為有效識別該區(qū)地震相帶邊界和優(yōu)質(zhì)砂礫巖儲集層，自2013年開始，在研究區(qū)實施了高密度寬方位地震資料采集。利用面向?qū)挿轿坏卣鹂碧降母叨颂幚砑夹g(shù)，能獲取炮檢距和方位角信息的疊前地震道集[4]，而這種資料需要采用特殊的處理技術(shù)，即OVT域處理技術(shù)。

OVT（Offset vector tile）通常翻譯成“炮檢距向量片技術(shù)”，最早由Vermeer在研究采集其研究區(qū)的最小數(shù)據(jù)集表達時提出[5-6]。由于原始OVT域道集的質(zhì)量受不穩(wěn)定能量、不均衡密度和較低的信噪比影響，其質(zhì)量并不穩(wěn)定，因此，為了充分發(fā)揮高密度寬方位地震勘探技術(shù)在儲集層預(yù)測中的效能，滿足瑪湖凹陷百口泉組油氣勘探的需求，通過對OVT域資料評價、預(yù)處理、動態(tài)分析和優(yōu)勢道集疊加等一系列處理后，建立了一套基于OVT域道集的疊前儲集層預(yù)測方法。

1 研究區(qū)概況

研究區(qū)位于瑪湖凹陷瑪東地區(qū)，主體位于夏鹽凸起和達巴松凸起上（圖1）。二疊紀至侏羅紀早期，瑪湖凹陷一直是盆地的沉降中心之一，沉積了巨厚的陸源碎屑巖。目的層下三疊統(tǒng)百口泉組埋藏深度為4 000～5 000 m，主要為灰綠色砂礫巖、含礫砂巖，夾灰色、綠灰色泥巖和砂質(zhì)泥巖。

圖1 研究區(qū)構(gòu)造位置

研究區(qū)三維地震資料源于2014年部署的高密度寬方位地震資料，采用正交觀測系統(tǒng)，面元為25 m×25 m，覆蓋次數(shù)為264，橫縱比為0.79.三維地震勘探面積408 km2，研究區(qū)內(nèi)已鉆探井有9口，目的層為百口泉組。

為了提高研究區(qū)“甜點”儲集層預(yù)測的精度，并充分發(fā)揮OVT域資料的優(yōu)勢，在研究區(qū)的“甜點”儲集層預(yù)測中采用如圖2所示的疊前彈性參數(shù)反演流程。

圖2 研究區(qū)“甜點”儲集層預(yù)測流程

2 優(yōu)質(zhì)儲集層成因

研究區(qū)三疊系百口泉組沉積相類型為扇三角洲相，以扇三角洲前緣亞相為主，主要發(fā)育河口壩和水下分流河道等微相；扇三角洲前緣亞相較扇三角洲平原亞相遠離物源，儲集層顆粒淘洗更為充分，儲集層物性更好，是有利的儲集相帶。水下分流河道砂礫巖與河口壩砂巖結(jié)構(gòu)較好，是目前百口泉組儲集層的主要巖石類型[7]。研究區(qū)儲集層的物性與其次生溶蝕孔隙的數(shù)量密切相關(guān)，而次生溶蝕孔隙的形成則主要由成巖溶蝕作用導致。因此，沉積相和成巖相是百口泉組優(yōu)質(zhì)儲集層的2個主要控制因素[8-9]。

2.1 扇三角洲前緣亞相是優(yōu)質(zhì)儲集層形成的物質(zhì)基礎(chǔ)

瑪東地區(qū)百口泉組儲集層中，次生溶蝕孔隙類型主要為顆粒長石及巖屑內(nèi)長石溶孔。根據(jù)儲集層巖石成分統(tǒng)計，扇三角洲平原亞相中長石的平均含量較低，為5.0%，而前緣亞相儲集層中長石顆粒平均含量較高，為14.0%；扇三角洲平原亞相中雜基的平均含量為5.2%，而前緣亞相中雜基的平均含量為2.9%，因此平原亞相的雜基含量明顯高于前緣亞相。雜基含量高，易于壓實減孔，同時在中成巖階段不利于長石溶蝕作用的持續(xù)發(fā)生。

壓汞資料顯示，百口泉組儲集層平均分選系數(shù)為1.64，平均變異系數(shù)為0.14，最大孔喉半徑為0.77 μm，孔喉體積比為1.67.總體來看，儲集層屬細喉—微細喉類儲集層。但不同相帶儲集層的結(jié)構(gòu)有所不同，對應(yīng)的儲集層孔隙也存在差異。扇三角洲前緣砂礫巖儲集層結(jié)構(gòu)最好，對應(yīng)儲集層物性也最好；其次為扇三角洲前緣河口壩砂巖相；扇三角洲平原辮狀河道砂礫巖儲集層結(jié)構(gòu)最差，對應(yīng)物性也最差。水下分流河道砂礫巖具有較低的排驅(qū)壓力（0.71 MPa），平均毛細管半徑為0.26 μm，但孔喉分選性差，非均質(zhì)性強；河口壩砂巖排驅(qū)壓力相對較高（1.12 MPa），平均毛細管半徑為0.25 μm，但孔喉分選性較水下分流河道好；扇三角洲平原砂礫巖具有較高的排驅(qū)壓力（1.43 MPa），平均毛細管半徑為0.12 μm，孔喉分選性差，非均質(zhì)性強。從儲集層結(jié)構(gòu)看，前緣亞相的儲集層結(jié)構(gòu)明顯好于平原亞相，有利于溶蝕作用的發(fā)生。

2.2 溶蝕作用是優(yōu)質(zhì)儲集層的主要發(fā)育機制

在前緣亞相儲集層中，長石的溶解是優(yōu)質(zhì)儲集層發(fā)育的主要控制作用。油氣的侵位與黏土成巖演化也能導致酸性成巖環(huán)境的形成?，敄|地區(qū)百口泉組存在2期油氣充注，分別是風城組在早侏羅世的成熟油和早白堊世的高熟油充注，2次油氣侵位為長石溶蝕提供了酸性成巖環(huán)境。此外，酸性成巖環(huán)境中長石發(fā)生溶解需要滿足80～140℃的地溫條件，該溫度條件對應(yīng)的成巖階段主要為中成巖階段，這與中成巖階段長石的溶蝕具有很好的對應(yīng)關(guān)系。因此，瑪東斜坡區(qū)百口泉組長石溶蝕主要發(fā)生在早白堊世的晚期高熟油侵位階段，這就決定了“甜點”儲集層主要發(fā)育在酸性流體的優(yōu)勢運移通道上。

3 道集分析與處理

3.1 原始OVT域道集分析

原始OVT域道集的質(zhì)量是不穩(wěn)定的，如不穩(wěn)定的能量、不均衡的密度和較低的信噪比。對OVT域道集進行部分疊加，可提高道集質(zhì)量及信噪比[10-11]。

通過OVT域道集動態(tài)分析技術(shù)，可選取適當?shù)牡兰B加參數(shù)，以改善成像質(zhì)量。此外，OVT域道集動態(tài)分析和疊加，可在滿足地質(zhì)目標分析精度的前提下，減小OVT域道集的數(shù)據(jù)量，提高OVT域道集地震屬性分析的效率。

道集的動態(tài)分析，主要利用了基于各向異性介質(zhì)理論的地震道集疊加成像技術(shù)。各向異性介質(zhì)的縱波反射系數(shù)的近似公式，是入射角和測線方位角的函數(shù)[12-14]：

式中 C11，C13，C33，C44，C55——各向異性參數(shù)，無量綱；

γ=(C55-C44)/2C44；

θ——入射角，（°）；

φS——裂縫走向的方位角，（°）；

φ——測線方位角，（°）；

α——垂直入射縱波速度，m/s；

β——垂直入射橫波速度，m/s；

ρ——密度，g/cm3.

當入射角較小時，（1）式可簡化為

式中 A=(1/2)(ΔZ/Z)；

因此，在直角坐標系中波形起伏不大的余弦曲線，在極坐標系中會近似為橢圓。

利用基于模板的地質(zhì)目標驅(qū)動的地震道集疊加成像技術(shù)，即根據(jù)儲集層特征，選擇不同的疊前道集疊加模板，增強“甜點”儲集層的主要特征，提高疊前地震資料的信噪比。所謂疊加模板是一種炮檢距（入射角）-方位角動態(tài)選擇工具，同時也是一種OVT域道集的炮檢距與方位角聯(lián)合分析工具。分區(qū)疊加成像是指根據(jù)地質(zhì)條件的不同將一個工區(qū)劃分為兩個或多個區(qū)域，每個區(qū)域采用不同的疊加模板。圖3就是利用動態(tài)分析技術(shù)在研究區(qū)能夠識別“甜點”儲集層特征的優(yōu)勢道集，圖中色標表示優(yōu)勢道集的方位角角度，即在該區(qū)域的優(yōu)勢方位角區(qū)間。從圖3中可知，在研究區(qū)的東部和西部，優(yōu)勢方位角大于120°；在北部和南部，優(yōu)勢方位角為60°～120°；而中部的優(yōu)勢方位角小于60°，因此將研究區(qū)分為5個區(qū)域，在不同區(qū)域采用不同的疊加參數(shù)，對地震道集做疊加處理。

圖3 研究區(qū)百口泉組OVT域資料優(yōu)勢方位角平面分布

3.2 OVT域疊前道集預(yù)處理

首先分析疊前共反射點（CRP）道集的質(zhì)量，根據(jù)資料情況進行CRP道集的優(yōu)化處理，基于巖性預(yù)測的CRP道集優(yōu)化處理內(nèi)容包括5個方面：與炮檢距有關(guān)的吸收補償、動校拉伸校正、紡錘狀反射校正、道集不平校正以及去噪和提頻處理。在基于巖性預(yù)測的CRP道集優(yōu)化處理中，重點是與炮檢距有關(guān)的吸收補償（針對保真度）和大角度道集的剩余時差校正（針對寬角度）處理，而基礎(chǔ)內(nèi)容是提高信噪比和分辨率的處理。典型的基于巖性預(yù)測的CRP道集優(yōu)化處理流程是“去噪—吸收補償—提頻—拉平校正”[15]。

與正演道集和原始道集對比，優(yōu)化處理后的道集質(zhì)量明顯提升，信噪比得到明顯提高，能量關(guān)系與正演道集有較好的對比性（圖4）。

圖4 過D13井正演道集(a)、原始道集(b)和處理后道集(c)對比

4 應(yīng)用效果分析

4.1 “甜點”儲集層敏感參數(shù)優(yōu)選

根據(jù)優(yōu)質(zhì)儲集層成因和研究區(qū)鉆井儲集層物性及儲集層結(jié)構(gòu)統(tǒng)計分析結(jié)果，對儲集層進行了分類。儲集層主要為扇三角洲前緣水下分流河道砂礫巖及河口壩砂巖，“甜點”孔隙度大于10%，一般儲集層孔隙度為8%～10%.非儲集層主要為湖相泥巖和扇三角洲平原亞相及前緣亞相的部分致密砂礫巖，孔隙度小于8%.

基于研究區(qū)儲集層的特殊性，通過巖石物理測試，明確了低滲透砂礫巖的巖石物理特征。研究區(qū)低滲透砂礫巖受到孔隙超壓和塑性碎屑含量的影響，與常規(guī)砂巖相比，具有低縱橫波速度、高縱橫波速度比、各向異性強和高動靜態(tài)彈性模量比的特點，將彈性參數(shù)與孔隙度做交會圖，經(jīng)分析和優(yōu)選，縱橫波速度比與儲集層“甜點”相關(guān)性較好。因此，選定縱橫波速度比為研究區(qū)“甜點”儲集層的敏感屬性?？v橫波速度比與孔隙度具有反相關(guān)關(guān)系，即孔隙度高的樣品具有較低的縱橫波速度比?！疤瘘c”儲集層的孔隙度大于10%，縱橫波速度比小于1.78（圖5）。

圖5 研究區(qū)縱橫波速度比與孔隙度交會圖

在此基礎(chǔ)上，通過研究區(qū)內(nèi)測井曲線的交會，進一步證實了縱橫波速度比能較好地識別儲集層“甜點”。圖6為縱橫波速度比與縱波波阻抗的巖石物理解釋模板，散點為已鉆井提供的資料點。用鉆井資料驗證模板準確與否，兩者相近說明該巖石物理模板準確可行。

圖6 研究區(qū)疊前敏感參數(shù)優(yōu)選

4.2 基于OVT域資料的疊前反演結(jié)果分析

從反演得到不同入射角下的彈性阻抗體出發(fā)，可計算出縱、橫波速度和密度等基本巖性參數(shù)。最后根據(jù)各巖性參數(shù)間的相互關(guān)系，可計算得到拉梅常數(shù)、剪切模量、泊松比等能直接反映巖性特征的參數(shù)[16]。

圖7為研究區(qū)疊前縱橫波速度比反演平面圖，圖中紅黃色為應(yīng)用圖6中巖石物理定量解釋模板解釋的“甜點”儲集層發(fā)育區(qū)。預(yù)測研究區(qū)內(nèi)“甜點”儲集層發(fā)育區(qū)面積近140 km2，從圖7可以看出，研究區(qū)西北部有望成為下一步井位部署的接替區(qū)域。針對百口泉組整體部署，多口井分步試驗。其中D15井獲得高產(chǎn)油氣流，日產(chǎn)油32.3 m3.

圖7 研究區(qū)疊前彈性參數(shù)反演平面分布

4.3 應(yīng)用效果對比

分別利用常規(guī)三維地震資料與基于OVT域資料進行疊前彈性參數(shù)反演，基于OVT域資料的疊前彈性參數(shù)反演的分辨率較高，能更好地識別低滲透砂礫巖中的儲集層“甜點”（圖8）。

圖8 不同資料疊前彈性參數(shù)反演的效果對比

5 結(jié)論

（1）OVT技術(shù)的出現(xiàn)為高精度疊前地震解釋提供了有效途徑，其豐富的炮檢距和方位角信息，提高了復雜巖性油氣藏的“甜點”儲集層預(yù)測精度。

（2）沉積相和成巖相是瑪東地區(qū)百口泉組優(yōu)質(zhì)儲集層的2個主要控制因素，扇三角洲前緣亞相是優(yōu)質(zhì)儲集層形成的物質(zhì)基礎(chǔ)，成巖溶蝕作用是優(yōu)質(zhì)儲集層的主要發(fā)育機制。

（3）利用動態(tài)分析和部分疊加處理的OVT資料，提高了疊前儲集層預(yù)測的精度，該技術(shù)在類似地區(qū)有很好的推廣應(yīng)用前景。

［1］ 瞿建華，張磊，吳俊，等.瑪湖凹陷西斜坡百口泉組砂礫巖儲集層特征及物性控制因素［J］.新疆石油地質(zhì)，2017，38（1）：1-6.

QU Jianhua，ZHANG Lei，WU Jun，et al.Characteristics of sandy conglomerate reservoirs and controlling factors on physical proper?ties of Baikouquan formation in the western slope of Mahu sag，Jung?gar basin［J］.Xinjiang Petroleum Geology，2017，38（1）：1-6.

［2］ 雷德文，阿布力米提，唐勇，等.準噶爾盆地瑪湖凹陷百口泉組油氣高產(chǎn)區(qū)控制因素與分布預(yù)測［J］.新疆石油地質(zhì)，2014，35（5）：495-499.

LEI Dewen，ABULIMITI，TANG Yong，et al.Controlling factors and occurrence prediction of high oil?gas production zones in Lower Tri?assic Baikouquan formation of Mahu sag in Junggar basin［J］.Xinji?ang Petroleum Geology，2014，35（5）：495-499.

［3］ 婁兵，姚茂敏，羅勇，等.高密度寬方位地震數(shù)據(jù)處理技術(shù)在瑪湖凹陷的應(yīng)用［J］.新疆石油地質(zhì)，2015，36（2）：208-213.

LOU Bing，YAO Maomin，LUO Yong，et al.Application of high den?sity and wide azimuth seismic data processing technology in Mahu sag，Junggar basin［J］.Xinjiang Petroleum Geology，2015，36（2）：208-213.

［4］ 劉依謀，印興耀，張三元，等.寬方位地震勘探技術(shù)新進展［J］.石油地球物理勘探，2014，49（3）：596-610.

LIU Yimou，YIN Xingyao，ZHANG Sanyuan，et al.Recent advances in wide?azimuth seismic exploration［J］.Oil Geophysical Prospect?ing，2014，49（3）：596-610.

［5］ 詹仕凡，陳茂山，李磊，等.OVT域?qū)挿轿化B前地震屬性分析方法［J］.石油地球物理勘探，2015，50（5）：956-966.

ZHAN Shifan，CHEN Maoshan，LI Lei，et al.OVT?domain wide?azi?muth prestack seismic attribute analysis［J］.Oil Geophysical Pros?pecting，2015，50（5）：956-966.

［6］ VILLIAMS M，JENNER E.Interpreting seismic data in the presence of azimuthal anisotropy or anisotropy in the presence of the seismic interpretation［J］.The Leading Edge，2002，21（8）：771-774.

［7］ 于興河，瞿建華，譚程鵬，等.瑪湖凹陷百口泉組扇三角洲礫巖巖相及成因模式［J］.新疆石油地質(zhì)，2014，35（6）：619-627.

YU Xinghe，QU Jianhua，TAN Chengpeng，et al.Conglomerate litho?facies and origin models of fan deltas of Baikouquan formation in Mahu sag，Junggar basin［J］.Xinjiang Petroleum Geology，2014，35（6）：619-627.

［8］ 楊曉萍，趙文智，鄒才能，等.川中氣田與蘇里格氣田“甜點”儲層對比研究［J］.天然氣工業(yè)，2007，27（1）：4-7.

YANG Xiaoping，ZHAO Wenzhi，ZOU Caineng，et al.Comparison of formation conditions of“sweet point”reservoirs in Sulige gasfield and Xiangxi group gasfield in the central Sichuan basin［J］.Natural Gas Industry，2007，27（1）：4-7.

［9］ 廖建波，劉華清，林衛(wèi)東.鄂爾多斯盆地山城-演武地區(qū)三疊系延長組長6-長8低滲儲層特征及成巖作用研究［J］.天然氣地球科學，2006，17（5）：682-687.

LIAO Jianbo，LIU Huaqing，LIN Weidong.Low permeable sand?stone reservoir characteristics and diagenesis of Chang 6，Chang 7 and Chang 8 formation of the Yanchang formation（Triassic），Shancheng?Yanwu area，Ordos basin［J］.Natural Gas Geoscience，2006，17（5）：682-687.

［10］ 段文勝，李飛，李世吉，等.OVT域疊前偏移衰減多次波［J］.石油地球物理勘探，2013，48（增刊1）：36-41.

DUAN Wensheng，LI Fei，LI Shiji，et al.Multiple attenuation by prestack migration in the OVT domain［J］.Oil Geophysical Pros?pecting，2013，48（Supp.1）：36-41.

［11］ 袁燎，張麗娟，溫鐵民.OVT域處理技術(shù)在準噶爾盆地東部致密油勘探中的應(yīng)用及效果［J］.長江大學學報（自然版），2015，12（35）：22-26.

YUAN Liao，ZHANG Lijuan，WEN Tiemin.Application of offset vector tile（OVT）processing technique in tight oil exploration in the east of Junggar area［J］.Journal of Yangtze University（Natural Science Edition），2015，12（35）：22-26.

［12］ RUGER A.Reflection coefficients and azimuthal AVO analysis in anisotropic media［D］.Golden Colorada，USA：Colorado School of Mines，1996.

［13］ RUGER A.Variation of P?wave reflectivity with offset and azimuth in anisotropic media［J］.Geophysics，1998，63（3）：935-947.

［14］ 陳天勝，魏修成，劉洋.一種新的各向異性彈性阻抗近似公式［J］.石油物探，2006，45（6）：563-569.

CHEN Tiansheng，WEI Xiucheng，LIU Yang.New approximation formula for calculation of elastic impedance in anisotropic media［J］.Geophysical Prospecting for Petroleum，2006，45（6）：563-569.

［15］ 許多年，尹路，瞿建華，等.低滲透砂礫巖“甜點”儲層預(yù)測方法及應(yīng)用——以準噶爾盆地瑪湖凹陷北斜坡區(qū)三疊系百口泉組為例［J］.天然氣地球科學，2015，26（增刊1）：154-161.

XU Duonian，YIN Lu，QU Jianhua，et al.Prediction method of the low permeability sandy?conglomerate“sweet point”reservoirs and it’s application：a case study of Mahu depression northern slope ar?ea in the Junggar basin［J］.Natural Gas Geoscience，2015，26（Supp.1）：154-161.

［16］ 古發(fā)明，李進步，鄒新寧，等.炮檢距向量片技術(shù)在蘇里格致密砂巖儲層預(yù)測中的應(yīng)用［J］.地球物理學進展，2017，32（2）：610-617.

GU Faming，LI Jinbu，ZOU Xinning，et al.Application of the offset vector tile technique in the Sulige tight sandstone reservoir predic?tion［J］.Progress in Geophysics，2017，32（2）：610-617.