劉少治，王 波，宋 兵，陳 昌，陳希光

（1.中國石油杭州地質(zhì)研究院，浙江 杭州 310023；2.中國石油遼河油田分公司勘探開發(fā)研究院，遼寧 盤錦 124010）

遼河西部凹陷北部曙北地區(qū)沙四段地層厚度大（平均226.7 m，最大773.0 m）、砂巖厚度薄（單層厚度一般2.0～6.0 m，累計厚度平均30.2 m）、橫向變化快。隨著增儲挖潛工作逐步向主儲量區(qū)外轉(zhuǎn)移，曙134井的高產(chǎn)薄砂巖（砂巖厚4.8 m／2層，日產(chǎn)油9.7 m3）日益受到重視，是下步巖性油氣藏或致密油氣勘探的潛在領(lǐng)域。但是由于區(qū)內(nèi)沙四段地震資料品質(zhì)差（主頻18 Hz，頻帶寬度5～30 Hz）、巖性復(fù)雜（既有泥巖、油頁巖、碳酸鹽巖等細(xì)粒沉積巖，又有下伏的玄武巖）、現(xiàn)今斷裂發(fā)育等因素影響，加之傳統(tǒng)的波阻抗反演、地質(zhì)統(tǒng)計學(xué)反演、隨機(jī)反演難以準(zhǔn)確預(yù)測薄砂巖儲層，制約了區(qū)內(nèi)油氣勘探。目前曲線重構(gòu)擬聲波反演［1-2］、波形指示反演［3-6］結(jié)合相控反演［7-9］在薄互層砂巖儲層預(yù)測中取得了較好的應(yīng)用效果，但是對復(fù)雜巖性區(qū)薄砂巖儲層預(yù)測仍需深入探索。本文在高頻等時層序格架的基礎(chǔ)上，充分吸收上述反演方法的優(yōu)點(diǎn)，對曙北地區(qū)沙四段復(fù)雜巖性區(qū)薄砂巖進(jìn)行儲層預(yù)測。

1 地質(zhì)概況

遼河西部凹陷是渤海灣盆地主要的富油氣凹陷之一，研究區(qū)位于凹陷中北部，曙光以北、高升以南、雷家以西地區(qū)，面積約為480 km2。截至目前已經(jīng)發(fā)現(xiàn)曙光潛山油氣藏，曙光、高升沙四段構(gòu)造油氣藏和雷家沙四段碳酸鹽巖油氣藏。西部凹陷新生代的演化經(jīng)歷了地殼拱張、裂陷和坳陷三個階段，其中裂陷階段又進(jìn)一步分為初陷期、深陷期、持續(xù)衰減期三個發(fā)育期，沙四段是遼河裂谷盆地初陷期的產(chǎn)物。曙北-高升地區(qū)沙四段巖性以灰色、深灰色泥巖和油頁巖為主，夾薄層細(xì)砂巖、泥質(zhì)粉砂巖，整體表現(xiàn)為厚層泥巖（厚度200.0～300.0 m）夾少量薄層砂巖（厚度2.0～6.0 m）；雷家地區(qū)沙四段巖性則為厚層泥巖夾薄層碳酸鹽巖，碳酸鹽巖主要為砂質(zhì)、泥質(zhì)、灰質(zhì)白云巖及泥灰?guī)r。對沙四段沉積分析認(rèn)為，曙北-高升地區(qū)為扇三角洲-湖泊過渡沉積體系［10-11］，雷家地區(qū)為扇三角洲-（鮞）粒屑灰?guī)r體系［12-13］，薄砂巖儲集體主要為遠(yuǎn)源水下分流砂體和濱淺湖灘壩砂體。

2 等時相控波形指示模擬反演

2.1 方法

在高精度層序格架基礎(chǔ)上，利用波形指示模擬反演軟件，在相控建模的基礎(chǔ)上，將研究區(qū)井旁道波形特征作為樣本，采用SMCMC算法，擬合有效樣本數(shù)、設(shè)定最佳截止頻率和濾波器參數(shù)，以達(dá)到最佳的反演效果。

2.2 流程

2.2.1 可行性分析

針對薄砂巖儲層預(yù)測前人做了大量的正演模型分析工作［2-4，9，14-18］，儲層類型包括河道砂巖、濱湖砂巖以及致密砂巖等。灘壩砂體的正演模型分析表明，利用波形特征可以反映灘壩砂巖儲層的發(fā)育特征，不同巖性組合的灘壩砂巖其敏感的地震屬性不同［14］；砂地比、砂泥巖互層頻率、是否含灰、是灘砂還是壩砂都會影響地震屬性，因此地震屬性預(yù)測應(yīng)當(dāng)分區(qū)分帶進(jìn)行［3］。研究區(qū)斷裂復(fù)雜、目的層埋深跨度大（500.0～4 500.0 m）、地震資料主頻低（18～20 Hz），沙四段巖性復(fù)雜、儲層類型較多、砂巖厚度薄、橫向變化快，地震縱向分辨率導(dǎo)致難以識別、追蹤并刻畫薄砂巖儲層。因此首先應(yīng)當(dāng)發(fā)揮不同反演軟件的各自優(yōu)勢進(jìn)行儲層預(yù)測，尤其是波形指示模擬反演；其次必須充分利用鉆井、地震屬性及屬性融合對反演結(jié)果進(jìn)行約束，以提高預(yù)測可信度。

