結(jié)果。
【具體實(shí)施方式】
[0040] 下面結(jié)合附圖和實(shí)例對(duì)本發(fā)明進(jìn)行詳細(xì)的描述。
[0041] 本發(fā)明從野外采集得到的地震記錄出發(fā)，經(jīng)過預(yù)處理后變?yōu)榀B后地震數(shù)據(jù)或偏移 后疊加地震數(shù)據(jù)，根據(jù)本發(fā)明方法需從剖面解釋和單道處理兩方面出發(fā)對(duì)地震數(shù)據(jù)進(jìn)行處 理。
[0042] 圖1示出本發(fā)明基于地震數(shù)據(jù)瞬時(shí)屬性的強(qiáng)振幅屏蔽層剔除方法的流程，包括以 下步驟：
[0043] 1)在疊后或偏移后疊加后數(shù)據(jù)上對(duì)強(qiáng)振幅屏蔽層同相軸進(jìn)行精細(xì)解釋；
[0044] 2)在疊后或偏移后疊加后數(shù)據(jù)上以步驟1)解釋好的屏蔽軸為中心，設(shè)計(jì)時(shí)窗，截 取時(shí)窗內(nèi)的地震數(shù)據(jù)；
[0045] 3)通過步驟1)解釋好的強(qiáng)振幅屏蔽層，獲取強(qiáng)振幅屏蔽層的時(shí)間和振幅信息，對(duì) 步驟2)截取的疊后或偏移疊加后的地震數(shù)據(jù)進(jìn)行Hilbert變換，逐道獲取構(gòu)成強(qiáng)振幅屏蔽 層子波的頻率和相位，通過精細(xì)井震標(biāo)定工作提取平均子波的寬度，在Morlet子波的基礎(chǔ) 上逐道構(gòu)建近似于強(qiáng)振幅屏蔽層的子波
[0046] Morlet小波是一種單頻復(fù)正弦調(diào)制高斯波，也是最常用的復(fù)值小波，其時(shí)頻域都 具有很好的局部性。Morlet小波的構(gòu)成取決于振幅、相位、頻率、中心時(shí)移、尺度參數(shù)，構(gòu)造 近似強(qiáng)振幅屏蔽軸的Morlet小波，具體步驟包括：
[0047] ①對(duì)通過步驟2)截取的疊后或偏移后疊加地震數(shù)據(jù)，逐道進(jìn)行Hilbert變換，對(duì) 第j道地震數(shù)據(jù)進(jìn)行Hilbert變換，構(gòu)建復(fù)地震道 Sj (t)
[0049] 式中，Xj⑴為第j道地震數(shù)據(jù)（j = 1，2. . . N)，巧〇：)為Xj⑴的Hilbert變換結(jié) 果；
[0050] ②逐道獲取步驟1)中解釋好的強(qiáng)振幅屏蔽軸的時(shí)間和振幅值，第j道地震數(shù)據(jù)振 幅值a,。作為構(gòu)建強(qiáng)振幅屏蔽層子波的振幅，以u(píng) ,。作為構(gòu)建強(qiáng)振幅屏蔽層子波的中心時(shí)移 時(shí)間；由于強(qiáng)振幅屏蔽層在地震疊后或偏移疊加后剖面上能量強(qiáng)、易拾取，因此若能利用強(qiáng) 振幅屏蔽層的解釋結(jié)果則能夠避免多子波分解所要進(jìn)行的復(fù)雜運(yùn)算
[0051] ③通過S] (t)，獲取第j道地震信號(hào)的瞬時(shí)頻率和瞬時(shí)相位，并記錄下U]。處的瞬時(shí) 頻率和瞬時(shí)相位

