專利名稱：抑制地表地震數(shù)據(jù)中的殘余水底能量的方法、程序產(chǎn)品和系統(tǒng)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地震數(shù)據(jù)(seismic data)處理領(lǐng)域。具體地，本發(fā)明涉及用于抑制地 震數(shù)據(jù)中不想要的能量，特別是以殘余水底能量(residual water bottom energy)的形式的
能量的系統(tǒng)、程序產(chǎn)品和相關(guān)的方法。
背景技術(shù)：
可以用部署在海床上的接收點陣列在海洋環(huán)境中收集地震反射/折射數(shù)據(jù)。但 是，這樣的數(shù)據(jù)獲取系統(tǒng)不僅記錄可用數(shù)據(jù)，還記錄大量截留在水層中的能量。這種能 量通過緊跟在來自更深處的事件(event)的每一個被記錄的信號之后而影響整個地震剖 面。這種不想要的能量相對于被記錄的反射通常是很高的幅度。源自于能量本身的特性 的這種效應(yīng)的另一術(shù)語是“虛反射(ghost)”能量。換句話說，在好的地震信號被記錄 的同時，水底信號以像虛反射一樣的方式緊跟在它們之后。這些虛反射事件可以構(gòu)造性 地結(jié)合達(dá)到幾乎無法檢測所需的儲集層信號的程度。特別地，在使用海底電纜的海洋勘 測中，由于電纜在海底的位置，截留在水層的能量至少部分地被記錄。這種不想要的能 量相對于被記錄的反射通常是很高的幅度。在這種類型的早期勘測期間，這些高的幅度 水平不能通過行業(yè)標(biāo)準(zhǔn)的單通道間隙去卷積來解決。這將海底電纜采集的使用限制于僅 在大部分淺的水深處，因為這有助于消除截留的水底能量。相同的問題使得海底電纜系統(tǒng)在過去20年間不斷發(fā)展。現(xiàn)有技術(shù)的地震數(shù)據(jù) 收集和處理系統(tǒng)已發(fā)展到嘗試處理這個問題。由于各種原因，這些系統(tǒng)沒有完全抑制虛 反射能量?，F(xiàn)代的電纜結(jié)構(gòu)具有最少兩個傳感器、水中聽音器(hydrophone)和速度檢波 器(velocity phone)，它們以一種極性記錄所需的反射信號而以相反的極性記錄不想要的 水底能量。在對這兩個傳感器求和時，理想的是水底能量相消而所需的信號增強。這種 雙傳感器設(shè)計允許電纜被部署至更深的水深處，這對于本行業(yè)是大的進步，因為其允許 在深水船先前無法接近的區(qū)域中進行采集。不幸的是，該求和過程的效率取決于許多因 素。這些因素包括電纜與海底耦合的良好程度、擱置電纜的沉積物的類型以及海底地形 的崎嶇程度。在求和之前，進行若干行業(yè)標(biāo)準(zhǔn)的處理步驟來對付這些問題。舉例來說，它們包括縮放、噪聲去除和波場分離。所有這些步驟都把水底能量作為目標(biāo)，但是沒有 完全抑制它。更糟的是，取決于所述求和的效率，殘余水平可能跨工程變化。行業(yè)標(biāo)準(zhǔn)的疊前(pre-stack)處理技術(shù)被用于把這種能量作為目標(biāo)，該技術(shù)包括 震源(source)、接收點(receiver)、中間點和偏移域中的地表一致性去卷積。但是，這樣 的技術(shù)不提供對留在所述數(shù)據(jù)中的殘余水底能量的測量。這些未知的、在空間上可變的 能量信號最終使得所述數(shù)據(jù)的判讀者(interpreter)對所述地震數(shù)據(jù)失去信任，因為觀察到 的可變幅度是由于殘余水底能量所引起的還是由于儲集層的重要的巖性特征所引起的將 是不清楚的。行業(yè)標(biāo)準(zhǔn)的疊后(post-stack)單通道去卷積也被采用。這樣的疊后單通道去卷積 通常被應(yīng)用在疊加數(shù)據(jù)上以減弱混響效應(yīng)(reverberation effect)，諸如殘余水底能量。在 這樣的技術(shù)中，根據(jù)疊加的地震道(seismic trace)(有時混有一個或兩個相鄰的道)上的數(shù) 據(jù)窗計算出的自相關(guān)通常被用在標(biāo)準(zhǔn)的預(yù)測去卷積(例如Weiner-Levinson)算法中以獲得 被應(yīng)用于所述數(shù)據(jù)的因果反向濾波器。但是，發(fā)明人認(rèn)識到由于殘余水底幅度水平跨勘 探區(qū)域變化，有時高于而有時低于被記錄的信號和背景噪聲水平，所以基于單一的自相 關(guān)的行業(yè)標(biāo)準(zhǔn)技術(shù)難以區(qū)分殘余水底能量與可能具有混響譜貌(reverberatory signature)的 其它所需的事件，諸如周期性的地質(zhì)反射物。結(jié)果是連同對所需信號的潛在抑制一起對 不想要的水底能量的次最佳抑制。相應(yīng)地，發(fā)明人認(rèn)識到現(xiàn)有的系統(tǒng)和技術(shù)/方法的缺陷包括雙傳感器海底電纜 數(shù)據(jù)的不恰當(dāng)?shù)那蠛?、低效的疊前或疊后預(yù)測去卷積算子設(shè)計以及缺少定位問題區(qū)域或 測量所述數(shù)據(jù)中的殘余幅度水平的質(zhì)量控制繪圖(quality control plot)。這些缺陷導(dǎo)致了 所述數(shù)據(jù)中的可變水平的殘余水底能量，這降低了所述數(shù)據(jù)的分辨能力和可判讀性。因此，發(fā)明人認(rèn)識到對可以從地震數(shù)據(jù)體(seismic data volume)中有效地找到和 抑制殘余水底能量，從而有效地提高地震數(shù)據(jù)的分辨能力一因此得到對從儲油層反射的 地震信號的改進的判讀(interpretation)的系統(tǒng)、程序產(chǎn)品和用計算機實現(xiàn)的方法的需要。 另外，所認(rèn)識到的是對有效地標(biāo)識殘余的多重能量的節(jié)省成本的系統(tǒng)、程序產(chǎn)品和用計 算機實現(xiàn)的方法的需要，其僅需要疊加數(shù)據(jù)，并且因此可以在臺式工作站上而不是在超 大型計算機系統(tǒng)上被執(zhí)行，依賴于疊前數(shù)據(jù)的多通道行業(yè)技術(shù)/方法通常需要超大型計 算機系統(tǒng)。發(fā)明人還認(rèn)識到對實現(xiàn)多通道疊后去卷積并且提供有效的質(zhì)量控制繪圖以允 許用戶在空間上分析地震數(shù)據(jù)中的殘余水底能量的水平的改進的系統(tǒng)、程序產(chǎn)品和用 計算機實現(xiàn)的方法，所述多通道疊后去卷積調(diào)節(jié)(leverage)縱測線(sub-line)、橫測線 (cross-line)和水底時間域中的冗余(redundancy)以得出用于地震數(shù)據(jù)的有效的預(yù)測去卷 積算子。所認(rèn)識到的是通過有效地抑制由殘余水底能量引起的反射事件上的旁瓣，這樣 的系統(tǒng)、程序產(chǎn)品和方法可以產(chǎn)生具有更高的分辨能力的地震體。將原始的地震數(shù)據(jù)轉(zhuǎn) 換為具有大體上更高的分辨能力的地震數(shù)據(jù)將為使用地震數(shù)據(jù)來估計儲集層性質(zhì)并且指 導(dǎo)水平鉆井作業(yè)的地震數(shù)據(jù)判讀者提供顯著的好處。

