基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法
【專利摘要】本發(fā)明涉及一種提高地震記錄分辨率的同時(shí)，壓制噪聲，以恢復(fù)反射系數(shù)為目標(biāo)的反褶積方法，尤其是一種基于時(shí)變子波的分頻段處理的反褶積方法，以解決常用反褶積方法提高分辨率的同時(shí)卻降低信噪比，和求解不穩(wěn)定的問(wèn)題。該方法包括：采用井控多道統(tǒng)計(jì)自相關(guān)法提取地震記錄中的時(shí)變子波；依據(jù)井資料提取高、低頻算子，用于構(gòu)造對(duì)應(yīng)原始地震記錄的高、低頻信息，作為約束條件；對(duì)分頻信息擴(kuò)散濾波處理；原始地震記錄及其高、低頻信息聯(lián)立反褶積。本發(fā)明使地震資料與測(cè)井資料有機(jī)的結(jié)合起來(lái)，避免了常規(guī)反褶積方法提取的不可控的反褶積因子，壓制噪音的同時(shí)能有效提高分辨率，考慮了地震信號(hào)的衰減，實(shí)現(xiàn)了時(shí)域點(diǎn)譜白化的功能。
【專利說(shuō)明】基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種提高地震記錄分辨率的同時(shí)，壓制噪聲，以恢復(fù)反射系數(shù)為目標(biāo)的反褶積方法，尤其是一種基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法。
【背景技術(shù)】
[0002]提高地震資料分辨率一直是地震資料數(shù)字處理工作者的追求目標(biāo)。高分辨率地震資料數(shù)據(jù)處理的關(guān)鍵環(huán)節(jié)就是壓縮地震子波，或者去除地震波在地下傳播過(guò)程中干涉、調(diào)諧等效應(yīng)對(duì)地下地層的影響，拓寬有效地震信號(hào)的頻帶范圍，特別是較為準(zhǔn)確地拓寬高頻成分。薄層的識(shí)別和厚度成為當(dāng)前地震石油勘探領(lǐng)域的主要研究方向。薄層反射系數(shù)反演方法也就應(yīng)運(yùn)而生。
[0003]譜反演的原理就是根據(jù)時(shí)間域褶積模型，從地震記錄中去除地震子波的影響，進(jìn)而得到反射系數(shù)序列。在時(shí)間域內(nèi)一個(gè)脈沖對(duì)的表達(dá)為:
[0004]δ (t-tD+i^ δ (t-t「T) (1)
[0005]其中Rι為層頂部反射系數(shù)，r2為層底部反射系數(shù)，t為時(shí)間位置，h為頂部反射的時(shí)間位置，T為層厚度。
[0006]將分析點(diǎn)放在層的中心點(diǎn)位置，進(jìn)行傅立葉變換后，用三角法則進(jìn)行化簡(jiǎn)約去t，得到:
[0007]g (f) =2recos ( π fT) +2r0sin ( π fT) (2)
[0008]基于褶積原理，s(t，f)是地震數(shù)據(jù)，w(t，f)是已知的子波。
[0009]
【權(quán)利要求】
1.一種基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，其特征在于，包括以下步驟:51、采用井控下多道統(tǒng)計(jì)自相關(guān)法提取時(shí)變地震子波；52、根據(jù)井資料構(gòu)造高、低頻濾波因子，提取地震記錄的高、低頻信息；53、對(duì)步驟S2各頻段信息進(jìn)行擴(kuò)散濾波，以得到高信噪比和增強(qiáng)同相軸連續(xù)性的高、中、低頻數(shù)據(jù)體；54、對(duì)步驟S3中的高、中、低頻數(shù)據(jù)體進(jìn)行反褶積，以得到提頻剖面。
2.根據(jù)權(quán)利要求1所述基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，其特征在于，在步驟S1中，采用井控下多道統(tǒng)計(jì)自相關(guān)法提取地震記錄中的時(shí)變地震子波的具體步驟為:(101)、提取井旁鄰近道，多道均衡處理；(102)、等間隔取時(shí)窗中心；(103)、各時(shí)窗信號(hào)逐次處理；(104)、各時(shí)窗子波提取結(jié)束，對(duì)子波做能量歸一化處理，第一個(gè)時(shí)窗直接歸一化，其余時(shí)窗，以第一個(gè)時(shí)窗為窗函數(shù)，加窗歸一化；(105)、各時(shí)窗子波光滑 處理。
3.根據(jù)權(quán)利要求2所述基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，其特征在于，在步驟(101)中，將井旁鄰近 道資料與通過(guò)檢測(cè)后得出的井資料進(jìn)行對(duì)比，以確定井旁鄰近道資料中時(shí)變地震子波是否準(zhǔn)確。
4.根據(jù)權(quán)利要求3所述基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，其特征在于，在步驟(103)中，各時(shí)窗信號(hào)逐次處理的具體步驟為:(a)、加窗處理；(b)、求自相關(guān)；(c)、該時(shí)窗信號(hào)處理，保證其中心點(diǎn)位最大振幅值；(d)、帶通濾波；(e)、濾波后的信號(hào)變換到頻率域，去噪處理；(f)、反變換到時(shí)間域，再做一次帶通濾波，得到該時(shí)窗的子波。
5.根據(jù)權(quán)利要求1所述基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，，其特征在于，在步驟S2中，井控提取高、低頻濾波因子Py PH，構(gòu)造高、低頻信息Sp SH，fy,ι?)^ ',Sll=Rl?Sn=PllS,Rll=PllRo(12)
6.根據(jù)權(quán)利要求1所述基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，其特征在于，在步驟S3中，對(duì)各頻段信息做擴(kuò)散濾波，去噪的同時(shí)增強(qiáng)同相軸的連續(xù)性，形成高信噪比和產(chǎn)狀保持良好的高、中、低頻三個(gè)數(shù)據(jù)體。
7.根據(jù)權(quán)利要求1所述基于對(duì)時(shí)變子波進(jìn)行分頻段處理的反褶積方法，，其特征在于，在步驟S4中，在去噪后的高、中、低頻三個(gè)數(shù)據(jù)體進(jìn)行反褶積，利用下式得出提頻剖面，mmJ = p￡ ~WlRl\^ +p-叫|2 +pH ~WhRh\^。(13)
【文檔編號(hào)】G01V1/36GK103713317SQ201210380227
【公開日】2014年4月9日 申請(qǐng)日期:2012年9月29日 優(yōu)先權(quán)日:2012年9月29日
【發(fā)明者】劉力輝 申請(qǐng)人:北京諾克斯達(dá)石油科技有限公司