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

【發(fā)明內容】

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