[0054] 式中0j〇為瞬時(shí)相位，U j。為瞬時(shí)頻率；
[0055] ④結(jié)合測(cè)井資料，進(jìn)行井震標(biāo)定，并提取平均子波，將平均子波與不同尺度的 Morlet小波進(jìn)行互相關(guān)，并進(jìn)一步確定〇。；
[0056] ⑤利用獲取中心時(shí)移時(shí)間Uj。以及u j。位置處所對(duì)應(yīng)的瞬時(shí)頻率ω j。和瞬時(shí)相位 0j〇，構(gòu)建Morlet子波，Morlet小波的表達(dá)式為
[0058] ⑥將步驟3)獲取的振幅值a,。與第⑤步得到的Morlet小波相乘，獲取第j道的近 似強(qiáng)振幅屏蔽軸子波
[0059] P j (t) = aj0Xmj (t)
[0060] 4)將疊加或偏移疊加后數(shù)據(jù)逐道與第3)步獲取的Morlet子波做差，將差值D逐 道輸出作為剔除強(qiáng)振幅屏蔽層后的結(jié)果
[0062] 圖2給出了不同尺度參數(shù)所對(duì)應(yīng)的Morlet子波的形態(tài)，用以說明尺度參數(shù)和 Morlet 子波的關(guān)系，亮藍(lán)色 σ e [0· 1，0.5];藍(lán)色 σ e [0.6，L0];紅色 σ e [1. 1，1.5]; 綠色 σ e [1. 6, 2. 0];粉色 σ e [2. 1，2. 5];黑色 σ e [2. 6, 3. 0];可見尺度參數(shù)對(duì)Morlet 小波形態(tài)的影響主要在尺度參數(shù)為I. 〇、l. 5、2. 0、2. 5、3. 0附近，尺度參數(shù)每增加0. 5,其 Morlet子波會(huì)增加一個(gè)明顯的旁瓣；
[0063] 圖3給出了井震標(biāo)定后提取的平均子波；
[0064] 圖4給出了圖3與圖2中每個(gè)子波進(jìn)行互相關(guān)后，其相關(guān)系數(shù)與尺度參數(shù)的關(guān)系 圖；
[0065] 圖5給出了經(jīng)過精細(xì)解釋后的實(shí)際地震剖面，解釋層位所在位置即為強(qiáng)振幅屏蔽 層的位置；
[0066] 圖6給出了剔除強(qiáng)振幅屏蔽層后的地震剖面，對(duì)比圖5和圖6可見，剔除強(qiáng)振幅 屏蔽層后的地震剖面，原強(qiáng)振幅屏蔽層所在位置下方的同相軸能量較原始地震剖面有所增 強(qiáng)，且地震剖面中同相軸形態(tài)與原始地震資料相比并無異常。
[0067] 圖7給出了基于圖5已解釋層位的基礎(chǔ)上所重構(gòu)的強(qiáng)振幅屏蔽層，可見子波的能 量隨原始實(shí)際地震剖面上強(qiáng)振幅屏蔽層的能量而變化，所重構(gòu)的強(qiáng)振幅屏蔽層子波能夠與 原始實(shí)際地震剖面中強(qiáng)振幅屏蔽層實(shí)現(xiàn)較好的匹配；
[0068] 圖8給出了圖5、圖6中標(biāo)注位置處的單道地震記錄，圖8 (a)是原始地震剖面中所 標(biāo)注位置處的單道記錄；(b)是剔除強(qiáng)振幅屏蔽層后的實(shí)際地震剖面中所標(biāo)注位置處的單 道記錄，對(duì)比圖8(a)和（b)可見，在3200ms附近，本發(fā)明方法可以有效削弱強(qiáng)振幅屏蔽層 的能量，突顯強(qiáng)振幅屏蔽層下方有效弱信號(hào)的能量，且對(duì)地震記錄其余位置處的能量并無 影響；
[0069] 圖9給出了圖8 (a)、圖8 (b)的時(shí)頻譜，圖9 (a)是圖8 (a)所對(duì)應(yīng)的時(shí)頻譜；圖9 (b) 是圖8(b)所對(duì)應(yīng)的時(shí)頻譜，對(duì)比圖9(a)和圖9(b)可見剔除強(qiáng)振幅屏蔽層后的地震道，在 原有強(qiáng)振幅屏蔽層位置下方（3200ms附近）的能量有所增強(qiáng)，在3200ms附近原強(qiáng)振幅屏蔽 層位置處頻率分辨率有所提高。
【主權(quán)項(xiàng)】
1. 一種基于地震數(shù)據(jù)瞬時(shí)屬性的強(qiáng)振幅屏蔽層剔除方法，所述方法其特征在于，包括 以下步驟： 1) 在疊后或偏移后疊加后數(shù)據(jù)上對(duì)強(qiáng)振幅屏蔽層同相軸進(jìn)行精細(xì)解釋； 2) 在疊后或偏移后疊加后數(shù)據(jù)上以步驟1)解釋好的屏蔽軸為中心，設(shè)計(jì)時(shí)窗，截取時(shí) 窗內(nèi)的地震數(shù)據(jù)； 3) 在第1)步解釋好的強(qiáng)振幅屏蔽層的基礎(chǔ)上，獲取強(qiáng)振幅屏蔽層的中心時(shí)間信息 和振幅信息，對(duì)第2)步截取的疊后或偏移疊加后的地震數(shù)據(jù)進(jìn)行Hilbert變換，逐道獲 取構(gòu)成強(qiáng)振幅屏蔽層子波的頻率和相位，通過精細(xì)井震標(biāo)定工作提取平均子波的寬度，在 Morlet子波的基礎(chǔ)上逐道構(gòu)建近似于強(qiáng)振幅屏蔽層的子波； 4) 將疊后或偏移疊加后的地震數(shù)據(jù)逐道與第3)步獲取的Morlet子波做差，將差值逐 道輸出作為剔除強(qiáng)振幅屏蔽層后的結(jié)果。2. 根據(jù)權(quán)利要求1所述的基于地震數(shù)據(jù)瞬時(shí)屬性的強(qiáng)振幅屏蔽層剔除方法，所述步驟 3) 中，構(gòu)造近似強(qiáng)振幅屏蔽軸的Morlet小波，具體步驟包括： ① 對(duì)通過步驟2)截取的疊加或偏移后疊加的地震數(shù)據(jù)逐道進(jìn)行Hilbert變換，對(duì)第j 道地震數(shù)據(jù)進(jìn)行Hilbert變換，構(gòu)建復(fù)地震道S] (t)式中，Xi⑴為第j道地震數(shù)據(jù)（j=1，2···Ν)，Xj(t)丨為Xj⑴的Hilbert變換結(jié)果； ② 逐道獲取步驟1)中解釋好的強(qiáng)振幅屏蔽軸的時(shí)間和振幅值，第j道地震數(shù)據(jù)振幅值 aj(]作為構(gòu)建強(qiáng)振幅屏蔽層子波的振幅，uj。作為構(gòu)建強(qiáng)振幅屏蔽層子波的中心時(shí)移時(shí)間； ③ 通過S]⑴，獲取第j道地震信號(hào)的瞬時(shí)頻率和瞬時(shí)相位，并記錄下U]。處的瞬時(shí)頻率 和瞬時(shí)相位式中0j〇為瞬時(shí)相位，《i。為瞬時(shí)頻率； ④ 結(jié)合測(cè)井資料，進(jìn)行井震標(biāo)定，并提取平均子波，根據(jù)平均子波與不同尺度Morlet 小波的相關(guān)系數(shù)，進(jìn)一步確定σ。； ⑤ 利用獲取中心時(shí)移時(shí)間+。以及uj。位置處對(duì)應(yīng)的瞬時(shí)頻率ωj。和瞬時(shí)相位0jo,構(gòu) 建Morlet子波，Morlet小波的表達(dá)式為⑥ 將步驟3)獲取的振幅值a]。與第⑤步得到的Morlet小波相乘，獲取第j道的近似強(qiáng) 振幅屏蔽軸子波 p.j(t) =a^Xm^t) 〇3. 根據(jù)權(quán)利要求1所述的基于地震數(shù)據(jù)瞬時(shí)屬性的強(qiáng)振幅屏蔽層剔除方法，所述步驟 4) 中，獲取剔除強(qiáng)振幅屏蔽層后的數(shù)據(jù)D，具體步驟包括：
【專利摘要】本發(fā)明涉及一種基于地震數(shù)據(jù)瞬時(shí)屬性的強(qiáng)振幅屏蔽層剔除方法，其特征在于，包括以下步驟：1)在疊后或偏移后數(shù)據(jù)上對(duì)強(qiáng)振幅屏蔽層進(jìn)行精細(xì)解釋；2)以解釋好的強(qiáng)振幅屏蔽層為中心，設(shè)計(jì)時(shí)窗，截取時(shí)窗內(nèi)的地震數(shù)據(jù)；3)獲取已解釋好的強(qiáng)振幅屏蔽層的時(shí)間和振幅信息，對(duì)時(shí)窗內(nèi)的地震數(shù)據(jù)進(jìn)行Hilbert變換，獲取構(gòu)成強(qiáng)振幅屏蔽層子波的頻率和相位，通過井震標(biāo)定提取平均子波，根據(jù)平均子波與不同尺度參數(shù)理論子波的相關(guān)系數(shù)，確定平均子波的尺度參數(shù)，逐道構(gòu)建近似于強(qiáng)振幅屏蔽層的Morlet子波；4)將第1)步的數(shù)據(jù)體與第3)步獲取的Morlet子波逐道做差，得到剔除強(qiáng)振幅屏蔽層后的地震記錄。本發(fā)明處理后的數(shù)據(jù)能夠有效消除地震剖面上強(qiáng)振幅屏蔽層的影響。
【IPC分類】G01V1/40, G01V1/28
【公開號(hào)】CN105259579
【申請(qǐng)?zhí)枴緾N201510758326
【發(fā)明人】王淑艷, 楊勇, 劉磊, 王軍, 王彥春, 劉學(xué)清, 楊博文, 劉鑫, 班麗, 李靚, 王瑋, 李忠新, 孫渡, 李超
【申請(qǐng)人】中國(guó)地質(zhì)大學(xué)（北京）, 中國(guó)石油化工股份有限公司勝利油田分公司勘探開發(fā)研究院
【公開日】2016年1月20日
【申請(qǐng)日】2015年11月10日