一種基于疊前瞬時(shí)頻率屬性分析的儲(chǔ)層烴類預(yù)測(cè)方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種石油地球物理勘探技術(shù)領(lǐng)域，具體地說是一種基于疊前瞬時(shí)頻率 屬性分析的儲(chǔ)層烴類預(yù)測(cè)方法。
【背景技術(shù)】
[0002] 地震資料解釋是地震數(shù)據(jù)處理的重要環(huán)節(jié)，它是將地震資料轉(zhuǎn)化成我們對(duì)勘探區(qū) 地下地質(zhì)情況的認(rèn)識(shí)，并從中找到與油氣有關(guān)的信息。在相當(dāng)長(zhǎng)的一段時(shí)間內(nèi)，振幅屬性以 及AVO特征都是地震資料解釋的重要信息，然而在實(shí)際地層中，由于上覆下浮地層介質(zhì)參 數(shù)的不同，含氣砂巖所表現(xiàn)出來的AVO異常也不盡相同，除了亮點(diǎn)異常外，還有暗點(diǎn)、相位 反轉(zhuǎn)等特征，從而導(dǎo)致了基于振幅屬性的資料信息在實(shí)際應(yīng)用中的多解性問題。對(duì)于某個(gè) 特定地層，需要有針對(duì)性的選擇能夠反映地質(zhì)現(xiàn)象或油氣信息的屬性加以分析，才能更準(zhǔn) 確的預(yù)測(cè)儲(chǔ)層、揭示地下異常。
[0003] 在疊前地震信息中，瞬時(shí)頻率獨(dú)立于振幅屬性，有著明確的物理意義，地層的吸收 作用使得反射波能量發(fā)生衰減，而高頻成分衰減速度大于低頻成分，使得主頻逐漸向低頻 移動(dòng)，特別是當(dāng)儲(chǔ)層空隙中含有烴類物質(zhì)時(shí)，品質(zhì)因素 Q會(huì)變大，從而導(dǎo)致地震反射波主頻 有更明顯的降低，這種頻率在空間的變化也已成為指示油氣藏存在的重要指標(biāo)。
[0004] Hilbert變換，也稱希爾伯特變換，是以著名數(shù)學(xué)家大衛(wèi)?希爾伯特（David Hilbert)來命名。希爾伯特變換（hilbert transform) -個(gè)連續(xù)時(shí)間信號(hào)x (t)的希爾伯 特變換等于該信號(hào)通過具有沖激響應(yīng)h(t) =1/Jit的線性系統(tǒng)以后的輸出響應(yīng)xh(t)。所 以信號(hào)經(jīng)希爾伯特變換后，在頻域各頻率分量的幅度保持不變，但相位將出現(xiàn)90°相移。即 對(duì)正頻率滯后π/2,對(duì)負(fù)頻率導(dǎo)前π/2,因此希爾伯特變換器又稱為90°移相器?，F(xiàn)在的 復(fù)地震道分析技術(shù)主要是建立在Hilbert變換基礎(chǔ)上，將地震波看作是一種解析信號(hào)，通 過Hilbert變換，獲得解析信號(hào)的虛部，進(jìn)而求得復(fù)地震道屬性，根據(jù)屬性特點(diǎn)判斷儲(chǔ)層巖 性變化，預(yù)測(cè)油氣。該方法在相當(dāng)長(zhǎng)的一段時(shí)間內(nèi)，被廣泛應(yīng)用于疊后地震資料中。如何基 于疊前瞬時(shí)頻率屬性進(jìn)行分析，以得到儲(chǔ)層烴類預(yù)測(cè)，是目前還空白的技術(shù)。

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

