羅釧江， 羅耀華， 魯建隆， 殷曄宗， 胡 鑫

(1.成都理工大學(xué) 能源學(xué)院，成都 610059；2.中國石油集團(tuán) 渤海鉆探井下作業(yè)分公司，任丘 062552)

0 研究區(qū)概況

洛帶遂寧組氣藏構(gòu)造以北東低南西高的鼻狀背斜為特征，在西南構(gòu)造軸上形成了低振幅局部高點(diǎn)，局部低鞍且閉合面積受限。氣藏類型一般是低-特低孔隙度、特低滲儲層、細(xì)孔-小喉孔隙型三角洲前緣砂巖、構(gòu)造-巖性復(fù)合型圈閉氣藏，該氣藏?zé)o邊緣、底部有水，是一種低豐度、中低產(chǎn)量的恒定容量封閉式彈性氣驅(qū)驅(qū)動(dòng)類型。該研究區(qū)沉積環(huán)境主要為三角洲前緣亞相，具有非均質(zhì)性強(qiáng)、垂向上具有多層疊置的正粒序、儲層薄、相變快等特征 (圖1) 。

圖1 區(qū)域構(gòu)造位置圖Fig.1 Location map of regional tectonics

圖2 Earthmodel模塊建立的低頻模型Fig.2 Low frequency model established by Earthmodel module

圖3 洛帶地區(qū)疊后波阻抗剖面Fig.3 Post stack impedance profile in Luodai area

1 地質(zhì)統(tǒng)計(jì)學(xué)反演原理

疊后地震反演的方法主要分兩類：①基于物性的波阻抗反演和孔隙度反演；②基于巖性的GR擬聲波反演和地質(zhì)統(tǒng)計(jì)學(xué)反演。

1.1 約束稀疏脈沖反演

約束稀疏脈沖反演是地質(zhì)統(tǒng)計(jì)學(xué)的基礎(chǔ), 波阻抗反演采用稀疏脈沖反演方法[2]。構(gòu)建初始波阻抗模型，將測井資料、鉆井?dāng)?shù)據(jù)和地震資料結(jié)合，在EarthModel模塊里面建立低頻初始模型。該模型能夠較好地反映沉積體地質(zhì)特征，在建立初始的模型過程中，將在縱向上具有高分辨率測井資料與在橫向上具有連續(xù)變化特征的地震界面結(jié)合。建立初始波阻抗模型的原因有兩方面：①在地震反演過程中控制解的范圍，保證結(jié)果的穩(wěn)定性和收斂性；②大大彌補(bǔ)了地震資料的低頻部分(圖2)。

稀疏脈沖反演的步驟：①遵從稀疏的原則在地震道中提取的反射系數(shù)與子波褶積合成新的地震道；②依據(jù)原始的地震道修改反射系數(shù)；③循環(huán)步驟①、步驟②，直至得到最逼近地震道的反射系數(shù)序列；④通過得到的反射系數(shù)序列，在測井資料的約束下得到一個(gè)全頻帶的絕對波阻抗[3]。反演結(jié)果如圖3所示。

1.2 疊后地質(zhì)統(tǒng)計(jì)學(xué)反演

疊后地質(zhì)統(tǒng)計(jì)學(xué)反演是在地震反演當(dāng)中引入了序貫擬合的思路[4]。首先從約束稀疏脈沖反演得到波阻抗，掌握儲層大致的展布特征，以此為依據(jù)求取差變函數(shù)；然后以測井資料為依據(jù)，以原始地震數(shù)據(jù)為約束，隨機(jī)模擬出井間的波阻抗并轉(zhuǎn)換成反射系數(shù)，再與子波褶積合成新的地震道，通過反復(fù)迭代，達(dá)到最佳逼近原始地震道的程度[3,8,11]。

2 地質(zhì)統(tǒng)計(jì)學(xué)反演過程

地質(zhì)統(tǒng)計(jì)學(xué)反演的過程主要包括：①巖性劃分；②反演參數(shù)分析；③巖性模擬；④孔隙度模擬4個(gè)部分[6]。

圖4 洛帶遂寧組氣藏GR和阻抗的交匯圖Fig.4 Intersection diagram of GR and impedance of Suining formation gas reservoir in Luodai

2.1 巖性劃分

本次項(xiàng)目巖性劃分是在對工區(qū)內(nèi)鉆井分析的基礎(chǔ)上開展的，通過對工區(qū)內(nèi)100余口井的巖性與GR和阻抗的交匯分析(圖4)，泥巖主要表現(xiàn)為高GR低阻抗特征，而砂巖主要表現(xiàn)為低GR高阻抗特征，當(dāng)砂巖含氣或空隙較好時(shí)則表現(xiàn)為低GR低阻抗特征，巖性的劃分則可以根據(jù)GR與阻抗雙重參數(shù)最終確定。

2.2 反演參數(shù)分析

地質(zhì)統(tǒng)計(jì)學(xué)反演的精確程度主要受反演參數(shù)約束，反演參數(shù)包括概率密度函數(shù)、差變函數(shù)等。其中，變差函數(shù)反映儲層的空間結(jié)構(gòu)，就地質(zhì)意義來講，縱向變程對縱向分辨率有影響，主要是通過測井?dāng)?shù)據(jù)來求??；橫向變程反映了目的層的橫向展布特征，需要對地質(zhì)沉積的規(guī)律并結(jié)合分布的泥巖或煤層反演的測試結(jié)果來確定。

變差函數(shù)的公式為：

(1)

式中：h為滯后距；Z為區(qū)域變量；N為距離為h的個(gè)數(shù)；R為變差函數(shù)。

圖5 變差函數(shù)參數(shù)示意圖Fig.5 Schematic diagram of variogram parameters

