一種疊后地震波阻抗反演方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種疊后地震波阻抗反演方法，屬于地球物理勘探技術(shù)領(lǐng)域。
【背景技術(shù)】
[0002]隨著石油勘探開發(fā)的不斷深入，基于地震同相軸追蹤進(jìn)行層位解釋的常規(guī)方法不能很好地滿足儲層識別和解釋的地質(zhì)需求。借助于地震波阻抗反演提高地震資料識別儲層的能力已經(jīng)成為儲層預(yù)測和油藏描述的重要技術(shù)手段。地震反演是利用地表觀測的地震資料，以已知地質(zhì)規(guī)律和鉆井、測井資料為約束，對地下巖層空間結(jié)構(gòu)和物理性質(zhì)進(jìn)行認(rèn)知的過程。波阻抗反演是指利用地震資料反演地層波阻抗的地震特殊處理解釋技術(shù)，其具有明確的物理意義，是儲集層巖性預(yù)測、油藏特征描述的確定性方法，在實(shí)際應(yīng)用中取得了顯著的地質(zhì)效果，因此地震反演通常特指波阻抗反演。
[0003]地震反演通常分為疊前和疊后反演兩大類。近30年來，疊后地震反演取得了巨大進(jìn)展，已形成了多種成熟技術(shù)，尤其是基于模型地震反演方法由于其精度和分辨率較高，目前得到了廣泛應(yīng)用。基于模型地震反演實(shí)質(zhì)上是地震-測井聯(lián)合反演，其結(jié)果的低、高頻信息來源于測井資料，構(gòu)造特征及中頻段取決于地震數(shù)據(jù)。多解性是基于模型地震反演的固有特性，即地震有效頻帶以外的信息不會影響合成地震資料的最終結(jié)果，減小基于模型方法多解性問題的關(guān)鍵在于正確建立初始模型。
[0004]地球物理學(xué)家近年來對于波阻抗反演結(jié)果中高頻分量的多解性探討較多，對于地層沉積壓實(shí)效應(yīng)帶來的極低頻分量的巖性識別的多解性則往往忽視。地震波阻抗是進(jìn)行巖性識別和儲層預(yù)測的主要依據(jù)，然而由于受地層沉積的壓實(shí)效應(yīng)影響，地震波阻抗不能真實(shí)反映巖性變化。

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