[0005] 本發(fā)明的技術(shù)任務(wù)是提供一種減小了解釋誤差，提高了烴類預(yù)測(cè)精度的一種基于 疊前瞬時(shí)頻率屬性分析的儲(chǔ)層烴類預(yù)測(cè)方法。
[0006] 本發(fā)明的技術(shù)任務(wù)是按以下方式實(shí)現(xiàn)的：一種基于疊前瞬時(shí)頻率屬性分析的儲(chǔ)層 烴類預(yù)測(cè)方法，先對(duì)疊前地震信號(hào)做加窗處理，消除干擾波對(duì)有效信號(hào)的影響；再利用復(fù) 地震道技術(shù)，通過Hilbert變換獲取地震信號(hào)的瞬時(shí)頻率屬性及其隨偏移距的變化信息， 通過地震信號(hào)的瞬時(shí)頻率屬性及其隨偏移距的變化信息得到FVO曲線（FVO-Frequency Versus Offset頻率隨偏移距）；將每條FVO曲線順時(shí)針旋轉(zhuǎn)90°，然后將旋轉(zhuǎn)后多個(gè)道集 的FVO曲線按順序排列在一起，掛在一個(gè)目的儲(chǔ)藏層位上，得到基于瞬時(shí)頻率屬性的疊前 偏移距體；結(jié)合模型和測(cè)井資料，找到與油氣相關(guān)的變化特征；在基于瞬時(shí)頻率的疊前偏 移距體中，烴類異常區(qū)的瞬時(shí)頻率值明顯小于同一偏移距非異常區(qū)的瞬時(shí)頻率，而且隨著 偏移距的增加而減小；從而進(jìn)行儲(chǔ)層烴類預(yù)測(cè)。
[0007] 實(shí)現(xiàn)步驟為： (1) 、采用人工激發(fā)并由檢波器記錄地震波，取得原始疊前地震數(shù)據(jù)，去噪保幅處理，抽 取疊前地震道集并對(duì)其做同相軸拉平處理； (2) 、通過Hilbert變換獲取復(fù)地震道信號(hào)X(t) = f(t)+ih(t)，繼而求得瞬時(shí)頻率
【主權(quán)項(xiàng)】
1. 一種基于疊前瞬時(shí)頻率屬性分析的儲(chǔ)層姪類預(yù)測(cè)方法，其特征在于先對(duì)疊前地震信 號(hào)做加窗處理，消除干擾波對(duì)有效信號(hào)的影響；再利用復(fù)地震道技術(shù)，通過化化ed變換獲 取地震信號(hào)的瞬時(shí)頻率屬性及其隨偏移距的變化信息，通過地震信號(hào)的瞬時(shí)頻率屬性及其 隨偏移距的變化信息得到FVO曲線；將每條FVO曲線順時(shí)針旋轉(zhuǎn)90°，然后將旋轉(zhuǎn)后多個(gè) 道集的FVO曲線按順序排列在一起，掛在一個(gè)目的儲(chǔ)藏層位上，得到基于瞬時(shí)頻率屬性的 疊前偏移距體；結(jié)合模型和測(cè)井資料，找到與油氣相關(guān)的變化特征；在基于瞬時(shí)頻率的疊 前偏移距體中，姪類異常區(qū)的瞬時(shí)頻率值明顯小于同一偏移距非異常區(qū)的瞬時(shí)頻率，而且 隨著偏移距的增加而減??；從而進(jìn)行儲(chǔ)層姪類預(yù)測(cè)。
2. 根據(jù)權(quán)利要求1所述的一種基于疊前瞬時(shí)頻率屬性分析的儲(chǔ)層姪類預(yù)測(cè)方法，其特 征在于實(shí)現(xiàn)步驟為： (1) 、采用人工激發(fā)并由檢波器記錄地震波，取得原始疊前地震數(shù)據(jù)，去噪保幅處理，抽 取疊前地震道集并對(duì)其做同相軸拉平處理； (2) 、通過化化6的變換獲取復(fù)地震道信號(hào)X(t) = f(t)+化(t)，繼而求得瞬時(shí)頻率 W似，拾取每個(gè)CMP道集的FVO曲線； (it (3) 、將每條FVO曲線順時(shí)針旋轉(zhuǎn)90°，此時(shí)數(shù)值方向是每個(gè)CMP道集的偏移距，水平方 向是偏移距對(duì)應(yīng)的瞬時(shí)頻率信息； (4) 、將旋轉(zhuǎn)后的多個(gè)CMP道集的FV0曲線順序排列在一起，掛在目的儲(chǔ)層上，由此得到 基于瞬時(shí)頻率屬性的疊前偏移距體，縱坐標(biāo)表示偏移距信息，橫坐標(biāo)表示CMP道號(hào)； 巧）、瞬時(shí)頻率的突變位置對(duì)應(yīng)了含姪類異常儲(chǔ)層的位置，當(dāng)儲(chǔ)層中含有某種姪類物質(zhì) 時(shí)，地層衰減因子變大，從而使反射波主頻有更明顯的降低，利用該種頻率在空間的變化可 指示姪類異常。
3. 根據(jù)權(quán)利要求1所述的一種基于疊前瞬時(shí)頻率屬性分析的儲(chǔ)層姪類預(yù)測(cè)方法，其特 征在于在基于高密觀測(cè)系統(tǒng)的S層地層模型疊前資料解釋上的應(yīng)用，具體步驟如下： (1) 、在基于高密觀測(cè)系統(tǒng)的S層地層模型中有四個(gè)區(qū)域，區(qū)域1縱波速度化=4600m/ S，橫波速度Vs =化/2,密度P = 2. 2g/cm3,品質(zhì)因數(shù)Q = 100 ;區(qū)域2中化=4500m/s， S =化/2,密度P = 2. 2g/cm3,品質(zhì)因數(shù)Q = 80 ;區(qū)域3中化=4200m/s，S =化/2,密度 P = 2. Og/cm3,品質(zhì)因數(shù)Q = 20,區(qū)域4中介質(zhì)參數(shù)與區(qū)域1相同； (2) 、該模型在2500m深度的中間位置有個(gè)長(zhǎng)600m，厚50m的姪類異常區(qū)，即區(qū)域3 ;在 地表設(shè)置了 91個(gè)激發(fā)炮點(diǎn)，炮間隔為50m ;每一炮對(duì)應(yīng)1000個(gè)檢波器，間距為5m ; (3) 、該模型的部分疊前地震記錄，共有91000道，采樣點(diǎn)數(shù)為3400,采樣間隔為500US， 共1900個(gè)CMP道集； (4) 、任取一道目的層地震信號(hào)，對(duì)目的層加窗，并通過HUbed變換得到有用信號(hào)的 瞬時(shí)頻率； 巧）、沿目的層同相軸拾取異常區(qū)CMP的頻率隨偏移距的變化曲線，即FV0曲線； 化）、縱坐標(biāo)為瞬時(shí)頻率信息，橫坐標(biāo)為偏移距；瞬時(shí)頻率在零偏位置值最大，隨著偏移 距的增加，逐漸減??； (7)、分析道集的地震記錄，拾取每個(gè)CMP道集的FV0曲線，經(jīng)變換，得到道集在正偏移 距范圍內(nèi)的基于瞬時(shí)頻率屬性的疊前偏移距體； 巧）、在該基于瞬時(shí)頻率屬性的疊前偏移距體中，瞬時(shí)頻率發(fā)生異常突變的區(qū)域，瞬時(shí) 頻率的突變位置對(duì)應(yīng)了含姪類異常儲(chǔ)層的位置，當(dāng)儲(chǔ)層中含有某種姪類物質(zhì)時(shí)，地層衰減 因子變大，從而使反射波主頻有更明顯的降低，利用該種頻率在空間的變化可指示姪類異 常。
【專利摘要】本發(fā)明公開了一種基于疊前瞬時(shí)頻率屬性分析的儲(chǔ)層烴類預(yù)測(cè)方法，屬于石油地球物理勘探技術(shù)領(lǐng)域，技術(shù)方案為：先對(duì)疊前地震信號(hào)做加窗處理；再利用復(fù)地震道技術(shù)，通過Hilbert變換獲取地震信號(hào)的瞬時(shí)頻率屬性及其隨偏移距的變化信息，將每條FVO曲線順時(shí)針旋轉(zhuǎn)90°，然后將旋轉(zhuǎn)后多個(gè)道集的FVO曲線按順序排列在一起，掛在一個(gè)目的儲(chǔ)藏層位上，得到基于瞬時(shí)頻率屬性的疊前偏移距體；結(jié)合模型和測(cè)井資料，找到與油氣相關(guān)的變化特征；在基于瞬時(shí)頻率的疊前偏移距體中，烴類異常區(qū)的瞬時(shí)頻率值明顯小于同一偏移距非異常區(qū)的瞬時(shí)頻率，而且隨著偏移距的增加而減??；從而進(jìn)行儲(chǔ)層烴類預(yù)測(cè)。本發(fā)明減小了解釋誤差，提高了烴類預(yù)測(cè)精度。
【IPC分類】G01V1-30
【公開號(hào)】CN104597497
【申請(qǐng)?zhí)枴緾N201510088511
【發(fā)明人】高亞力
【申請(qǐng)人】浪潮電子信息產(chǎn)業(yè)股份有限公司
【公開日】2015年5月6日
【申請(qǐng)日】2015年2月26日