黃玉 郭娟娟 卞保力 王彥君 劉海磊
摘? ? 要:準(zhǔn)噶爾盆地瑪東斜坡區(qū)深層石炭系溢流相火山巖是主要儲(chǔ)層之一。由于火山巖儲(chǔ)層埋藏深、巖性巖相變化快、縱橫向非均質(zhì)性強(qiáng)、深層探井少等特征,導(dǎo)致常規(guī)火山巖巖相預(yù)測識(shí)別難度大,精度不高。為更精確預(yù)測溢流相火山巖的分布,首次在深層火山巖領(lǐng)域應(yīng)用地震波形指示反演技術(shù)進(jìn)行巖相預(yù)測。采用井震結(jié)合,挖掘相似波形對(duì)應(yīng)的測井曲線的共性結(jié)構(gòu)信息,充分利用地震波形的橫向變化,開展高頻成分估計(jì),獲得對(duì)儲(chǔ)層巖性、物性反映較敏感的高分辨率反演數(shù)據(jù)體。經(jīng)鉆井證實(shí),該方法解決了研究區(qū)資料反演結(jié)果不理想、模型化嚴(yán)重、精確度不高的問題。提高了高頻成分的確定性,準(zhǔn)確性更高,使反演結(jié)果從完全隨機(jī)走向逐步確定,為井位部署和儲(chǔ)量提交提供可靠依據(jù)。
關(guān)鍵詞:地震波形指示反演;深層火山巖;溢流相;巖相預(yù)測;瑪東斜坡區(qū)
準(zhǔn)噶爾盆地石炭系深層火山巖分布廣泛,緊鄰多個(gè)生烴中心,整體勘探潛力大,是盆地深層油氣發(fā)現(xiàn)的潛在領(lǐng)域[1]。前人通過多年研究,主要在盆地邊緣、中淺層領(lǐng)域發(fā)現(xiàn)多個(gè)大中型火山巖油氣田[2-3]。受地震采集處理、儲(chǔ)層預(yù)測技術(shù)等條件限制,深層火山巖領(lǐng)域勘探程度較低,鉆井少且分布不均。火山巖本身巖性復(fù)雜,火山巖巖性識(shí)別困難,多解性強(qiáng)、精度不高。
前人對(duì)深層火山巖巖相進(jìn)行了稀疏脈沖反演,地質(zhì)統(tǒng)計(jì)學(xué)反演等一系列反演工作[4-7],反演結(jié)果不理想、模型化嚴(yán)重、精確度不高。究其原因,常規(guī)反演方法具有一定的適用條件。稀疏脈沖反演方法以地震為主,分辨率低,僅能反演波阻抗,無法對(duì)地質(zhì)形態(tài)復(fù)雜的火山巖進(jìn)行有效預(yù)測[4-7]。地質(zhì)統(tǒng)計(jì)學(xué)反演方法反演結(jié)果隨機(jī)性較強(qiáng)、橫向分辨率較低[8-10],對(duì)井的數(shù)量及井點(diǎn)分布均勻要求較高,適應(yīng)于評(píng)價(jià)-開發(fā)階段,對(duì)勘探程度低、井點(diǎn)少且分布不均的研究區(qū)反演效果差。地震波形指示反演為新興的地震儲(chǔ)層預(yù)測方法,該方法的特點(diǎn)在于利用地震波形橫向變化表征儲(chǔ)層空間的變異性,對(duì)井點(diǎn)分布不做要求,并充分結(jié)合地震波形的橫向變化開展高頻成分估計(jì),建立的插值模型更符合地質(zhì)沉積規(guī)律。杜偉維、高君、韓長城等在薄砂層、灘壩砂、砂礫巖儲(chǔ)層預(yù)測中運(yùn)用該方法[11-15],提高了縱向識(shí)別精度,降低了反演結(jié)果的隨機(jī)性。本次研究首次在深層溢流相火山巖領(lǐng)域應(yīng)用該反演技術(shù),進(jìn)行巖相預(yù)測,期望提高深層溢流相火山巖預(yù)測精度,為井位部署和儲(chǔ)量提交提供可靠依據(jù)。
1? 方法原理
1.1? 原理及優(yōu)勢性分析
地震數(shù)據(jù)體反映沉積環(huán)境和巖性組合的空間變化,地震波形包含的振幅、相位、頻率等信息可較為全面的代表儲(chǔ)層特征。波形指示反演充分利用地震數(shù)據(jù)與測井?dāng)?shù)據(jù)建立初始地質(zhì)模型,通過相關(guān)參數(shù)優(yōu)選,不斷地進(jìn)行波形指示反演,將反演結(jié)果QC質(zhì)控與實(shí)際地震結(jié)果不斷對(duì)比,直至獲得最優(yōu)的反演結(jié)果[15],特別適用于高精度預(yù)測非均質(zhì)性強(qiáng)的火山巖儲(chǔ)層。公示為:
[Yx0= k=1nλkYxk]…(1)
式中:[Yx0]——未知點(diǎn)的值;
[Yxk]——波形優(yōu)選的已知樣本點(diǎn)的值;
[λk]——第k個(gè)已知樣本點(diǎn)對(duì)未知樣點(diǎn)的
權(quán)重;
[k]——優(yōu)選樣本點(diǎn)的個(gè)數(shù)。
反演過程中反演參數(shù)的選擇直接影響反演結(jié)果,主要參數(shù)包括:有效樣本數(shù)、最佳截止頻率。有效樣本數(shù)表征地震波形空間變化對(duì)儲(chǔ)層的影響程度,該參數(shù)的確定參照對(duì)已知井的統(tǒng)計(jì)結(jié)果。隨著樣本數(shù)增加,相關(guān)指數(shù)逐漸增大,在相關(guān)指數(shù)不再增大時(shí)所對(duì)應(yīng)的樣本數(shù)即為確定的有效樣本數(shù)。實(shí)際應(yīng)用中,有效樣本數(shù)的優(yōu)選除參考質(zhì)控窗口的統(tǒng)計(jì)結(jié)果外,還要考慮沉積環(huán)境、井?dāng)?shù)的影響。如儲(chǔ)層橫向變化快,可適當(dāng)減小有效樣本數(shù),如井?dāng)?shù)少應(yīng)適當(dāng)減小樣本數(shù)。
最佳截止頻率值不同,反演結(jié)果分辨率不同。隨著最佳截止頻率的增加,反演結(jié)果的分辨率逐漸增大,剖面的隨機(jī)性也逐漸增強(qiáng)。該參數(shù)與有效樣本數(shù)相關(guān),在相關(guān)指數(shù)不再增大時(shí)所對(duì)應(yīng)的頻率即為最大有效頻率,頻率成分以外的都是隨機(jī)的。其值的確定除了參考質(zhì)控結(jié)果,還應(yīng)考慮兩方面內(nèi)容:一是如偏向于反演的確定性,可適當(dāng)減小該參數(shù)值;若偏向反演的分辨率,可適當(dāng)增大該參數(shù)值;二是參考目標(biāo)體厚度,目標(biāo)體厚度薄,應(yīng)設(shè)置較高的截止頻率,如目標(biāo)體厚度小,可適當(dāng)減小該參數(shù)。
地震波形指示反演是一種基于相控模型的高精度反演方法,在地震波形特征指導(dǎo)下對(duì)反射系數(shù)組合尋優(yōu),反演結(jié)果具“低頻確定、高頻隨機(jī)”特點(diǎn)[5]。該反演方法的優(yōu)勢體現(xiàn):一是加入了地震波形作為先驗(yàn)信息,更好地體現(xiàn)了相控的思想;二是針對(duì)深層火山巖地震資料成像信噪比低、分辨率低的特點(diǎn),利用高頻測井曲線,建立地震波形結(jié)構(gòu)與高頻測井曲線結(jié)構(gòu)的映射關(guān)系,提高了反演結(jié)果的縱向分辨率和橫向分辨率,使反演結(jié)果從完全隨機(jī)到逐步確定;三是對(duì)井位分布均勻性的要求遠(yuǎn)低于常規(guī)的地震統(tǒng)計(jì)學(xué)反演,很大程度上提高了儲(chǔ)層預(yù)測的精度和適用領(lǐng)域。
1.2? 實(shí)現(xiàn)流程
結(jié)合前人地震波形指示反演在其他領(lǐng)域的研究思路及流程[14-15],據(jù)地震波形指示反演原理方法,考慮深層溢流相火山巖儲(chǔ)層巖性以玄武巖、安山巖、流紋巖為主,電阻率較爆發(fā)相、火山沉積相高,地震反射外形呈楔狀、丘狀、席狀和透鏡狀,內(nèi)部反射特征是傾斜層狀反射,連續(xù)性好,強(qiáng)-中振幅,平行-亞平行反射結(jié)構(gòu)等特征。本文制定了針對(duì)深層溢流相火山巖波形指示反演技術(shù)流程(圖1)。關(guān)鍵技術(shù)步驟有4步:第一,據(jù)地震波形特征對(duì)已知井進(jìn)行分析,與待判別道地震波形進(jìn)行類比并排序,優(yōu)選出波形相似的井作為樣本井建立初始模型;第二,據(jù)不同火山巖巖性巖相測井曲線的差異性,對(duì)樣本井曲線進(jìn)行分頻分析,統(tǒng)計(jì)溢流相火山巖測井曲線結(jié)構(gòu)特征;第三,利用樣本井波阻抗結(jié)構(gòu)特征建立初始模型,參照鄰井對(duì)確定性以外的高頻成分進(jìn)行模擬,使模擬結(jié)果的結(jié)構(gòu)關(guān)系既符合地震中頻阻抗特征又符合井曲線結(jié)構(gòu)特征;第四,以地震中頻為依據(jù),在貝葉斯框架下將初始模型的高頻成分予以保留,不符合的進(jìn)行過濾,通過反復(fù)隨機(jī)模擬,確保隨機(jī)模擬中頻成分與地震一致,同時(shí)保留各井樣本的確定性結(jié)果成分,輸出最終模擬結(jié)果。
2? 應(yīng)用實(shí)例
2.1? 實(shí)例區(qū)概況
瑪東斜坡區(qū)石炭—二疊系發(fā)育繼承性古隆起,形成一系列構(gòu)造鼻?。▓D2)?,敽枷菔桥璧赜蜌飧患潭茸罡叩纳鸁N凹陷,具有充足的油氣源供應(yīng)[17-18],同時(shí)溝通有緣的深部斷裂的發(fā)育,使實(shí)例區(qū)鼻隆成為油氣藏聚集成藏有利區(qū)。但實(shí)例區(qū)目前現(xiàn)狀是勘探目的層深,一般超過5 000 m,石炭系三維地震資料主頻為28 Hz,有效頻寬4~52 Hz?;鹕綆r輪廓明顯,構(gòu)造特征清楚,頂界及內(nèi)部層序界面清晰。研究區(qū)內(nèi)井的測井項(xiàng)目有雙側(cè)向、自然伽馬、巖石密度、補(bǔ)償中子、聲波等,采集及應(yīng)用質(zhì)量為優(yōu),基本達(dá)到儲(chǔ)層預(yù)測要求。已鉆井11口,其中僅1口井獲工業(yè)油流,勘探程度極低。結(jié)合前人研究成果[19],分析石炭系新鉆井巖心、孔滲數(shù)據(jù),認(rèn)為實(shí)例區(qū)石炭系火山巖優(yōu)質(zhì)儲(chǔ)層為溢流相安山巖,次為爆發(fā)相火山角礫巖和凝灰?guī)r、溢流相流紋巖。
2.2? 關(guān)鍵技術(shù)步驟
利用已鉆井巖心、分析化驗(yàn)、測井、地震波形特征等資料,據(jù)不同巖性對(duì)應(yīng)的不同測井響應(yīng)值,建立不同巖性測井曲線交匯圖,分析各類巖性所對(duì)應(yīng)的敏感曲線特征,溢流相的測井響應(yīng)表現(xiàn)為明顯的高電阻率特征。在此基礎(chǔ)上,優(yōu)選能區(qū)分溢流相火山巖的敏感測井參數(shù)。通過精細(xì)時(shí)深標(biāo)定,分析溢流相火山巖地震波形特征。通過反演參數(shù)優(yōu)選,進(jìn)行地震波形指示反演,預(yù)測深層溢流相火山巖。
2.2.1? 敏感測井參數(shù)優(yōu)選
測井資料信息豐富,能反映巖性、巖相特征[20-21]。據(jù)分析資料較全的已鉆井取芯段不同巖性對(duì)應(yīng)測井曲線響應(yīng)特征,建立石炭系各類巖性聲波時(shí)差(AC)-密度(DEN)交匯圖和波阻抗(IMP)-深側(cè)向電阻率(RT)交匯圖(圖3),明確不同巖相巖性測井響應(yīng)特征。由AC-DEN 交匯圖可看出,溢流相英安巖、流紋巖、爆發(fā)相火山角礫巖、凝灰?guī)rAC數(shù)據(jù)區(qū)間重疊較多,難以區(qū)分。從IMP-RT交匯圖看出,深側(cè)向電阻率(RT)對(duì)不同巖相巖性特征較敏感,溢流相火山巖與上下圍巖電阻率差異較大,可較好地區(qū)分巖相巖性。RT>100 Ω·m為溢流相安山巖,40 Ω·m