發(fā)明內(nèi)容
鑒于前述內(nèi)容，本發(fā)明的各實施例提供了解決地震數(shù)據(jù)處理中的問題的來源的系統(tǒng)、程序產(chǎn)品和用計算機實現(xiàn)的方法，并且提供了通過使所述過程適應(yīng)特定數(shù)據(jù)而解 決和克服殘余水底能量的有害的副作用的系統(tǒng)、程序產(chǎn)品和方法。本發(fā)明的各實施例包 括例如被設(shè)計用于有效地定位和抑制來自地震數(shù)據(jù)體內(nèi)的殘余水底能量的系統(tǒng)、程序產(chǎn) 品和方法。本發(fā)明的各實施例實現(xiàn)了被應(yīng)用于疊后數(shù)據(jù)集(post-stackdataset)的數(shù)據(jù)處理 進程。本發(fā)明的各實施例將多通道去卷積與新穎的分選鍵(sorting key)結(jié)合以用節(jié)省成 本的方式有效地標(biāo)識和抑制共深度點(CDP)疊加的地震數(shù)據(jù)中的殘余水底能量。本發(fā)明 的各實施例有利地提高了地震數(shù)據(jù)的分辨能力(resolving power)—因此得到對從儲油層反 射的地震信號的改進的判讀。具體地說，本發(fā)明的各實施例提供了抑制地震數(shù)據(jù)中的殘余水底能量的用計算 機實現(xiàn)/協(xié)助的方法。根據(jù)這樣的方法的實施例的示例，方法可以包括獲得地震的疊后 數(shù)據(jù)并且將多域去卷積應(yīng)用于所述地震的疊后數(shù)據(jù)以由此抑制地震數(shù)據(jù)體中的殘余水底 能量的步驟。這個步驟可以包括將共同水底時間與縱測線號(subline number)的結(jié)合用 作所述多域去卷積的自相關(guān)求平均基準(zhǔn)(autocorrelation averaging reference)以由此抑制 殘余水底能量檢測。在所述方法的實施例中，所述步驟可以進一步包括將地震數(shù)據(jù)疊加 到共深度點面元(bin)中以形成包括多個地震道的三維疊加數(shù)據(jù)體，在所述多個疊加的 地震道中的每一個上執(zhí)行多個自相關(guān)，在多個域中對所述相關(guān)數(shù)據(jù)求平均，響應(yīng)于所述 多個域中的平均相關(guān)數(shù)據(jù)而確定多個預(yù)測去卷積逆算子(predictive deconvolution inverse operator),將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)，以及響應(yīng)于將所述預(yù)測去卷 積逆算子應(yīng)用于所述地震數(shù)據(jù)而標(biāo)識存在于所述地震數(shù)據(jù)中的水底能量。所述方法還可 以包括計算機或判讀者使用質(zhì)量控制繪圖在空間上分析所述數(shù)據(jù)中的殘余水底能量的水 平。根據(jù)本發(fā)明的實施例的另一示例，所述方法可以包括(a)讀入在首部(header) 中有橫測線和縱測線號的疊加的地震數(shù)據(jù)；(b)讀入每一道的雙程水底時間(two-way water bottom time)并且將其放在所述首部中；(c)讀入每一道的地表下方層位時間 (sub-surface horizon time)并且將其存儲在所述首部中；(d)拉平(flatten)所述地震數(shù)據(jù) 使得基準(zhǔn)層位對于所述體中的所有道都處于恒定時間；(e)從拉平的地震體中提取時間 窗；(f)根據(jù)被提取的時間窗計算所述體中的每一道的自相關(guān)；(g)測量在所述雙程水底 時間處每一道上的自相關(guān)幅度水平并且將該值繪制成質(zhì)量控制圖；(h)使用來自步驟(f) 的輸出導(dǎo)出多個平均自相關(guān)；⑴將所述平均自相關(guān)保存到磁盤；(j)使用步驟(h)的平 均自相關(guān)中的至少一個將所述數(shù)據(jù)疊加到共同水底時間面元中并且將按道格式(in trace format)的結(jié)果顯示為質(zhì)量控制繪圖；(k)測量步驟(h)的自相關(guān)中的至少一個以得到在 對應(yīng)于所述雙程水底時間的時間樣點處的幅度值并且將結(jié)果繪制成質(zhì)量控制圖；⑴獲得 兩個因果反向濾波器以應(yīng)用于所述疊加體中的每一道；(m)將來自步驟⑴的算子應(yīng)用于 來自步驟(c)的數(shù)據(jù)并且將輸出保存到磁盤；(n)使用來自步驟(m)的輸出，重復(fù)步驟 (d)到(k)以評估步驟(m)的有效性；以及(0)改變步驟⑴中的算子設(shè)計參數(shù)并且根據(jù) 需要重新運行步驟(m)。本發(fā)明的實施例還可以包括被存儲在可由計算機讀取的有形的計算機可讀存儲 器或其它介質(zhì)中的程序產(chǎn)品，并且所述程序產(chǎn)品包括指令集使得在該指令集被所述計算 機執(zhí)行時使所述計算機執(zhí)行各種操作以抑制殘余水底能量。根據(jù)所述程序產(chǎn)品的實施例的示例，所述操作可以包括獲得地震數(shù)據(jù)以及將共同水底時間與縱測線號的結(jié)合用作多 域去卷積的自相關(guān)求平均基準(zhǔn)以由此抑制殘余水底能量。程序產(chǎn)品的另一實施例可以包 括在被執(zhí)行時使計算機操作以顯示可用于在空間上分析所述數(shù)據(jù)中的殘余水底能量的水 平的質(zhì)量控制繪圖的指令集。根據(jù)所述程序產(chǎn)品的另一實施例，所述操作可以包括將多域去卷積應(yīng)用于地震 的疊后數(shù)據(jù)以由此抑制地震數(shù)據(jù)體中的殘余水底能量的那些操作。這種操作可以對應(yīng)地 包括將共同水底時間與縱測線號的結(jié)合用作所述多域去卷積的自相關(guān)求平均基準(zhǔn)以由此 抑制殘余水底能量檢測。在所述程序產(chǎn)品的實施例中，所述操作可以進一步包括將地震 數(shù)據(jù)疊加到共深度點面元中以形成包括多個地震道的三維疊加數(shù)據(jù)體，在所述多個疊加 的地震道中的每一個上執(zhí)行多個自相關(guān)，在多個域中對所述相關(guān)數(shù)據(jù)求平均，響應(yīng)于所 述多個域中的平均相關(guān)數(shù)據(jù)而確定多個預(yù)測去卷積逆算子，將所述預(yù)測去卷積逆算子應(yīng) 用于所述地震數(shù)據(jù)，以及響應(yīng)于將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)而標(biāo)識存 在于所述地震數(shù)據(jù)中的水底能量。所述操作還可以包括顯示質(zhì)量控制繪圖以在空間上分 析所述數(shù)據(jù)中的殘余水底能量的水平。本發(fā)明的實施例還可以包括可由計算機讀取以分析地震數(shù)據(jù)的計算機可讀介 質(zhì)。根據(jù)本發(fā)明的實施例的示例，所述計算機可讀介質(zhì)可以包括指令集，該指令集在被 計算機執(zhí)行時使所述計算機執(zhí)行將多域去卷積應(yīng)用于地震的疊后數(shù)據(jù)的操作以由此抑制 地震數(shù)據(jù)體中的殘余水底能量的操作。這種操作可以包括將共同水底時間與縱測線號的 結(jié)合用作所述多域去卷積的自相關(guān)求平均基準(zhǔn)以由此抑制殘余水底能量檢測。所述操作 可以進一步包括將地震數(shù)據(jù)疊加到共深度點面元中以形成所述三維疊加數(shù)據(jù)體，在所述 疊加的數(shù)據(jù)體中的多個疊加的地震道中的每一個上執(zhí)行多個自相關(guān)，在多個域中對所述 相關(guān)數(shù)據(jù)求平均，響應(yīng)于所述多個域中的平均相關(guān)數(shù)據(jù)而確定多個預(yù)測去卷積逆算子， 將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)，以及響應(yīng)于將所述預(yù)測去卷積逆算子應(yīng) 用于所述地震數(shù)據(jù)而標(biāo)識存在于所述地震數(shù)據(jù)中的水底能量。所述操作還可以包括顯示 質(zhì)量控制繪圖以在空間上分析所述數(shù)據(jù)中的殘余水底能量的水平。本發(fā)明的各實施例不同于其它的行業(yè)處理技術(shù)。例如，不同于采用基于多重 (multi-fold)自相關(guān)的間隙去卷積的方法來解決水層能量問題的、炮點-接收點空間中的 預(yù)測去卷積，本發(fā)明的各實施例不使用疊前地表一致性方法，更不用說在炮點(shot)、 接收點、共深度點或偏移域中的疊前地表一致方法。另外，本發(fā)明的各實施例不需要使 用Split-Backus混響響應(yīng)來估計在每個震源(source)和接收點以下的混響時間。相反地， 本發(fā)明的各實施例有利地是在疊后的、大致零偏移的地表下方的域中，其使用與水底時 間結(jié)合的共同的縱測線、僅使用共同的縱測線和/或使用共同的橫測線。所述縱測線和 橫測線不同于中間點域。另外，在本發(fā)明的實施例中，有利的是水層時間估計來自于由 震源艇導(dǎo)航系統(tǒng)所記錄的水深信息。本發(fā)明的各實施例與地表一致性去卷積(疊前多通道去卷積)明顯不同之處在于 例如這樣的實施例通過將地表下方的基準(zhǔn)(sub-surface reference)與水底時間的結(jié)合用作 分解域來處理海洋數(shù)據(jù)(marine data)以及截留在水層中的能量的特性。另外，本發(fā)明的 各實施例不同于Log/Fourier域中的地表一致性去卷積，因為它們不在震源、接收點、結(jié) 構(gòu)和偏移域中假定疊前地表一致性的前提。相反地，本發(fā)明的各實施例在由縱測線與水底時間的結(jié)合、所述縱測線和所述橫測線所定義的疊后的地表下方的域中操作。本發(fā)明的各實施例也不同于去卷積地震濾波。舉例來說，雖然這樣的實施例從 每一地震道中提取自相關(guān)，但是所述去卷積算子不是從相同的自相關(guān)中導(dǎo)出的。相反 地，它是在通過求平均技術(shù)處理所述自相關(guān)，接著為每個地震道基于表示對于該道位置 唯一的縱測線(或橫測線)和水底時間的自相關(guān)的結(jié)合構(gòu)造濾波器之后被導(dǎo)出的。根據(jù)行業(yè)的處理技術(shù)，在疊前運行多通道去卷積，而在疊后運行單通道去卷 積。典型的疊前多通道去卷積(也被稱為地表一致性去卷積)使用不同的分選鍵，這些 分選鍵通常是震源、接收點、中間點和偏移的結(jié)合。有利的是本發(fā)明的各實施例不需要 通過這些鍵(key)中的任一個來引用(referto)所述數(shù)據(jù)。另外，有利的是通過使用疊后 多通道去卷積，本發(fā)明的各實施例可以利用所述數(shù)據(jù)中的某些冗余以由此暴露所述殘余 水底能量并且增強對其的抑制。