2.2.2 擬聲波曲線重構(gòu)

擬聲波曲線重構(gòu)的目的是通過將聲波曲線的低頻信息和對巖性變化比較敏感的非聲波測井曲線的高頻信息融合，從而能夠較好地區(qū)分巖性變化［1，19］。在對多口井的波阻抗、自然電位、自然伽馬和電阻率等測井曲線分別進(jìn)行泥巖、油頁巖和砂巖統(tǒng)計（圖1）的基礎(chǔ)上，發(fā)現(xiàn)波阻抗和自然伽馬對砂巖、泥巖都有一定的區(qū)分能力，但是波阻抗不能很好地區(qū)分砂巖和油頁巖，自然電位曲線對厚層砂巖識別較好。因此通過融合波阻抗低頻信息和自然伽馬的高頻信息形成GR擬聲波曲線，可以很好地區(qū)分薄砂巖與油頁巖（圖2），從而提高后期儲層反演的分辨率和精度。在可行性分析的基礎(chǔ)上，對聲波和自然伽馬曲線分別進(jìn)行井環(huán)境校正、巖電校正和標(biāo)準(zhǔn)化處理，隨后利用等時相控波形指示模擬反演軟件特征曲線重構(gòu)功能進(jìn)行GR擬聲波曲線重構(gòu)。

圖1 不同測井曲線對砂泥巖及油頁巖概率統(tǒng)計圖

圖2 曙134井GR擬聲波曲線可行性分析

2.2.3 等時相控建模

在前人對西部凹陷三級層序研究的基礎(chǔ)上，通過測井曲線降頻法結(jié)合區(qū)域?qū)Ρ葮?biāo)志層（房身泡組玄武巖、高升油層頂部低平泥巖段和杜家臺油層頂部特殊巖性段）將沙四段三級層序自下而上劃分為高升油層和杜家臺油層兩個準(zhǔn)層序組，并進(jìn)一步細(xì)分為高Ⅰ砂組、高Ⅱ砂組、杜Ⅰ砂組、杜Ⅱ砂組、杜Ⅲ砂組5個準(zhǔn)層序（圖3），各砂組地層厚度平均為40.0～60.0 m。三級層序底界面為一個區(qū)域性角度不整合面，在地震剖面上表現(xiàn)為一個強(qiáng)振幅、連續(xù)性好的地震反射同相軸，可以進(jìn)行全區(qū)追蹤。三級層序的頂界面在地震剖面上不同地區(qū)具有差異，在鉆測井資料上界面上下巖電特征發(fā)生突變。

準(zhǔn)層序（高頻層序）以測井曲線降頻后的極值點(diǎn)或拐點(diǎn)為界（圖3），不整合面上下也意味著沉積環(huán)境的變化，由水退到水進(jìn)，或者由上一層序的高位域變?yōu)橄乱粚有虻牡臀挥?，在測井曲線上也有明顯的響應(yīng)，一般也是測井曲線的變化拐點(diǎn)，因此在降頻曲線上也可以體現(xiàn)。這樣不同級次降頻曲線的極值點(diǎn)可以反映高精度層序界面，其降頻曲線旋回反映高精度層序的旋回性變化。通過對多條曲線的多次降頻處理，然后根據(jù)降頻曲線的旋回性，并結(jié)合縱向巖性組合和沉積旋回，就可以實現(xiàn)真正意義上的將細(xì)粒沉積的高頻旋回識別出來。

圖3 曙112井測井降頻高精度層序劃分方案

在井上高頻層序劃分的基礎(chǔ)上，通過制作地震合成記錄，在連井地震剖面上進(jìn)行精細(xì)層位標(biāo)定，為了質(zhì)保，研究區(qū)共建連井井震剖面11條，標(biāo)定探井64口；在井震標(biāo)定的基礎(chǔ)上開展高頻層序剖面解釋工作，最終完成高頻層序格架建立；然后進(jìn)行初始模型的建立和模型插值，通過模型插值檢驗高頻層序格架建立的正確性，并反復(fù)修改閉合，最終使得插值模型達(dá)到縱橫向能量均衡、過渡自然、符合地質(zhì)規(guī)律。

2.2.4 波形指示模擬反演

地震波形指示反演的基本思想是在等時格架約束下，利用地震波形的橫向變化表征儲層的空間變異性，其優(yōu)勢是地震的橫向密集性［4］。