圖5中,α為r(h)達(dá)到平穩(wěn)時(shí)滯后距離的值(變程),其表示空間中的最大相關(guān)距離;C0為一常數(shù),表示區(qū)域化變量隨機(jī)變化; 滯后值α表示在空間分布上區(qū)域化變量的平均尺度，滯后值越大，在空間分布上區(qū)域化變量相關(guān)性越大，變化的速度也就越慢，隨機(jī)性越小。C為拱高,代表區(qū)域化變量的結(jié)構(gòu)性變化;C0+C為基臺值,代表著區(qū)域化變量的最大變化范圍。當(dāng)C0+C保持不變時(shí),C0變大,則C變小,表示在空間分布上區(qū)域化變量的隨機(jī)性變強(qiáng),則區(qū)域化變量的結(jié)構(gòu)性越差;反之如果C0+C保持不變C變大,則C0變小,結(jié)果也將相反[7]。

筆者對研究區(qū)內(nèi)100余口測井曲線質(zhì)量好的井進(jìn)行分析，阻抗概率分布嚴(yán)格服從高斯分布，如圖6所示。中值為11 069.77(kg/cm3*m/s)，偏差為538.185 4(kg/ cm3*m/s)2；變差函數(shù)的求取來自于井震擬合，在擬合的過程中需要確定兩個(gè)方向的變程，即縱向變程以及橫向變程。由于井的數(shù)量較多，在縱向上的變程也就比較精確，再通過比較高斯型和指數(shù)型擬合，該地區(qū)指數(shù)型擬合較好，縱向變程為6 m。

圖6 洛帶遂寧組氣藏波阻抗概率分布圖Fig.6 Probability distribution of wave impedance in Suining formation gas reservoir, Luodai

圖7 洛帶遂寧組氣藏連井巖性反演剖面Fig.7 Lithology inversion section of Suining formation gas reservoir in Luodai

圖8 洛帶遂寧組氣藏砂組孔隙度與縱波阻抗交匯關(guān)系圖Fig.8 Relationship between porosity and P-wave impedance of Jsn10 and Jsn11-1 sand formations in Suining formation gas reservoir, Luodai

圖9 洛帶遂寧組氣藏砂組孔隙度與縱波阻抗交匯關(guān)系圖Fig.9 Relationship between porosity and P-wave impedance of Jsn11-2 and Jsn12 sand formations in Suining formation gas reservoir, Luodai

圖10 洛帶地區(qū)連井孔隙度反演剖面Fig.10 Inversion profile of continuous well porosity in Luodai area

圖11 地震巖性反演砂體厚度預(yù)測與實(shí)鉆對比Fig.11 Seismic lithology inversion sand body thickness prediction and real drilling contrast

2.3 巖性模擬

地質(zhì)統(tǒng)計(jì)學(xué)反演具有很強(qiáng)的識別薄層和精確預(yù)測儲層的能力[1,3,9]。該反演技術(shù)在水平方向上使用多網(wǎng)格蒙特卡洛模擬，在縱向時(shí)間域中使用馬爾科夫鏈蒙特卡洛模擬方法。通過對地震資料、地震子波、巖性分布、地層格架、差變函數(shù)，統(tǒng)計(jì)學(xué)參數(shù)的優(yōu)選、新井檢驗(yàn)等進(jìn)行反演預(yù)測[5,13]。通過20次的巖性模擬，再根據(jù)最新的鉆井成果與模擬出的結(jié)果對比分析，從中優(yōu)選模擬結(jié)果進(jìn)行極大似然巖性體求取，最終得到本區(qū)的極大似然巖性體，為后續(xù)的砂體厚度求取和儲層厚度求取提供借鑒(圖7)。

2.4 孔隙度模擬

3 反演結(jié)果校驗(yàn)

通過進(jìn)行地質(zhì)統(tǒng)計(jì)學(xué)反演，最終得到極大似然巖性體，在巖性體上進(jìn)行各砂體頂、底面的精細(xì)解釋。為了驗(yàn)證反演的平面展布特征的準(zhǔn)確度，開展了預(yù)測砂體與實(shí)際鉆井砂巖厚度的對比分析，通過圖11的對比來看，砂巖厚度在正負(fù)1.5 m誤差范圍內(nèi)到達(dá)81.3%，這表明儲層預(yù)測結(jié)果可靠，可以應(yīng)用到后期的各項(xiàng)的研究。

4 結(jié)論

本次研究旨在刻畫精確有利儲層，優(yōu)化了地震資料，使儲層的地質(zhì)模型在空間分布上更加合理。也以此了解到要優(yōu)化地震資料，只用一種地震反演方法遠(yuǎn)遠(yuǎn)不能達(dá)到精細(xì)的預(yù)測結(jié)果。雖然通過約束稀疏脈沖反演將界面形地震剖面轉(zhuǎn)換為地層形剖面，提出了圍巖變化對于儲層地震響應(yīng)特征的干擾，但是在縱向分辨率上得不到精確的效果。地質(zhì)統(tǒng)計(jì)學(xué)反演是一種概率隨機(jī)反演，通過隨機(jī)模擬和優(yōu)化隨機(jī)模擬的結(jié)果組成，并且利用了隨機(jī)建模和地震反演的特點(diǎn)[6]，極大地提高了垂向分辨率。但是它同樣需要以約束稀疏脈沖反演方法得到的波阻抗來校正，通過不斷完鉆的新井，來檢驗(yàn)預(yù)期的成果，才能更有效地識別厚度薄的砂體，更精確地細(xì)刻畫地震相帶[12]，建立的地質(zhì)模型才能更加接近實(shí)際地質(zhì)情況。