為了更詳細(xì)地理解本發(fā)明的特征和優(yōu)點以及將變得顯而易見的其它特征和優(yōu) 點，對上文所簡要概括的發(fā)明的更具體的說明可以通過參考其實施例來得到，所述實施 例在附圖中被示意，附圖構(gòu)成本說明書的一部分。但還應(yīng)注意的是附圖僅示意了本發(fā)明 的各實施例并且因此不應(yīng)被認(rèn)為是限制本發(fā)明的范圍，因為它還可以包括其它有效的實 施例。圖1是根據(jù)本發(fā)明的實施例的用于抑制殘余水底能量的系統(tǒng)的一部分的示意框 圖；圖2是根據(jù)本發(fā)明的實施例的抑制殘余水底能量的過程的流程框圖；圖3是根據(jù)本發(fā)明的實施例的從在抑制殘余水底能量的方法之前被示意的輸入 共深度點(CDP) 3D疊加數(shù)據(jù)中被提取的、雙程時間與共深度點垂直剖面相對的圖表；圖4是根據(jù)本發(fā)明的實施例的以毫秒為單位的雙程水底時間與以秒為單位的自 相關(guān)時間相對的圖表，其示意所有縱測線的共同水底時間的自相關(guān)疊加；圖5A是根據(jù)本發(fā)明的實施例的在雙程水底時間處的單道自相關(guān)幅度的圖形映像 (graphical map)；圖5B是根據(jù)本發(fā)明的實施例的以毫秒為單位的雙程水底時間的圖形映像；圖6A是根據(jù)本發(fā)明的實施例的縱測線位置與以毫秒為單位的雙程水底時間相對 的質(zhì)量控制圖表，其示意研究(study)區(qū)域中的每個縱測線的在共同水底時間處的自相關(guān) 疊加幅度；圖6B是根據(jù)本發(fā)明的實施例的圖6A中的標(biāo)繪值的幅度直方圖；圖7A是根據(jù)本發(fā)明的實施例的以毫秒為單位的雙程水底線與以秒為單位的自相 關(guān)滯后時間(lagtime)相對的圖表，其示意輸入數(shù)據(jù)的與整個研究的水底時間相對的自相 關(guān)疊加；圖7B是根據(jù)本發(fā)明的實施例的以毫秒為單位的雙程水底線與以秒為單位的自相 關(guān)滯后時間相對的圖表，其示意在應(yīng)用殘余水底能量抑制方法的實施例之后輸入數(shù)據(jù)的 與整個研究的水底時間相對的自相關(guān)疊加；圖8A是根據(jù)本發(fā)明的實施例的在應(yīng)用抑制方法的實施例之前在水底時間處的自相關(guān)幅度的圖形映像；圖8B是根據(jù)本發(fā)明的實施例的在應(yīng)用抑制方法的實施例之后在水底時間的自相 關(guān)幅度的圖形映像；圖9A是以毫秒為單位的雙程水底時間與以秒為單位的自相關(guān)滯后時間相對的圖 表，其示意輸出數(shù)據(jù)的跨整個研究區(qū)域與水底時間相對的自相關(guān)疊加；圖9B是根據(jù)本發(fā)明的實施例的以毫秒為單位的雙程水底時間與以秒為單位的自 相關(guān)滯后時間相對的圖表，其示意在應(yīng)用抑制方法的實施例之后的輸出數(shù)據(jù)的跨整個研 究區(qū)域與水底時間相對的自相關(guān)疊加；圖9C是以毫秒為單位的雙程水底時間與以秒為單位的自相關(guān)滯后時間相對的圖 表，其示意在使用標(biāo)準(zhǔn)單道間隙去卷積的情況下的輸出數(shù)據(jù)的跨整個研究區(qū)域與水底時 間相對的自相關(guān)疊加；圖10是根據(jù)本發(fā)明的實施例在應(yīng)用殘余水底能量抑制方法的實施例之后的輸出 雙程時間垂直剖面的圖表；圖11A是根據(jù)本發(fā)明的實施例的在應(yīng)用抑制殘余水底能量的方法的實施例之前 的輸入剖面；以及圖11B是根據(jù)本發(fā)明的實施例的在應(yīng)用抑制殘余水底能量的方法的實施例之后 的輸出剖面。
具體實施例現(xiàn)在將在下文中參考附圖更全面地描述本發(fā)明，附圖示意了本發(fā)明的實施例。 但是，本發(fā)明可以用許多不同形式來具體化而不應(yīng)當(dāng)被看作受限于在此所闡述的被示意 的實施例；相反地，這些實施例被提供使得本公開內(nèi)容將是詳細(xì)而完整的，并且將向本 領(lǐng)域的技術(shù)人員全面地傳達(dá)本發(fā)明的范圍。相似的標(biāo)號通篇都指的是相似的元件。注 意，在此所使用的術(shù)語“虛反射分析和終止”和縮寫“GHOAT”應(yīng)被理解為是對根據(jù)本 發(fā)明的一個或多個實施例的過程不完整的引用(short-handedreference)。在此所使用的術(shù) 語“首部(header)”包括與疊加的地震道相關(guān)聯(lián)的數(shù)據(jù)的存儲位置。該數(shù)據(jù)可以指反應(yīng) 地震道在工程(project)內(nèi)的位置的編號(number)或者可以指與該具體的道位置相關(guān)聯(lián)的 值，諸如水底時間。其它的道數(shù)據(jù)也可以被存儲在該首部中。圖1-11B示意了根據(jù)本發(fā)明的示范性實施例的被設(shè)計用于從地震數(shù)據(jù)體中有效 地找到并且抑制殘余水底能量的、示范性的地震數(shù)據(jù)分析系統(tǒng)12、地震數(shù)據(jù)分析程序產(chǎn) 品13以及方法。如在圖1中可能是最佳地被示出的那樣，系統(tǒng)12可以包括地震數(shù)據(jù)分 析計算機14，該地震數(shù)據(jù)分析計算機具有處理器15和與處理器15通信的存儲器16以在 其中存儲運行指令(operating instruction)并且存儲地震數(shù)據(jù)分析程序產(chǎn)品13。系統(tǒng)12還 包括數(shù)據(jù)庫17，該數(shù)據(jù)庫能被地震數(shù)據(jù)分析計算機14的處理器15訪問并且具有包括地震 數(shù)據(jù)集的多個數(shù)據(jù)記錄，所述地震數(shù)據(jù)集可以包括疊加和未疊加的地震道兩者，每個地 震道與包含相應(yīng)的地震道的屬性的首部相關(guān)聯(lián)和/或包含該首部。每個地震道的這樣的 首部數(shù)據(jù)可以包括縱測線號和橫測線號，并且可以包括或使其包括雙程水底時間，以及 其它變量或標(biāo)識符。首部數(shù)據(jù)還可以包括或使其包括地表下方基準(zhǔn)層位時間。系統(tǒng)12 還可以包括計算機系統(tǒng)領(lǐng)域的技術(shù)人員所知道的各種輸入設(shè)備、顯示器、網(wǎng)絡(luò)連接設(shè)備等等(未示出)。系統(tǒng)12還可以包括本領(lǐng)域的技術(shù)人員所知道的各種裝置(未示出)以 收集或者產(chǎn)生地震數(shù)據(jù)。例如，這樣的數(shù)據(jù)可以通過記錄來自傳感器對(舉例來說諸如 水中聽音器傳感器)和來自速度檢波器的反射信號而產(chǎn)生，它們接收例如從海底電纜所 產(chǎn)生的反射信號。注意，如計算機系統(tǒng)領(lǐng)域的技術(shù)人員還知道并且理解的那樣，存儲器16可以包 括本領(lǐng)域的技術(shù)人員所知道的易失和非易失存儲器，包括例如RAM、ROM和磁盤或光 盤，這僅是舉幾個例子。還要注意，地震數(shù)據(jù)分析程序產(chǎn)品13可以是以提供特定的一組 或多組有序操作的微代碼、程序、例程和符號語言的形式的，所述有序操作控制硬件的 運行并且指導(dǎo)其操作，如本領(lǐng)域的技術(shù)人員所知道并且理解的那樣。另外注意，根據(jù)本 發(fā)明的實施例的示例，地震數(shù)據(jù)分析程序產(chǎn)品13不需要全部存在于易失存儲器中，而可 以根據(jù)本領(lǐng)域的技術(shù)人員所知道并且理解的各種方法按需要選擇性地被加載。另外，還 應(yīng)理解的是被示意的計算機配置是通過示例的方式給出的，而根據(jù)本領(lǐng)域的技術(shù)人員所 知道的各種其它方法所配置的其它類型的計算機可以被使用。但是，計算機14優(yōu)選地采 用計算機系統(tǒng)領(lǐng)域的技術(shù)人員所知道的臺式、膝上型或其它半便攜式或便攜式計算機的 形式，而不是采用大型計算機系統(tǒng)的形式，依賴于疊前數(shù)據(jù)的多通道的行業(yè)技術(shù)/方法 通常需要大型計算機系統(tǒng)。本發(fā)明的以下示范性實施例被描述，它們被應(yīng)用于疊后數(shù)據(jù)集。圖3是從輸入 的3D疊加體中被提取的雙程時間垂直剖面的示例。反向箭頭61突出了具有差的分辨率 和高的旁瓣的事件，這可能是由于數(shù)據(jù)中的殘余水底能量的“虛反射”效應(yīng)所引起的。 本發(fā)明的各實施例將多通道去卷積與唯一的分選鍵唯一地結(jié)合以有效地標(biāo)識并且抑制共 深度點(CDP)疊加的地震數(shù)據(jù)中的殘余水底能量。例如三個鍵中的第一個是地震道首部 的結(jié)合(例如共同縱測線乘以100加上以毫秒為單位的共同水底時間)，其跨該工程以可 調(diào)整的方式有效地隔離了水底能量。另外兩個鍵是縱測線和橫測線。本發(fā)明的各實施例 有利地不需要將共同中間點用作分選鍵。如將在下面更詳細(xì)地被描述的那樣，在導(dǎo)出每 個鍵的平均自相關(guān)之后，通常只有兩個隨后被用于去卷積。這兩個鍵包括與共同水底時 間項結(jié)合的、結(jié)合的共同縱測線以及單獨的共同縱測線的項(term)。常規(guī)地，僅在CDP疊加前運行多通道去卷積，而在疊后運行單通道去卷積。但 是，通過使用根據(jù)本發(fā)明的實施例的疊后多通道去卷積，用戶可以利用這樣的疊后數(shù)據(jù) 中的某些冗余，它們可以被用來暴露要被抑制的殘余水底能量。因此，現(xiàn)在有可能的是 有效地并且快速地隔離這種不想要的能量的周期性和幅度。利用多域去卷積來抑制這種 不定的能量而同時在目標(biāo)層位保持結(jié)構(gòu)和子波的完整性。如在圖2中可能最佳地顯示的那樣，根據(jù)本發(fā)明的實施例的示例，虛反射分析 和終止(GHOAT)方法被實現(xiàn)為一系列步驟和/或操作，其可以采用能夠被本領(lǐng)域的技 術(shù)人員所知道的計算機系統(tǒng)執(zhí)行的、市場上可買到的軟件程序和模塊。這樣的軟件的 示例包括可從在得克薩斯州的休斯敦有辦事處的Paradigm Geotechnology BV買到的名為
“Disco”的軟件。特別地，圖2示意了根據(jù)本發(fā)明的用于抑制殘余水底能量的方法、程 序產(chǎn)品和系統(tǒng)的實施例的示例的流程圖，并且可以采用以下步驟在框21處，計算機14(或其它分析器)讀入在首部中有縱測線和橫測線號的疊 加體(例如三維的疊加體)，并且作為連續(xù)的步驟或者作為隨后執(zhí)行的獨立的步驟，計算機14讀入每一道的地表下方基準(zhǔn)層位時間并且將其存儲在首部中。在框23處，計算機 14讀入/加載每一道的雙程水底時間并且將其放在首部中。在框27處，地震數(shù)據(jù)被拉平 使得基準(zhǔn)層位對于所述體中的所有道都處于恒定時間，并且從被拉平的地震體中提取時 間窗。在框31處，一旦地震數(shù)據(jù)已被疊加到共深度點面元中，就根據(jù)被提取的時間窗 計算所述體中的每一道的自相關(guān)。注意，所述步驟/操作還可以包括根據(jù)縱測線和橫測 線首部創(chuàng)建虛炮點(dummy shot point)和接收點首部，并且根據(jù)共同縱測線號和共同水底 時間創(chuàng)建虛偏移(dummy offset)首部。在框33處，測量在雙程水底時間處每一道上的自相關(guān)幅度水平，并且按顏色將 該值繪制成質(zhì)量控制圖，其中縱測線在垂直軸上而橫測線在水平軸上。圖5A提供了這種 繪圖的示例。注意，水底幅度中大的變化對應(yīng)于圖5B所示的水底時間的改變，圖5提供 了雙傳感器沒有有效地抑制水底能量的指示。在框35處，使用框31處的自相關(guān)分析的輸出和/或虛炮點和接收點首部以及虛 偏移首部導(dǎo)出平均自相關(guān)。來自每個疊加的道的自相關(guān)被收集(gather)并且在多個域中被 求平均，這允許用戶有效地標(biāo)識可能存在于所述數(shù)據(jù)中的水底能量。例如，可以基于地 震道首部輸出三個平均數(shù)據(jù)集/域，所述自相關(guān)在這三個平均數(shù)據(jù)集/域中被分解。舉 例來說，這些平均數(shù)據(jù)集/域被定義為(1)共同縱測線號乘以100加上以毫秒為單位的共 同水底時間，(2)共同縱測線號，以及(3)共同橫測線號。換句話說，所述自相關(guān)集合 (autocorrelationgather)根據(jù)計算平均自相關(guān)的子步驟被形成，其具有相同的(1)共同縱測 線號乘以100加上以毫秒為單位的共同水底時間，(2)共同縱測線號，以及(3)共同橫測 線號。假定這三“項”是卷積的，則輸入可以在對數(shù)頻譜域中由這三個項的線性組合表 示，其中每個項最后通過最小二乘誤差最小化被導(dǎo)出。根據(jù)本發(fā)明的實施例，已發(fā)現(xiàn)該 三項自相關(guān)的解產(chǎn)生優(yōu)于現(xiàn)有技術(shù)的對殘余水底能量更有效的標(biāo)識。平均自相關(guān)被保存 到磁盤。如在框49、50和51處所指示的那樣，在導(dǎo)出每個域/鍵的平均自相關(guān)之后，舉 例來說所述“項”或“鍵”中的兩個或更多個可以隨后被用于去卷積例如(1)縱測線 與水底時間結(jié)合的項；以及(2)只有縱測線的項。這類去卷積被稱為地表下方一致性去 卷積。在框45處，使用來自在框35處被標(biāo)識的步驟的部分/子步驟(1)的平均自相 關(guān)，所述數(shù)據(jù)被疊加到共同水底時間面元中，并且結(jié)果按道格式被顯示為質(zhì)量控制繪 圖。圖4提供了通過將單重(single fold)自相關(guān)加到共同水底時間面元中(例如對于整個 勘測)所形成的自相關(guān)分析的示例。不想要的水底能量的首次和第二次到達(dá)(arrival)清 晰可見，如虛線71、75和箭頭73、77所指示的那樣。通常不預(yù)期剛好高于感興趣區(qū)域 的這些水平的水底能量如此易于檢測。此處，自相關(guān)的零滯后時間在0.12秒處被顯示。 使用縱測線與水底時間結(jié)合的鍵的質(zhì)量控制繪圖可以展示這些鍵或域的冗余優(yōu)點。在框47處，來自在框35處被標(biāo)識的步驟的部分/子步驟(1)的自相關(guān)被測量以 得到時間樣點處的幅度值，所述時間樣點例如對應(yīng)于雙程水底時間。按顏色將結(jié)果繪制 成質(zhì)量控制圖，其中縱測線在一個軸上(例如Y軸)而共同水底時間在另一軸上(例如 X軸)。圖6A提供了這種繪圖的示例，其示意了對于共同水底時間跨縱測線的可變幅度(箭頭)，從而指示需要在空間上可變的去卷積算子。圖6A-B對于顯示相同的殘留能量 跨勘測的空間分布是有用的。在框49和50處，對于所述疊加體中的每一道，來自在框35處被標(biāo)識的步驟的 部分/子步驟(1)和(2)的兩個平均自相關(guān)(即共同縱測線和水底時間分解)可以被用作 預(yù)測去卷積(例如Weiner-Levinson)算法的輸入以獲得兩個因果反向濾波器而應(yīng)用于每一 道。同樣地，來自在框45處被標(biāo)識的步驟的繪圖可以被用于確定/規(guī)定長度和間隙算子 參數(shù)。在框51處，在框49處被標(biāo)識的步驟處計算出的算子可以被應(yīng)用于每一道的地表 下方基準(zhǔn)層位時間并且輸出的三維疊加體可以被保存到磁盤。在框53處，使用來自框51 的三維疊加體輸出，可以重復(fù)相對于框27到49被標(biāo)識的步驟以評估在框51處被標(biāo)識的 去卷積算子的應(yīng)用的有效性。如果有必要，可以改變(重新規(guī)定)在框49處被標(biāo)識的算 子設(shè)計參數(shù)并且重新運行在框51處被標(biāo)識的步驟。如先前所提到的，根據(jù)現(xiàn)有技術(shù)，在共深度點(CDP)疊加之前或在共深度點 (CDP)疊加前專有地運行多通道去卷積，而在疊后運行單通道去卷積。根據(jù)本發(fā)明的各 實施例，通過代替地使用疊后多通道去卷積，已發(fā)現(xiàn)這樣的方法有利地標(biāo)識所述數(shù)據(jù)中 的某些冗余，這可以被用來更完全地暴露想要抑制的殘余水底能量而不需要大量資源。 表示共同縱測線(例如被乘以100)加上共同水底時間(例如以毫秒為單位)的疊加自相 關(guān)可以用多種方式來繪制以跨勘測用空間方式觀察幅度水平。另外，上述基于自相關(guān)分 解的質(zhì)量控制繪圖允許水底能量以有用的方式的可視化。這些繪圖可以幫助預(yù)測去卷積 逆算子的設(shè)計。在應(yīng)用所述算子之后，可以重新產(chǎn)生所述繪圖以觀察所述過程如何有效地起作 用。例如，圖7A和7B示意了所述數(shù)據(jù)的與整個研究的水底時間相對的自相關(guān)疊加的前 后比較。圖7A是在應(yīng)用參考圖4所論述的用于抑制水底能量的過程的實施例之前。而 圖7B是在應(yīng)用根據(jù)本發(fā)明的實施例的用于抑制水底能量的過程的實施例之后。如箭頭 91、93所指示的那樣，首次和第二次水底到達(dá)被有效地被抑制。圖8A和8B是在水底時間處的自相關(guān)幅度的圖表。圖8A是在應(yīng)用抑制方法的實 施例之前，而圖8B是在應(yīng)用抑制方法的實施例之后在水底時間處的自相關(guān)幅度的映像。 注意，通過抑制，本發(fā)明的被示意的實施例均衡了跨研究區(qū)域的不想要的能量(參見相 應(yīng)的直方圖)，除了最淺的時間(箭頭101)之外，在該處對背景地震子波的干擾導(dǎo)致了調(diào) 諧效應(yīng)(tuning effect)。圖9A、9B和9C是地震數(shù)據(jù)的跨整個研究區(qū)域與水底時間相對的自相關(guān)疊加的 比較。水平軸是向右增加的以毫秒為單位的雙程水底時間。垂直軸是從下到上增加的以 秒為單位的自相關(guān)滯后時間。該面板是顯示標(biāo)準(zhǔn)化的以對于每個面板的最大零滯后值都 具有相同的道偏差(trace deflection)。圖9A是使用輸入數(shù)據(jù)的比較。圖9B是使用應(yīng)用 抑制過程的實施例之后的數(shù)據(jù)的比較。圖9C是使用通過標(biāo)準(zhǔn)的單道間隙去卷積的數(shù)據(jù)的 比較。注意，殘余水底抑制過程的各實施例(參見例如圖9B)在使水底能量衰減方面顯 著地更加有效，如由圖9C所示的箭頭111所標(biāo)識的那樣。圖10示意了在應(yīng)用殘余水底能量抑制的實施例之后的輸出雙程時間垂直剖面， 這在圖3中被介紹。注意，被示意的反向箭頭61'突出了具有改進的分辨率和低旁瓣能
量的事件。
圖11A和11B進一步示意了由于應(yīng)用殘余水底能量抑制方法的實施例而帶來的 改進的分辨率和低旁瓣能量。換句話說，圖11A提供了在應(yīng)用抑制殘余水底能量的方法 的實施例之前的輸入剖面的示例。比較起來，圖11B提供了在應(yīng)用抑制殘余水底能量的 方法的實施例之后的輸出剖面。根據(jù)本發(fā)明的實施例的各種數(shù)據(jù)繪圖有利地向地震數(shù)據(jù)的判讀者提供了更高的 置信度，因為它們跨勘測在空間上清晰地表示了殘余水底能量的最終水平，允許他或她 確定觀察到的幅度變化是由于殘余水底能量所引起的還是由于儲集層的重要特征所引起 的?？傊?，通過將共同水底時間與縱測線號的結(jié)合用作多域去卷積的自相關(guān)求平均 基準(zhǔn)，連同新穎的質(zhì)量控制繪圖一起，可以展示與被用于疊后數(shù)據(jù)的標(biāo)準(zhǔn)的單通道去卷 積相比對殘余水底能量更有效的抑制。除了有效地標(biāo)識殘余的多重能量以外，所述方法 的各實施例已被證明是格外節(jié)省成本的，因為只需要疊加的數(shù)據(jù)。臺式工作站而不是大 得多的系統(tǒng)可以容易地實現(xiàn)所述方法的各實施例。比較起來，大多數(shù)多通道的行業(yè)方法 依賴于疊前數(shù)據(jù)，該疊前數(shù)據(jù)只能用超大型計算機系統(tǒng)來管理。如本領(lǐng)域的普通技術(shù)人員將理解的那樣，由本發(fā)明的系統(tǒng)、程序產(chǎn)品和方法的 實施例產(chǎn)生的在圖7A-11B中被呈現(xiàn)的比較結(jié)果顯示了明顯的改進并且示意了殘余水底能 量已有效率地并且有效果地被抑制。本發(fā)明的方法的實施例可以通過對現(xiàn)有軟件程序進 行微小的修改而有利地被實現(xiàn)，并且通過組裝所需要的軟件模塊，其可以利用來自市場 上可買到的地震數(shù)據(jù)處理軟件庫的實用程序?？蓮纳虡I(yè)軟件供應(yīng)商買到供本發(fā)明的各實 施例使用的合適的地震處理模塊，舉例來說諸如可從在得克薩斯州的休斯敦有辦事處的 Paradigm Geotechnology BV 買至Ll的 Disco。抑制地震數(shù)據(jù)中的殘余水底能量檢測的方法的實施例例如可以包括(a)讀入在 首部中有橫測線和縱測線號的疊加的地震數(shù)據(jù)，(b)讀入每一道的雙程水底時間并且將其 放在所述首部中，(c)讀取每一道的地表下方層位時間并且將其存儲在所述首部中，(d) 拉平地震數(shù)據(jù)使得基準(zhǔn)層位對于所述體中的所有道都處于恒定時間，(e)從被拉平的地震 體中提取時間窗，(f)根據(jù)被提取的時間窗計算所述體中的每一道的自相關(guān)，(g)測量在 雙程水底時間處每一道上的自相關(guān)幅度水平并且將該值繪制成質(zhì)量控制圖，以及(h)使 用來自步驟(f)的輸出導(dǎo)出平均自相關(guān)。步驟(h)可以包括展開例如三個數(shù)據(jù)集。這三 個平均數(shù)據(jù)集可以基于道首部被輸出，舉例來說被定義為(1)共同縱測線乘以100加上以 毫秒為單位的共同水底時間，(2)共同縱測線號，以及(3)共同橫測線號。所述方法還可以包括⑴將平均自相關(guān)保存到磁盤，(j)使用來自步驟(h)的 部分1的平均自相關(guān)將所述數(shù)據(jù)疊加到共同水底時間面元中并且將按道格式的結(jié)果顯示為 質(zhì)量控制繪圖，(k)測量來自步驟(h)的部分⑴的自相關(guān)以得到在對應(yīng)于雙程水底時間 的時間樣點處的幅度值并且將結(jié)果繪制成質(zhì)量控制圖，⑴獲得兩個因果反向濾波器以應(yīng) 用于所述疊加體中的每一道，(m)將來自步驟⑴的算子應(yīng)用于來自步驟(c)的數(shù)據(jù)并且 將輸出保存到磁盤，(n)使用來自步驟(m)的輸出，重復(fù)步驟(d)到(k)以評估步驟(m) 的有效性，以及(o)改變步驟⑴中的算子設(shè)計參數(shù)并且根據(jù)需要重新運行步驟(m)。對于具有疊加的三維地震數(shù)據(jù)的地震數(shù)據(jù)，步驟(h)可以包括基于道首部被輸 出的三個平均數(shù)據(jù)集，其被定義為(1)共同縱測線乘以100加上以毫秒為單位的共同水底時間，(2)共同縱測線號，以及(3)共同橫測線號，而步驟⑴可以包括將來自步驟(h) 中的部分(1)和(2)的兩個平均自相關(guān)用作標(biāo)準(zhǔn)的預(yù)測去卷積算法的輸入并且使用來自步 驟(j)的繪圖來確定算子的長度和間隙參數(shù)。抑制地震數(shù)據(jù)中的殘余水底能量的方法的實施例的示例可以包括將共同水底時 間與縱測線號的結(jié)合用作多域去卷積的自相關(guān)求平均基準(zhǔn)以由此抑制殘余水底能量檢 測。方法的實施例的另一示例包括使用質(zhì)量控制繪圖在空間上分析所述數(shù)據(jù)中的殘余水 底能量的水平?？捎捎嬎銠C讀取的被存儲在有形的計算機可讀存儲器中的程序產(chǎn)品的實施例被 提供。根據(jù)本發(fā)明的實施例的示例，所述程序產(chǎn)品可以包括指令集使得在該指令集被計 算機讀取時使所述計算機執(zhí)行將共同水底時間與縱測線號的結(jié)合用作多域去卷積的自相 關(guān)求平均基準(zhǔn)的操作以由此抑制殘余水底能量檢測。程序產(chǎn)品的實施例的另一示例包括 使所述計算機執(zhí)行產(chǎn)生并且顯示質(zhì)量控制繪圖以在空間上分析所述數(shù)據(jù)中的殘余水底能 量的水平的操作的指令集。重要的是應(yīng)注意雖然本發(fā)明的實施例已在全功能系統(tǒng)的背景下被描述，但是本 領(lǐng)域的技術(shù)人員將理解本發(fā)明和/或其方面的至少部分的機制能夠以指令的計算機可讀 介質(zhì)的形式分布，所述指令以各種形式用于在處理器、多個處理器或諸如此類上的執(zhí) 行，并且無論被用于實際執(zhí)行所述分布的信號承載介質(zhì)的特定類型如何，本發(fā)明的實施 例都同樣適用。計算機可讀介質(zhì)的示例包括但不限于非易失的、硬編碼型介質(zhì)，諸 如只讀存儲器(ROM)、CD-ROM和DVD-ROM，或者可擦除的、電可編程的只讀存 儲器(EEPROM)、可記錄型介質(zhì)，諸如軟盤、硬盤驅(qū)動器、CD-R/RW、DVD-RAM、 DVD-R/RW、DVD+R/RW,閃存驅(qū)動器和其他更新的類型的存儲器，以及諸如數(shù)字和 模擬通信鏈路的傳輸型介質(zhì)。例如，這樣的介質(zhì)可以包括與系統(tǒng)軟件、地震數(shù)據(jù)分析程 序產(chǎn)品13以及在上文中所描述的 方法和過程步驟有關(guān)的運行指令和操作指令兩者。本申請涉及Burnstad等人的、申請?zhí)枮?1/032,637、申請日為2008年2月29 日的、題為"Method, Program, Product, and System for Suppression of Residual Water Bottom Energy in Surface Seismic Data”的美國臨時申請，以及Burnstadt等人的、申請?zhí)?為 61/033,012、申請日為 2008 年 3 月 2 日的、題為 ‘‘Method，Program, Product, and System for Suppression of Residual Water Bottom Energy in Surface Seismic Data” 的美國臨 時申請，通過引用將其每個都整體并入。以下參考文獻(xiàn)每個都通過引用在此整體被并入Morley, L.禾口 Claerbout, J. (1983 年 5 月)Predictive Evolution in Shot-Receiver Space ； Geophysics, 48, pp.515-531.Levin, S.A. (1989 年 9 月)Surface—Consistent Deconvolution ; Geophysics, 54, pp.1123-1133.Cambois, G.禾口 Stoffaa, P. (1993 年 8 月)Surface—Consistent Phase Decomposition in the Log/Fourier Domain ； Geophysics, 58, pp.1099-1111.Duren 的、專利號為 6,678,207、題為 “ Trapped Water Bottom Multiple and Peg-Laden Multiple Suppression for Ocean Bottom Seismic Data” 的美國專禾丨J。Burg 的、專利號為 3,512,127、題為 “Deconvolution Seismic Filtering” 的美國專利。Dobrin 等人的、專利號為 3,370,268、題為 “Method of Processing Geological and Geophysical Data” 的美國專利。Starr 的、專利號為 6,263,285、題為 “Amplitude Spectra Estimation，，的美國專利。Starr 的、專利號為 6,151,275、題為 “Method of Dual Wavefield Reinforcement”
的美國專利。Starr 的、專利號為 6,026,059、題為 “Method of Creating Common-Offset/ Common-Azimuth Gathers in 3-D Seismic Surveys and Method of Conducting Reflection Attribute Variation Analysis ” 的美國專利。McCormick 等人的、專利號為 6,154J05、題為 “ System for Attenuating High Order Free Surface Multiples from a Seismic Shot Record Using a Genetic Procedure” 的美國專利。在附圖和說明書中已公開了本發(fā)明的典型的優(yōu)選實施例，而雖然采用了特定的 術(shù)語，但這些術(shù)語僅在描述的意義上被使用而不是為了限制的目的。已通過具體參考這 些被示意的實施例非常詳細(xì)地描述了本發(fā)明。但是，將顯而易見的是可以在前面的說明 書所描述的本發(fā)明的精神和范圍內(nèi)進行各種修改和變化。
權(quán)利要求
1.一種抑制地震數(shù)據(jù)中的殘余水底能量檢測的用計算機實現(xiàn)的方法，所述方法包括 以下步驟獲得定義疊后數(shù)據(jù)體的共深度點疊加的三維地震數(shù)據(jù)，所述疊后數(shù)據(jù)體包括 多個地震道，每個所述地震道包括首部，該首部包括橫測線和縱測線號；獲得所述多個 地震道中的每一個的雙程水底時間數(shù)據(jù)；以及獲得所述多個地震道中的每一個的地表下 方層位時間數(shù)據(jù)，所述方法的特征在于以下步驟將所述多個地震道中的每一個的所獲得的相應(yīng)的雙程水底時間數(shù)據(jù)存儲在相應(yīng)的首 部中；將所述多個地震道中的每一個的所獲得的相應(yīng)的地表下方層位時間數(shù)據(jù)存儲在相應(yīng) 的首部中；拉平所述疊加的地震數(shù)據(jù)以定義拉平的地震數(shù)據(jù)體使得基準(zhǔn)層位對于所述疊后數(shù)據(jù) 體中的多個道中的每一個都處于恒定時間； 從拉平的地震數(shù)據(jù)體中提取時間窗；響應(yīng)于所述體中相應(yīng)的地震道的相應(yīng)的被提取的時間窗而計算所述體中的多個地震 道中的每個單獨的地震道的自相關(guān)；響應(yīng)于所述多個自相關(guān)而在多個不同的域上執(zhí)行多域自相關(guān)分析；以及 使用來自所述多個不同的域中的至少兩個的定義分選鍵的項將多域去卷積應(yīng)用于所 述地震的疊后數(shù)據(jù)以由此標(biāo)識和抑制所述殘余水底能量。
2.如權(quán)利要求1所定義的方法，其中所述執(zhí)行多域自相關(guān)分析的步驟包括 響應(yīng)于所述多個自相關(guān)而導(dǎo)出多個平均自相關(guān)，所述導(dǎo)出包括在多個不同的域的每一個中為所述多個自相關(guān)中的每一個確定平均值。
3.如權(quán)利要求2所定義的方法，其中所述在多個不同的域的每一個中為所述多個自相 關(guān)中的每一個確定平均值的步驟包括響應(yīng)于所述多個地震道中的每一個的首部數(shù)據(jù)而輸出三個平均數(shù)據(jù)集，所述數(shù)據(jù)集 包括(1)共同縱測線號和共同水底時間，(2)所述共同縱測線號，以及(3)共同橫測線 號。
4.如權(quán)利要求2或3中的任意一項所定義的方法，其特征還在于以下步驟 響應(yīng)于所述多個導(dǎo)出的平均自相關(guān)中的至少一個而將所述疊后地震數(shù)據(jù)疊加到共同水底時間面元中；以及響應(yīng)于所述將所述疊后地震數(shù)據(jù)疊加到共同水底時間面元中的步驟而顯示按道格式 的疊后地震數(shù)據(jù)的質(zhì)量控制圖。
5.如權(quán)利要求4所定義的方法，其中所述將多域去卷積應(yīng)用于所述地震的疊后數(shù)據(jù)的 步驟包括響應(yīng)于按道格式的疊后地震數(shù)據(jù)的質(zhì)量控制圖而接收去卷積算子的參數(shù)； 響應(yīng)于所述多個導(dǎo)出的平均自相關(guān)中的至少兩個和所述去卷積算子的參數(shù)而獲得確 定一個或多個去卷積算子的兩個因果反向濾波器來以應(yīng)用于所述疊后數(shù)據(jù)體中的多個地 震道中的每一個；以及將所述去卷積算子應(yīng)用于所述多個地震道中的每一個的地表下方層位時間數(shù)據(jù)。
6.如權(quán)利要求1至5中的任意一項所定義的方法，其特征還在于以下步驟 將兩個導(dǎo)出的平均自相關(guān)用作預(yù)測去卷積算法的輸入而獲得兩個因果反向濾波器以應(yīng)用于所述疊后數(shù)據(jù)體中的多個地震道中的每一個，所述兩個平均自相關(guān)各自相對于不 同的數(shù)據(jù)集，所述數(shù)據(jù)集包括(1)共同縱測線號和共同水底時間，以及(2)所述共同縱測線號。
7.如權(quán)利要求1至6中的任意一項所定義的方法，其特征還在于以下步驟 測量在所述雙程水底時間處所述多個地震道中的每一個的自相關(guān)幅度水平以定義水底時間幅度；以及在計算機顯示器上顯示所述多個地震道中的每一個的水底時間幅度的質(zhì)量控制圖。
8.如權(quán)利要求1至7中的任意一項所定義的方法，其特征還在于以下步驟 響應(yīng)于所述多個導(dǎo)出的平均自相關(guān)中的至少一個而為所述多個不同的域中的至少一個確定導(dǎo)出的平均自相關(guān)在對應(yīng)于所述雙程水底時間的時間樣點處的幅度值；以及 顯示在雙程水底時間處的導(dǎo)出的平均自相關(guān)幅度值的質(zhì)量控制圖。
9.一種抑制地震數(shù)據(jù)中的殘余水底能量的用計算機實現(xiàn)的方法，所述方法包括獲得 地震的疊后數(shù)據(jù)的步驟，所述方法的特征還在于將多域去卷積應(yīng)用于所述地震的疊后數(shù)據(jù)以由此標(biāo)識存在于所述地震數(shù)據(jù)中的殘余 水底能量；以及響應(yīng)于對存在于所述地震數(shù)據(jù)中的殘余水底能量的標(biāo)識而抑制地震數(shù)據(jù)體中的殘余 水底能量。
10.如權(quán)利要求9所定義的方法，其中所述將多域去卷積應(yīng)用于所述地震的疊后數(shù) 據(jù)的步驟包括將共同水底時間與縱測線號的結(jié)合用作所述多域去卷積的自相關(guān)求平均基 準(zhǔn)。
11.如權(quán)利要求9或10中的任意一項所定義的方法，其中所述地震的疊后數(shù)據(jù)包括三 維疊加數(shù)據(jù)體，該三維疊加數(shù)據(jù)體包括多個地震道，并且其中所述將多域去卷積應(yīng)用于 所述地震的疊后數(shù)據(jù)的步驟包括以下步驟訪問所述疊加的地震數(shù)據(jù)； 訪問每個地震道的雙程水底時間； 訪問每個地震道的地表下方層位時間；拉平所述地震數(shù)據(jù)以形成對于所述體中的每個地震道都處于恒定時間的基準(zhǔn)層位；以及從拉平的地震體中提取時間窗。
12.如權(quán)利要求11所定義的方法，其中所述將多域去卷積應(yīng)用于所述地震的疊后數(shù)據(jù) 的步驟還包括以下步驟根據(jù)被提取的時間窗計算所述地震數(shù)據(jù)體中的多個地震道中的每一個的自相關(guān)； 測量在所述雙程水底時間處每個地震道上的自相關(guān)幅度水平；以及 將所述幅度水平值繪制成質(zhì)量控制圖。
13.如權(quán)利要求11或12中的任意一項所定義的方法，其中所述將多域去卷積應(yīng)用于 所述地震的疊后數(shù)據(jù)的步驟還包括以下步驟響應(yīng)于被提取的時間窗而計算所述地震數(shù)據(jù)體中的多個地震道中的每一個的自相 關(guān)；以及在多個域中對計算出的自相關(guān)求平均以由此標(biāo)識存在于所述地震數(shù)據(jù)中的水底能量。
14.如權(quán)利要求13所定義的方法，其中所述疊加的地震數(shù)據(jù)被存儲在數(shù)據(jù)記錄中，每 個所述數(shù)據(jù)記錄具有包含橫測線和縱測線號的首部，并且其中所述在多個域中對計算出 的自相關(guān)求平均的步驟包括以下步驟計算具有共同縱測線和共同水底時間的第一平均自相關(guān)； 計算具有共同縱測線號的第二平均自相關(guān)；以及 計算具有共同橫測線號的第三平均自相關(guān)。
15.如權(quán)利要求14所定義的方法，其中所述將多域去卷積應(yīng)用于地震的疊后數(shù)據(jù)的步 驟還包括以下步驟響應(yīng)于計算出的第一平均自相關(guān)而為每個縱測線號將所述地震數(shù)據(jù)疊加到共同水底 時間面元中；顯示每個縱測線號的共同水底時間的自相關(guān)疊加；為每個縱測線號確定所述地震數(shù)據(jù)在共同水底時間處的幅度值；以及顯示所述地震數(shù)據(jù)對于每個縱測線號在共同水底時間處的幅度值。
16.如權(quán)利要求14或15中的任意一項所定義的方法，其中所述將多域去卷積應(yīng)用于 地震的疊后數(shù)據(jù)的步驟包括以下步驟應(yīng)用預(yù)測去卷積算法以響應(yīng)于所述第一和所述第二平均自相關(guān)而獲得定義多個算子 的多個因果反向濾波器以應(yīng)用于所述地震疊加體中的多個地震道中的每一個；以及將所述多個算子應(yīng)用于所述多個地震道中的每一個以由此抑制地震數(shù)據(jù)體中的殘余 水底能量。
17.如權(quán)利要求16所定義的方法，其中所述將多域去卷積應(yīng)用于地震的疊后數(shù)據(jù)的步 驟還包括以下步驟響應(yīng)于計算出的第一平均自相關(guān)而為每個縱測線號將所述地震數(shù)據(jù)疊加到共同水底 時間面元中；為每個縱測線號確定所述地震數(shù)據(jù)在共同水底時間處的幅度值； 顯示所述地震數(shù)據(jù)對于每個縱測線號在共同水底時間處的幅度值； 響應(yīng)于所顯示的幅度值而確定對在空間上可變的去卷積算子的需要； 響應(yīng)于計算出的第一平均自相關(guān)而為每個縱測線號將所述地震數(shù)據(jù)疊加到共同水底 時間面元中；顯示每個縱測線號的共同水底時間的自相關(guān)疊加；以及響應(yīng)于所顯示的每個縱測線號的共同水底時間的自相關(guān)疊加而確定所述多個算子的 長度和間隙參數(shù)。
18.如權(quán)利要求9至17中的任意一項所定義的方法，其中所述應(yīng)用多域去卷積的步驟 還包括以下步驟將地震數(shù)據(jù)疊加到共深度點面元中以形成所述三維疊加數(shù)據(jù)體，所述疊加的地震數(shù) 據(jù)包括多個地震道；在所述多個疊加的地震道中的每一個上執(zhí)行多個自相關(guān)； 在多個域中對所述相關(guān)數(shù)據(jù)求平均；響應(yīng)于所述多個域中的平均相關(guān)數(shù)據(jù)而確定多個預(yù)測去卷積逆算子；將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)；以及響應(yīng)于將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)而標(biāo)識存在于所述地震數(shù)據(jù)中 的水底能量。
19.一種用于抑制地震數(shù)據(jù)中的殘余水底能量檢測的地震數(shù)據(jù)分析程序產(chǎn)品(13)，其 被存儲在可由計算機(14)讀取的有形的計算機可讀介質(zhì)(16)上，所述程序產(chǎn)品(13)包 括指令集，該指令集在被所述計算機(14)執(zhí)行時使所述計算機(14)執(zhí)行以下操作獲 得定義疊后數(shù)據(jù)體的共深度點疊加的三維地震數(shù)據(jù)，所述疊后數(shù)據(jù)體包括多個地震道， 每個所述地震道包括首部，該首部包括橫測線和縱測線號；獲得所述多個地震道中的每 一個的雙程水底時間數(shù)據(jù)；以及獲得所述多個地震道中的每一個的地表下方層位時間數(shù) 據(jù)，所述操作的特征在于拉平所述疊加的地震數(shù)據(jù)以定義拉平的地震數(shù)據(jù)體使得基準(zhǔn)層位對于所述疊后數(shù)據(jù) 體中的多個道中的每一個都處于恒定時間； 從拉平的地震數(shù)據(jù)體中提取時間窗；響應(yīng)于所述體中相應(yīng)的地震道的相應(yīng)的被提取的時間窗而計算所述體中的多個地震 道中的每個單獨的地震道的自相關(guān)；響應(yīng)于所述多個自相關(guān)而在多個不同的域上執(zhí)行多域自相關(guān)分析；以及 使用來自所述多個不同的域中的至少兩個的定義分選鍵的項將多域去卷積應(yīng)用于所 述地震的疊后數(shù)據(jù)以由此標(biāo)識和抑制所述殘余水底能量。
20.如權(quán)利要求19所定義的程序產(chǎn)品(13)，其中所述執(zhí)行多域自相關(guān)分析的操作包括響應(yīng)于所述多個自相關(guān)而導(dǎo)出多個平均自相關(guān)，所述導(dǎo)出包括在多個不同的域的每 一個中為所述多個自相關(guān)中的每一個確定平均值。
21.如權(quán)利要求20所定義的程序產(chǎn)品(13)，其中所述在多個不同的域的每一個中為 所述多個自相關(guān)中的每一個確定平均值的操作包括響應(yīng)于所述多個地震道中的每一個的首部數(shù)據(jù)而輸出三個平均數(shù)據(jù)集，所述數(shù)據(jù)集 包括(1)共同縱測線號和共同水底時間，(2)所述共同縱測線號，以及(3)共同橫測線 號。
22.如權(quán)利要求20或21中的任意一項所定義的程序產(chǎn)品(13)，所述操作的特征還在于響應(yīng)于所述多個導(dǎo)出的平均自相關(guān)中的至少一個而將所述疊后地震數(shù)據(jù)疊加到共同 水底時間面元中；以及響應(yīng)于所述將所述疊后地震數(shù)據(jù)疊加到共同水底時間面元中的操作而顯示按道格式 的疊后地震數(shù)據(jù)的質(zhì)量控制圖。
23.如權(quán)利要求22所定義的程序產(chǎn)品(13)，其中所述將多域去卷積應(yīng)用于所述地震 的疊后數(shù)據(jù)的操作包括響應(yīng)于按道格式的疊后地震數(shù)據(jù)的質(zhì)量控制圖而接收去卷積算子的參數(shù)； 響應(yīng)于所述多個導(dǎo)出的平均自相關(guān)中的至少兩個和所述去卷積算子的參數(shù)而獲得確 定一個或多個去卷積算子的兩個因果反向濾波器以應(yīng)用于所述疊后數(shù)據(jù)體中的多個地震 道中的每一個；以及將所述去卷積算子應(yīng)用于所述多個地震道中的每一個的地表下方層位時間數(shù)據(jù)。
24.如權(quán)利要求19至23中的任意一項所定義的程序產(chǎn)品(13)，所述操作的特征還在于將兩個導(dǎo)出的平均自相關(guān)用作預(yù)測去卷積算法的輸入而獲得兩個因果反向濾波器以 應(yīng)用于所述疊后數(shù)據(jù)體中的多個地震道中的每一個，所述兩個平均自相關(guān)各自相對于不 同的數(shù)據(jù)集，所述數(shù)據(jù)集包括(1)共同縱測線號和共同水底時間，以及(2)所述共同縱測 線號。
25.如權(quán)利要求19至24中的任意一項所定義的程序產(chǎn)品(13)，所述操作的特征還在于測量在所述雙程水底時間處所述多個地震道中的每一個的自相關(guān)幅度水平以定義水 底時間幅度；以及在計算機顯示器上顯示所述多個地震道中的每一個的水底時間幅度的質(zhì)量控制圖。
26.如權(quán)利要求19至25中的任意一項所定義的程序產(chǎn)品(13)，所述操作的特征還在于響應(yīng)于所述多個導(dǎo)出的平均自相關(guān)中的至少一個而為所述多個不同的域中的至少一 個確定導(dǎo)出的平均自相關(guān)在對應(yīng)于所述雙程水底時間的時間樣點處的幅度值；以及 顯示在雙程水底時間處的導(dǎo)出的平均自相關(guān)幅度值的質(zhì)量控制圖。
27.—種用于抑制地震數(shù)據(jù)中的殘余水底能量檢測的地震數(shù)據(jù)分析程序產(chǎn)品(13)，其 被存儲在可由計算機(14)讀取的有形的計算機可讀介質(zhì)(16)上，所述程序產(chǎn)品(13)包 括指令集，該指令集在被所述計算機(14)執(zhí)行時使所述計算機(14)執(zhí)行包括獲得地震的 疊后數(shù)據(jù)的操作，所述操作的特征在于將多域去卷積應(yīng)用于所述地震的疊后數(shù)據(jù)以由此標(biāo)識存在于所述地震數(shù)據(jù)中的殘余 水底能量；以及響應(yīng)于對存在于所述地震數(shù)據(jù)中的殘余水底能量的標(biāo)識而抑制地震數(shù)據(jù)體中的殘余 水底能量。
28.如權(quán)利要求27所定義的程序產(chǎn)品(13)，其中所述將多域去卷積應(yīng)用于所述地震 的疊后數(shù)據(jù)的操作包括將共同水底時間與縱測線號的結(jié)合用作所述多域去卷積的自相關(guān) 求平均基準(zhǔn)。
29.如權(quán)利要求27或28中的任意一項所定義的程序產(chǎn)品(13)，其中所述地震的疊后 數(shù)據(jù)包括三維疊加數(shù)據(jù)體，該三維疊加數(shù)據(jù)體包括多個地震道，并且其中所述將多域去 卷積應(yīng)用于所述地震的疊后數(shù)據(jù)的操作包括以下操作訪問所述疊加的地震數(shù)據(jù)； 訪問每個地震道的雙程水底時間； 訪問每個地震道的地表下方層位時間；拉平所述地震數(shù)據(jù)以形成對于所述體中的每個地震道都處于恒定時間的基準(zhǔn)層位；以及從拉平的地震體中提取時間窗。
30.如權(quán)利要求29所定義的程序產(chǎn)品(13)，其中所述將多域去卷積應(yīng)用于所述地震 的疊后數(shù)據(jù)的操作還包括以下操作根據(jù)被提取的時間窗計算所述地震數(shù)據(jù)體中的多個地震道中的每一個的自相關(guān)； 測量在所述雙程水底時間處每個地震道上的自相關(guān)幅度水平；以及 將所述幅度水平值繪制成質(zhì)量控制圖。
31.如權(quán)利要求29或30中的任意一項所定義的程序產(chǎn)品(13)，其中所述將多域去卷 積應(yīng)用于所述地震的疊后數(shù)據(jù)的操作還包括以下操作響應(yīng)于被提取的時間窗而計算所述地震數(shù)據(jù)體中的多個地震道中的每一個的自相 關(guān)；以及在多個域中對計算出的自相關(guān)求平均以由此標(biāo)識存在于所述地震數(shù)據(jù)中的水底能量。
32.如權(quán)利要求31所定義的程序產(chǎn)品(13)，其中所述疊加的地震數(shù)據(jù)被存儲在數(shù)據(jù) 記錄中，每個所述數(shù)據(jù)記錄具有包含橫測線和縱測線號的首部，并且其中所述在多個域 中對計算出的自相關(guān)求平均的操作包括以下操作計算具有共同縱測線和共同水底時間的第一平均自相關(guān)； 計算具有共同縱測線號的第二平均自相關(guān)；以及 計算具有共同橫測線號的第三平均自相關(guān)。
33.如權(quán)利要求32所定義的程序產(chǎn)品(13)，其中所述將多域去卷積應(yīng)用于地震的疊 后數(shù)據(jù)的操作還包括以下操作響應(yīng)于計算出的第一平均自相關(guān)為每個縱測線號將所述地震數(shù)據(jù)疊加到共同水底時 間面元中；顯示每個縱測線號的共同水底時間的自相關(guān)疊加；為每個縱測線號確定所述地震數(shù)據(jù)在共同水底時間處的幅度值；以及顯示所述地震數(shù)據(jù)對于每個縱測線號在共同水底時間處的幅度值。
34.如權(quán)利要求32或33中的任意一項所定義的程序產(chǎn)品(13)，其中所述將多域去卷 積應(yīng)用于地震的疊后數(shù)據(jù)的操作還包括以下操作應(yīng)用預(yù)測去卷積算法以響應(yīng)于所述第一和所述第二平均自相關(guān)而獲得定義多個算子 的多個因果反向濾波器以應(yīng)用于所述地震疊加體中的多個地震道中的每一個；以及將所述多個算子應(yīng)用于所述多個地震道中的每一個以由此抑制地震數(shù)據(jù)體中的殘余 水底能量。
35.如權(quán)利要求34所定義的程序產(chǎn)品(13)，其中所述將多域去卷積應(yīng)用于地震的疊 后數(shù)據(jù)的操作還包括以下操作響應(yīng)于計算出的第一平均自相關(guān)而為每個縱測線號將所述地震數(shù)據(jù)疊加到共同水底 時間面元中；為每個縱測線號確定所述地震數(shù)據(jù)在共同水底時間處的幅度值； 顯示所述地震數(shù)據(jù)對于每個縱測線號在共同水底時間處的幅度值； 響應(yīng)于所顯示的幅度值而確定對在空間上可變的去卷積算子的需要； 響應(yīng)于計算出的第一平均自相關(guān)而為每個縱測線號將所述地震數(shù)據(jù)疊加到共同水底 時間面元中；顯示每個縱測線號的共同水底時間的自相關(guān)疊加；以及響應(yīng)于所顯示的每個縱測線號的共同水底時間的自相關(guān)疊加而確定所述多個算子的長度和間隙參數(shù)。
36.如權(quán)利要求27至35中的任意一項所定義的程序產(chǎn)品(13)，其中所述應(yīng)用多域去 卷積的操作還包括以下操作將地震數(shù)據(jù)疊加到共深度點面元中以形成所述三維疊加數(shù)據(jù)體，所述疊加的地震數(shù) 據(jù)包括多個地震道；在所述多個疊加的地震道中的每一個上執(zhí)行多個自相關(guān)； 在多個域中對所述相關(guān)數(shù)據(jù)求平均；響應(yīng)于所述多個域中的平均相關(guān)數(shù)據(jù)而確定多個預(yù)測去卷積逆算子； 將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)；以及響應(yīng)于將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)而標(biāo)識存在于所述地震數(shù)據(jù)中 的水底能量。
37.—種可由計算機(14)讀取以分析地震數(shù)據(jù)的計算機可讀介質(zhì)，所述計算機可讀 介質(zhì)包括用于抑制地震數(shù)據(jù)中的殘余水底能量檢測的指令集，該指令集在被所述計算機 (14)執(zhí)行時使所述計算機(14)執(zhí)行包括接收地震的疊后數(shù)據(jù)的操作，所述操作的特征在 于將多域去卷積應(yīng)用于被接收的地震的疊后數(shù)據(jù)以由此抑制地震數(shù)據(jù)體中的殘余水底能量°
38.如權(quán)利要求37所定義的計算機可讀介質(zhì)，其中所述將多域去卷積應(yīng)用于被接收的 地震的疊后數(shù)據(jù)的操作包括將共同水底時間與縱測線號的結(jié)合用作所述多域去卷積的自 相關(guān)求平均基準(zhǔn)。
39.如權(quán)利要求37或38中的任意一項所定義的計算機可讀介質(zhì)，其中所述應(yīng)用多域 去卷積的操作還包括以下操作將地震數(shù)據(jù)疊加到共深度點面元中以形成三維疊加數(shù)據(jù)體，所述疊加的地震數(shù)據(jù)包 括多個地震道；在所述多個疊加的地震道中的每一個上執(zhí)行多個自相關(guān)； 在多個域中對所述相關(guān)數(shù)據(jù)求平均；響應(yīng)于所述多個域中的平均相關(guān)數(shù)據(jù)而確定多個預(yù)測去卷積逆算子； 將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)；以及響應(yīng)于將所述預(yù)測去卷積逆算子應(yīng)用于所述地震數(shù)據(jù)而標(biāo)識存在于所述地震數(shù)據(jù)中 的水底能量。
40.如權(quán)利要求37至39中的任意一項所定義的計算機可讀介質(zhì)，其中所述疊加的地 震數(shù)據(jù)包括多個地震道，并且其中所述操作的特征還在于計算所述地震數(shù)據(jù)體中的多個地震道中的每個單獨的地震道的自相關(guān)； 響應(yīng)于所述多個自相關(guān)而在多個不同的域上執(zhí)行多域自相關(guān)分析，以及 使用來自所述多個不同域中的至少兩個的定義分選鍵的項將所述多域去卷積應(yīng)用于 所述地震的疊后數(shù)據(jù)以由此標(biāo)識和抑制所述殘余水底能量。
全文摘要
本發(fā)明提供了抑制地震數(shù)據(jù)中的殘余水底能量的系統(tǒng)(12)、程序產(chǎn)品(13)和方法。系統(tǒng)(12)、程序產(chǎn)品(13)和方法的示例可以被應(yīng)用于疊后數(shù)據(jù)集并且可以將多通道去卷積與新穎的分選鍵結(jié)合以有效地標(biāo)識和抑制共深度點(CDP)疊加的地震數(shù)據(jù)中的殘余水底能量，由此提高地震數(shù)據(jù)的分辨能力，從而得到對從儲油層反射的地震信號的改進的判讀。
文檔編號G01V1/38GK102016642SQ200980115887
公開日2011年4月13日 申請日期2009年2月27日 優(yōu)先權(quán)日2008年2月29日
發(fā)明者M·E·赫德法, R·M·伯恩斯塔 申請人:沙特阿拉伯石油公司