[0005]本發(fā)明的目的是針對現(xiàn)有疊后地震波阻抗反演方法受地層沉積壓實(shí)效應(yīng)影響大、利用波阻抗識別巖性具有多解性等缺陷，提供一種疊后地震波阻抗反演方法。
[0006]本發(fā)明包括以下步驟:
[0007]1、獲取待反演工區(qū)的地震處理成果疊后純波數(shù)據(jù)、地質(zhì)分層數(shù)據(jù)、測井曲線，提取地震子波；
[0008]1.1測井曲線標(biāo)準(zhǔn)化處理:依據(jù)該工區(qū)標(biāo)志層測井曲線特征響應(yīng)值對該工區(qū)所有鉆井的測井曲線進(jìn)行標(biāo)準(zhǔn)化編輯處理。
[0009]1.2提取地震子波:利用井旁地震處理成果疊后純波數(shù)據(jù)和步驟1.1得到的測井資料提取地震子波。
[0010]2、應(yīng)用步驟1.2提取的地震子波，生成待反演工區(qū)所有井的合成地震記錄，建立待反演工區(qū)內(nèi)所有井的井-震時深關(guān)系，標(biāo)定各巖性界面在地震剖面上的反射位置，生成所有井的時間域的波阻抗曲線。
[0011]3、對步驟2生成的所有井的時間域的波阻抗曲線進(jìn)行去壓實(shí)校正:對步驟2中得到的時間域波阻抗曲線進(jìn)行濾波處理，去除波阻抗曲線中反映壓實(shí)效應(yīng)的0-0.5Hz極低頻段成分。
[0012]4、根據(jù)步驟1獲取的地質(zhì)分層數(shù)據(jù)和步驟2建立的待反演工區(qū)內(nèi)所有井的井-震時深關(guān)系，在地震處理成果疊后純波數(shù)據(jù)體上，對待反演工區(qū)內(nèi)所有目標(biāo)層的層位進(jìn)行精細(xì)追蹤解釋，建立待反演工區(qū)的地質(zhì)構(gòu)造框架模型。
[0013]5、在步驟4中建立的待反演工區(qū)的地質(zhì)構(gòu)造框架模型約束下，利用步驟3獲取的經(jīng)去壓實(shí)校正后的所有井的時間域的波阻抗曲線，建立待反演工區(qū)的地震波阻抗反演初始模型。
[0014]6、對步驟5建立的待反演工區(qū)的地震波阻抗反演初始模型進(jìn)行疊后波阻抗反演，得到去壓實(shí)波阻抗反演結(jié)果:
[0015]6.1利用步驟1.2獲取的地震子波，對步驟5建立的待反演工區(qū)的地震波阻抗反演初始模型進(jìn)行正演，形成合成地震道數(shù)據(jù)體；
[0016]6.2將步驟6.1形成的合成地震道數(shù)據(jù)體與步驟1中獲取的待反演工區(qū)的地震處理成果疊后純波數(shù)據(jù)進(jìn)行比較產(chǎn)生剩余誤差道；
[0017]6.3采用模型優(yōu)選迭代撓動算法，對步驟5建立的待反演工區(qū)的地震波阻抗反演初始模型進(jìn)行變換、正演，使其與步驟1中獲取的待反演工區(qū)的地震處理成果疊后純波數(shù)據(jù)進(jìn)行比較產(chǎn)生的剩余誤差道逐步變小；
[0018]6.4重復(fù)上述步驟6.3，選取對應(yīng)于剩余誤差道誤差最小的變換模型，為疊后波阻抗反演結(jié)果。
[0019]本發(fā)明的有益效果在于對測井波阻抗曲線進(jìn)行壓實(shí)校正，消除了受地層沉積壓實(shí)效應(yīng)影響的極低頻段波阻抗，建立不受地層壓實(shí)影響的反演初始模型，構(gòu)建了去壓實(shí)效應(yīng)的迭后地震波阻抗反演方法，有效解決了地層沉積壓實(shí)作用對地震波阻抗產(chǎn)生的巖性識別多解性影響。經(jīng)反演得到的波阻抗數(shù)據(jù)體縱向上真實(shí)反映了地層沉積韻律變化，橫向上真實(shí)反映了巖性變化，能夠較好地指導(dǎo)油田勘探開發(fā)部署。
【附圖說明】
[0020]圖1是本發(fā)明技術(shù)方案流程框圖。
[0021]圖2是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼測井曲線一致性校正前直方圖。
[0022]圖3是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼測井曲線一致性校正后直方圖。
[0023]圖4是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼波阻抗反演中提取的地震子波。
[0024]圖5是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼波阻抗反演中合成記錄標(biāo)定圖。
[0025]圖6是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼胡99井測井曲線去壓實(shí)校正圖。
[0026]圖7是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼波阻抗反演中建立的初始模型。
[0027]圖8是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼波阻抗反演中得到的去壓實(shí)波阻抗剖面圖。
[0028]圖9是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼常規(guī)技術(shù)波阻抗波反演剖面圖。
[0029]圖10是本發(fā)明實(shí)施例中東濮凹陷柳屯次洼胡114-1井位部署圖。
【具體實(shí)施方式】
[0030]下面結(jié)合東濮凹陷柳屯次洼沙三中去壓實(shí)波阻抗反演及儲層預(yù)測實(shí)例和附圖，對本發(fā)明實(shí)施方式做進(jìn)一步詳細(xì)說明。由圖1可知，本發(fā)明具體步驟如下:
[0031]1、獲取東濮凹陷柳屯次洼沙三中的地震處理成果疊后純波數(shù)據(jù)、地質(zhì)分層數(shù)據(jù)、測井曲線。
[0032]1.1測井曲線標(biāo)準(zhǔn)化處理:利用直方圖法聲波曲線標(biāo)準(zhǔn)化編輯處理，分兩步進(jìn)行，首先選取可對比追蹤的標(biāo)準(zhǔn)層，選取條件為巖性、電性特征明顯、沉積穩(wěn)定、具有一定厚度的單層或?qū)咏M、工區(qū)內(nèi)多數(shù)井鉆遇、且靠近或在反演的目的層段上的泥巖層段；其次是利用交會分析技術(shù)作出所有井標(biāo)準(zhǔn)層總的測井響應(yīng)頻率直方圖和單井的標(biāo)準(zhǔn)層測井響應(yīng)頻率直方圖。如圖2為東濮凹陷柳屯次洼測井曲線一致性校正前直方圖，從圖2上讀出特征峰值所對應(yīng)的曲線值，以此為標(biāo)準(zhǔn)對測井曲線進(jìn)行標(biāo)準(zhǔn)化處理，以此得到圖3所示東濮凹陷柳屯次洼測井曲線一致性校正后直方圖。
[0033]1.2提取地震子波:首先利用步驟1獲取的東濮凹陷柳屯次洼的地震處理成果疊后純波數(shù)據(jù)和步驟1.1獲取的測井曲線提取地震子波的振幅譜和相位譜，然后利用振幅譜和相位譜的信息合成如圖4所示的東濮凹陷柳屯次洼波阻抗反演中提取的地震子波。上述子波的長度根據(jù)資料的信噪比、頻率特征等因素，通過反復(fù)試驗(yàn)來確定。依據(jù)圖4東濮凹陷柳屯次洼波阻抗反演中提取的地震子波長度一般為100ms左右，計算時窗至少是地震子波長度的3倍以上。
[0034]2、建立井-震時深關(guān)系，生成時間域波阻抗曲線。
[0035]2.1建立東濮凹陷柳屯次洼工區(qū)內(nèi)所有井的井-震時深關(guān)系。用步驟1.2提取的地震子波和步驟1.1獲取標(biāo)準(zhǔn)化測井曲線制作如圖5所示的東濮凹陷柳屯次洼波阻抗反演中合成記錄標(biāo)定圖