[發(fā)明專利]一種海上自適應(yīng)壓制鬼波的寬頻逆時偏移成像方法

有效

• 申請?zhí)枺?/td>: 201610949368.1
• 申請日: 2016-10-26
• 公開（公告）號: CN106526677B
• 公開（公告）日: 2018-12-21
• 發(fā)明（設(shè)計）人: 孫文博;張金淼;朱振宇;翁斌;王艷冬;郝振江;王小六;糜芳;薛東川;黃小剛
• 申請（專利權(quán)）人: 中海石油（中國）有限公司;中海石油（中國）有限公司北京研究中心
• 主分類號: G01V1/38
• 分類號: G01V1/38;G01V1/36
• 代理公司: 北京紀凱知識產(chǎn)權(quán)代理有限公司 11245
• 代理人: 關(guān)暢;劉美麗

本發(fā)明涉及一種海上自適應(yīng)壓制鬼波的寬頻逆時偏移成像方法，其特征在于包括以下內(nèi)容：1)去噪后的疊前炮集地震記錄；2)將疊前炮集地震記錄進行時間域速度分析，并進行深度域速度建模，得到深度域速度體；對疊前炮集地震記錄進行數(shù)據(jù)解析，拆分為數(shù)據(jù)和道頭兩部分信息，并利用道頭信息生成炮集數(shù)據(jù)索引文件；3)讀取炮集數(shù)據(jù)，將讀取的炮集數(shù)據(jù)作為寬頻逆時偏移輸入數(shù)據(jù)體，并將深度域速度體作為寬頻逆時偏移輸入?yún)?shù)，進行三維鬼波自適應(yīng)壓制，完成單炮寬頻逆時偏移成像得到單炮集偏移結(jié)果；4)對單炮集偏移結(jié)果進行噪音壓制處理，并根據(jù)單炮偏移結(jié)果的實際空間位置坐標，對所有噪音壓制后的單炮偏移結(jié)果進行疊加，獲得最終的寬頻成像結(jié)果。

技術(shù)領(lǐng)域

本發(fā)明涉及一種海上自適應(yīng)壓制鬼波的寬頻逆時偏移成像方法，涉及地震成像技術(shù)領(lǐng)域。

背景技術(shù)

海上地震資料采集時，氣槍和檢波器沉放于水面以下一定深度，因而地震記錄中必然出現(xiàn)由這種自由表面產(chǎn)生的虛反射干擾——鬼波，鬼波主要包括震源鬼波、檢波點鬼波和震源-檢波點混合鬼波。由于鬼波的存在導(dǎo)致地震記錄中出現(xiàn)嚴重的陷頻問題，影響高頻和低頻成分，鬼波特征與地震記錄中一次反射波相似，消除比較困難，嚴重干擾有效波，直接影響地震數(shù)據(jù)后續(xù)的反演和解釋。

近幾年鬼波壓制在海上地震記錄處理中非常重要，在采集、處理等方面研發(fā)了一系列壓制鬼波的寬頻地震技術(shù)，也成為業(yè)界研究的熱點問題?，F(xiàn)有技術(shù)的采集技術(shù)主要有：變深度纜、上下纜、雙檢等；處理方法為反褶積類、t-p變換類、基于Green理論的鬼波壓制方法等，鬼波是隨出射角發(fā)生變化，即與偏移距相關(guān)的三維算子，因而需要三維算法才能實現(xiàn)鬼波有效壓制。常規(guī)鬼波壓制處理方法都是基于地震記錄的疊前道集數(shù)據(jù)，另外實現(xiàn)全三維鬼波壓制相對困難，壓制效果也會受到一定影響。

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

針對上述問題，本發(fā)明的目的是提供一種能夠?qū)崿F(xiàn)三維鬼波壓制的海上自適應(yīng)壓制鬼波的寬頻逆時偏移成像方法。

為實現(xiàn)上述目的，本發(fā)明采取以下技術(shù)方案：一種海上自適應(yīng)壓制鬼波的寬頻逆時偏移成像方法，包括以下內(nèi)容：1)對地震記錄進行疊前噪音壓制處理，獲得去噪后的疊前炮集地震記錄；2)將疊前炮集地震記錄進行時間域速度分析，并進行深度域速度建模，得到深度域速度體；對疊前炮集地震記錄進行數(shù)據(jù)解析，拆分為數(shù)據(jù)和道頭兩部分信息，并利用道頭信息生成炮集數(shù)據(jù)索引文件；3)根據(jù)炮集數(shù)據(jù)索引文件記錄的數(shù)據(jù)索引關(guān)系，并行讀取炮集數(shù)據(jù)，將讀取的炮集數(shù)據(jù)作為寬頻逆時偏移輸入數(shù)據(jù)體，并將深度域速度體作為寬頻逆時偏移輸入?yún)?shù)，進行三維鬼波的自適應(yīng)壓制，完成單炮寬頻逆時偏移成像得到單炮集偏移結(jié)果；4)對每個單炮集偏移結(jié)果進行噪音壓制處理，并根據(jù)每個單炮偏移結(jié)果的實際空間位置坐標，對所有噪音壓制后的單炮偏移結(jié)果進行疊加，獲得最終的寬頻成像結(jié)果。

進一步地，所述步驟2)得到深度域速度體的具體過程為：首先，將得到的疊前炮集地震記錄抽成共中心點道集數(shù)據(jù)，并對共中心點道集數(shù)據(jù)，利用時間域速度分析方法進行多輪速度分析迭代，得到時間域速度體；然后，以時間域速度體為基礎(chǔ)，采用深度建模方法進行多次迭代，得到深度域速度體。

進一步地，數(shù)據(jù)索引文件用于記錄每炮信息情況，包括每炮道數(shù)和采樣點數(shù)。

進一步地，所述步驟3)的具體計算過程為：根據(jù)地震波傳播理論，鬼波表達式為：

式中，θ為鬼波初射角，ω為角頻率，z為震源或檢波點的沉放深度，v是海水速度，f是鬼波陷頻點，將含鬼波新邊值條件應(yīng)用于逆時偏移波場延拓所用波動方程中，以實現(xiàn)成像過程中壓制鬼波，具體形式如下：

a)實際檢測點為某個高度，根據(jù)逆時偏移震源波場正演所用新的修改邊值條件公式(1)，利用時間方向為四階差分精度的偽譜法計算正演波場P_F：