亚洲狠狠干,亚洲国产福利精品一区二区,国产八区,激情文学亚洲色图

一種基于地震數(shù)據(jù)瞬時(shí)屬性的強(qiáng)振幅屏蔽層剔除方法_2

文檔序號(hào):9505451閱讀:來源:國(guó)知局
結(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日
當(dāng)前第2頁(yè)1 2 
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評(píng)論。精彩留言會(huì)獲得點(diǎn)贊!
1