將連井剖面波形指示模擬反演結(jié)果與實際鉆井對比，反復(fù)調(diào)節(jié)反演參數(shù)，直到達(dá)到最佳效果，然后開展全區(qū)波形指示模擬反演。從過曙134井反演剖面（圖4）中可以看出，反演效果較好。首先相對于普通波阻抗反演剖面（圖5），在一定程度上提高了縱向分辨率，且與鉆井砂體吻合性較好；其次提高了薄層砂體的預(yù)測精度，平面上可以實現(xiàn)儲層初步預(yù)測。反演體的平面刻畫對儲層預(yù)測的效果至關(guān)重要，本文對采取阻抗體小閥值、多種子點(diǎn)追蹤方法取得較好預(yù)測效果。為了能更好地將砂泥巖在縱向上分離，并使同一砂體在橫向上保持連續(xù)，在對已有井不同巖性阻抗值分析的基礎(chǔ)上，最終確定閥值為200；通過閥值滑動進(jìn)行種子點(diǎn)追蹤，同時將追蹤出來的阻抗體分布與實際井巖性進(jìn)行比對，保留能反應(yīng)真實砂體的阻抗種子點(diǎn)，有效地提高砂巖識別的準(zhǔn)確率；最后對篩選的砂體時間厚度進(jìn)行時深轉(zhuǎn)化，使其達(dá)到預(yù)測砂巖分布的目的。

2.3 應(yīng)用效果

在上述等時相控波形指示模擬反演的基礎(chǔ)上，根據(jù)區(qū)帶評價結(jié)果，針對地質(zhì)目標(biāo)再次進(jìn)行砂體定量刻畫［9］，并與地質(zhì)規(guī)律相對比，最終實現(xiàn)地質(zhì)和地震綜合落實薄砂巖空間展布。通過與完鉆探井驗證發(fā)現(xiàn)，研究區(qū)薄砂巖儲層預(yù)測符合率約為80%。當(dāng)砂體厚度小于8.0 m時，其準(zhǔn)確率略低，平均為70%（表1）；當(dāng)砂體厚度大于8.0 m時，其準(zhǔn)確率較高，平均為83%（表1）。

曙134井區(qū)杜Ⅱ砂組預(yù)測砂巖分布在曙134井東南部（圖6），平面上呈北東-南西向展布，預(yù)測砂巖厚度一般在10.0～25.0 m，其中曙118井、曙127井、曙111井、曙112井、曙134井厚度吻合率大于85%（表1）。預(yù)測砂巖厚度高值區(qū)呈斑狀分布，且呈多排斜列分布特征，可能反映砂組內(nèi)部的灘壩遷移；每排內(nèi)又有2～3個次級主體，因此將灘壩砂體沉積微相按照砂體厚度、主次、微地貌劃分為灘壩主體、灘壩側(cè)翼和灘壩間3類（圖7）。根據(jù)砂地比給予定量描述，灘壩主體砂地比一般大于15%，灘壩側(cè)翼砂地比一般5%～15%，灘壩間砂地比小于5%。曙134井區(qū)杜Ⅱ砂組發(fā)育4個灘壩主體，其中兩個灘壩主體在構(gòu)造高部位有兩口井獲得工業(yè)油流，6口井油氣顯示良好，因此曙134井所在灘壩主體和曙111井西側(cè)灘壩主體是下步挖潛的有利區(qū)，累計面積11.5 km2。

圖7 曙134井區(qū)杜Ⅱ砂組沉積微相

表1 曙北地區(qū)等時相控SM I反演儲層預(yù)測與實際鉆遇砂巖厚度對比

3 結(jié)論

（1）等時相控波形指示模擬反演是資料品質(zhì)低、巖性復(fù)雜、斷裂較發(fā)育地區(qū)薄砂巖儲層預(yù)測的有效手段，技術(shù)思路包括正演模擬分析、敏感曲線重構(gòu)、等時相控建模、波形指示模擬反演和種子點(diǎn)追蹤等。其中，等時相控建模是關(guān)鍵，高精度層序劃分、明確的地質(zhì)認(rèn)識和精細(xì)的三維地震解釋可以使儲層預(yù)測結(jié)果更為可靠、精度更高。

（2）曙北地區(qū)沙四段薄砂巖儲層預(yù)測結(jié)果表明，等時相控波形指示模擬反演薄砂巖預(yù)測符合率約為80%，當(dāng)砂體厚度小于8.0 m時，其準(zhǔn)確率略低；當(dāng)砂體厚度大于8.0 m時，其準(zhǔn)確率較高。預(yù)測結(jié)果與鉆井及地質(zhì)規(guī)律吻合度高，可有效指導(dǎo)該區(qū)勘探部署。

（3）井震結(jié)合分析灘壩砂體發(fā)育規(guī)律表明，濱淺湖灘壩沉積可細(xì)分為灘壩主體、灘壩側(cè)翼和灘壩間三類，灘壩砂體縱向疊置、遷移可形成大面積砂巖分布，構(gòu)造高部位的灘壩主體是最有利的勘探領(lǐng)域。