[發(fā)明專利]基于地質(zhì)體連續(xù)性建模的高精度井拓反演方法

有效

• 申請?zhí)枺?/td>: 201811269424.2
• 申請日: 2018-10-29
• 公開（公告）號: CN109471163B
• 公開（公告）日: 2020-07-21
• 發(fā)明（設(shè)計(jì)）人: 田立新;王波;周東紅;明君;夏同星;張志軍;劉壘;沈東義;李福強(qiáng);薛明星;周學(xué)鋒
• 申請（專利權(quán)）人: 中國海洋石油集團(tuán)有限公司;中海石油（中國）有限公司天津分公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 天津市北洋有限責(zé)任專利代理事務(wù)所 12201
• 代理人: 杜文茹

一種基于地質(zhì)體連續(xù)性建模的高精度井拓反演方法，基于地質(zhì)體連續(xù)性變化的認(rèn)識，使用地質(zhì)體連續(xù)性模型建模方法和基于貝葉斯理論的波形反演方法，首先根據(jù)井點(diǎn)處測井資料計(jì)算得到阻抗數(shù)據(jù)，井點(diǎn)處的地震道和建立的精確阻抗模型作為輸入，使用貝葉斯波形反演方法得到井點(diǎn)處地震道的阻抗，反演結(jié)果根據(jù)層位解釋數(shù)據(jù)作相應(yīng)時移后又作為相鄰道的初始模型，結(jié)合相鄰道的地震數(shù)據(jù)，得到相鄰道的反演結(jié)果，這樣逐道反演下去，得到最終反演結(jié)果。該技術(shù)提高了初始輸入模型的準(zhǔn)確度，反演結(jié)果精度高，信噪比高，可用于薄互層砂體、厚砂層，以及地震資料品質(zhì)較差處的砂體的刻畫。

技術(shù)領(lǐng)域

本發(fā)明涉及一種用于石油勘探地震數(shù)據(jù)處理與解釋的反演方法。特別是涉及一種基于地質(zhì)體連續(xù)性建模的高精度井拓反演方法。

背景技術(shù)

地震反演作為儲層預(yù)測的核心技術(shù)，能夠?qū)y井縱向上的高分辨率和地震橫向上的密集采樣這兩者的優(yōu)勢結(jié)合起來，估算出地層的巖性、物性特征信息橫向上的變化，被廣泛應(yīng)用于石油勘探開發(fā)的各個階段。目前常用的方法有疊后地震波阻抗反演、疊前彈性參數(shù)反演、地質(zhì)統(tǒng)計(jì)學(xué)反演、擬測井曲線反演等。疊后波阻抗反演和疊前彈性參數(shù)反演首先利用地震解釋層面和斷層資料建立構(gòu)造模型，在構(gòu)造模型上利用測井資料插值建立屬性模型，利用屬性模型結(jié)合地震資料反演得到波阻抗資料和彈性資料，建立的屬性模型提供了反演結(jié)果的低頻信息，實(shí)踐看來，建模提供的低頻信息并不能反應(yīng)薄互儲層的劇烈變化，反演結(jié)果不理想，擬測井曲線反演求取地震資料與測井資料(比如GR)的映射關(guān)系，再利用映射關(guān)系將地震資料轉(zhuǎn)化為屬性資料(比如GR體)，反演得到的結(jié)果有明確的意義，但是，由于測井是有限的，因此，映射關(guān)系的準(zhǔn)確性不確定。地質(zhì)統(tǒng)計(jì)學(xué)反演在建立的構(gòu)造模型的基礎(chǔ)上，對測井資料進(jìn)行插值，同時應(yīng)用地震資料進(jìn)行控制，得到高分辨率的資料，該方法需要井網(wǎng)密度大，同時要求時深標(biāo)定特別準(zhǔn)確，應(yīng)用中工作難度大。綜上所述，目前地震反演方法普遍存在著建模精度不夠，反演結(jié)果精度不高，不能滿足目前油氣生產(chǎn)中對儲層精細(xì)刻畫的要求的問題。

發(fā)明內(nèi)容

本發(fā)明所要解決的技術(shù)問題是，提供一種通過提高初始模型建模的精度，提高反演資料的精度，實(shí)現(xiàn)儲層精細(xì)刻畫的基于地質(zhì)體連續(xù)性建模的高精度井拓反演方法。

本發(fā)明所采用的技術(shù)方案是：一種基于地質(zhì)體連續(xù)性建模的高精度井拓反演方法，包括如下步驟：

1)將井點(diǎn)處的地震道作為當(dāng)前反演道，在井點(diǎn)處利用測井資料計(jì)算得到縱波阻抗模型，把井點(diǎn)處的地震道和得到的縱波阻抗模型作為輸入，使用貝葉斯波形反演方法計(jì)算得到井點(diǎn)處地震道的縱波阻抗；

2)根據(jù)反演目標(biāo)層段的層位解釋成果讀取井點(diǎn)處地震道的時間值和與井點(diǎn)處地震道最近的地震道的時間值，將井點(diǎn)處地震道的時間值和與所述最近地震道的時間值作差得到時移量，根據(jù)得到的該時移量對當(dāng)前反演道的縱波阻抗作相應(yīng)移動；

3)將步驟2)中所述與井點(diǎn)處地震道最近的地震道作為當(dāng)前反演道，將步驟2)中移動后的縱波阻抗作為當(dāng)前反演道的初始模型，結(jié)合采集到的當(dāng)前反演道的地震數(shù)據(jù)，使用貝葉斯波形反演方法計(jì)算得到當(dāng)前反演道的縱波阻抗；

4)根據(jù)反演目標(biāo)層段的層位解釋成果讀取當(dāng)前反演道的時間值和與當(dāng)前反演道最近的地震道的時間值，將當(dāng)前反演道的時間值和與當(dāng)前反演道最近的地震道的時間值作差得到時移量，根據(jù)得到的該時移量對當(dāng)前反演道的縱波阻抗作相應(yīng)移動；

5)將步驟4)中所述與當(dāng)前反演道最近的地震道作為新的當(dāng)前反演道，將步驟4)中移動后縱波阻抗作為該新的當(dāng)前反演道的初始模型，結(jié)合采集到的該新的當(dāng)前反演道的地震數(shù)據(jù)，使用貝葉斯波形反演方法計(jì)算得到該新的當(dāng)前反演道的縱波阻抗；

6)重復(fù)執(zhí)行步驟4)和步驟5)，直到工區(qū)內(nèi)所有地震道計(jì)算完畢，分別得到所有地震道的縱波阻抗。

步驟1)、步驟2)和步驟5)所述的使用貝葉斯波形反演方法計(jì)算包括：

(1)建立正演模型