利用疊前地震道集計算地層品質因數的方法及裝置製造方法
2023-09-21 03:35:30
利用疊前地震道集計算地層品質因數的方法及裝置製造方法
【專利摘要】本發明實施例提供了一種利用疊前地震道集計算地層品質因數的方法,包括:獲取疊前地震數據,並進行預處理;對疊前地震數據進行疊加速度分析和解釋,生成疊加速度;從疊加剖面上識別和拾取反射層,利用所拾取的反射層時間和疊加速度提取反射波,得到反射波隨炮檢距變化的波形,生成反射波波形;對反射波波形進行傅立葉變換,計算反射波振幅譜;根據反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量δn;根據層對數衰減量,利用公式計算得到地層品質因數Qn。本發明實施例的利用疊前地震道集計算地層品質因數的方法與裝置,對震源子波的振幅譜和相位譜未作任何假設,具有計算效率高、精度較高以及適應性強的優點。
【專利說明】利用疊前地震道集計算地層品質因數的方法及裝置
【技術領域】
[0001]本發明涉及地震數據處理領域,尤其涉及一種利用疊前地震道集計算地層品質因數的方法及裝置。
【背景技術】
[0002]地下介質廣泛存在衰減特性,地層的吸收衰減主要表現為地震波在傳播過程中振幅發生衰減、相位發生畸變、頻率變低,且高頻部分比低頻部分衰減更快,在淺層比深層衰減更快,導致地震波主頻降低、頻帶變窄,降低了地震資料的解析度,從而影響到地震勘探的精度。
[0003]為了提高地震資料的解析度,可以採用反褶積方法壓縮子波長度,可以補償地層對地震波的吸收衰減。吸收衰減的補償方法主要包括反Q濾波方法和基於衰減補償的偏移方法,其關鍵之一在於Q值的求取。
[0004]Rainner Tonn (1989)對早期計算Q值的七種計算方法進行了比較,時間域方法有:振幅衰減法、上升時間法、子波模擬法、解析信號法;頻率域的方法有:譜比法、匹配法、譜模擬法。他得出的結論是:沒有哪一種方法適用於任何情況,它們的效果依賴於記錄的質量。在這些方法中,時間域的上升時間法和頻率域的譜比法最為引起人們的關注。
[0005]早期計算Q值方法主要基於疊後數據或零偏VSP數據。相對於這些數據而言,由於疊前地震數據包含更多信息,因此利用疊前數據進行反演能獲得更高精度的Q值。Dasgupta等(1998)給出了基於CMP道集的計算Q值流程。Zhang等(2002)給出了基於水平層狀介質模型的CMP道集反射波衰減方程,並利用這個方程進行Q值反演和反Q濾波。由於這個方程是基於直射線傳播路徑假設推導出來的,因此影響了其精度。為了提高CMP道集反射波吸收方程精度,劉洋等(2009)提出了基於曲射線路徑和斯奈爾定理,推導出了更高精度的吸收方程,並發展了基於質心頻率的Q值反演方法。
[0006]目前利用疊前道集估算地層品質因數研究不多,均存在一些假設,對震源子波進行了一定的假設,因而限制了現有的方法的廣泛應用。
【發明內容】
[0007]本發明的目的是提供一種利用疊前反射波計算層品質因數的方法,對震源子波的振幅譜和相位譜未做任何假設,具有易於實現、精度較高、適用性強的優點。
[0008]為了達到上述目的,本發明實施例提供一種利用疊前地震道集計算地層品質因數的方法,包括:獲取疊前地震數據,並進行預處理;對所述疊前地震數據進行疊加速度分析和解釋,生成疊加速度;從疊加剖面上識別和拾取反射層,利用所拾取的反射層時間和所述疊加速度提取反射波,得到所述反射波隨炮檢距變化的波形,生成反射波波形;對所述反射波波形進行傅立葉變換,計算反射波振幅譜;根據所述的反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量δ η ;根據所述層對數衰減量,利用公式計算得到地層品質因數Qn。
[0009]進一步地,在一實施例中,所述獲取疊前地震數據,並進行預處理,包括:對所述疊前地震數據進行觀測系統加載、道編輯和保幅去噪處理。
[0010]進一步地,在一實施例中對所述疊前地震數據進行疊加速度分析和解釋,生成疊加速度,包括:對所述疊前地震數據進行疊加速度分析,生成疊加速度譜;對所述疊加速度譜進行解釋,生成疊加速度。
[0011]進一步地,在一實施例中在從疊加剖面上識別和拾取反射層之前,還包括:對所述疊前地震數據進行動校正和水平疊加,獲得所述疊加剖面。
[0012]進一步地,在一實施例中根據所述的反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量Sn,包括:對於疊前共中心點道集,所述反射
波振幅譜為
【權利要求】
1.一種利用疊前地震道集計算地層品質因數的方法,其特徵在於,所述方法包括: 獲取疊前地震數據,並進行預處理; 對所述疊前地震數據進行疊加速度分析和解釋,生成疊加速度; 從疊加剖面上識別和拾取反射層,利用所拾取的反射層時間和所述疊加速度提取反射波,得到所述反射波隨炮檢距變化的波形,生成反射波波形; 對所述反射波波形進行傅立葉變換,計算反射波振幅譜; 根據所述的反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量δ η ; 根據所述層對數衰減量,利用公式
2.根據權利要求1所述的利用疊前地震道集計算地層品質因數的方法,其特徵在於,所述獲取疊前地震數據,並進行預處理,包括: 對所述疊前地震數據進行觀測系統加載、道編輯和保幅去噪處理。
3.根據權利要求1所述的利用疊前地震道集計算地層品質因數的方法,其特徵在於,對所述疊前地震數據進行疊加速度分析和解釋,生成疊加速度,包括: 對所述疊前地震數據進行疊加速度分析,生成疊加速度譜; 對所述疊加速度譜進行解釋,生成疊加速度。
4.根據權利要求1所述的利用疊前地震道集計算地層品質因數的方法,其特徵在於,在從疊加剖面上識別和拾取反射層之前,還包括: 對所述疊前地震數據進行動校正和水平疊加,獲得所述疊加剖面。
5.根據權利要求4所述的利用疊前地震道集計算地層品質因數的方法,其特徵在於,根據所述的反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量δ η,包括: 對於疊前共中心點道集,所述反射波振幅譜為
6.一種利用疊前地震道集計算地層品質因數的裝置,其特徵在於,所述裝置包括: 預處理單元,用於獲取疊前地震數據,並進行預處理; 速度分析解釋單元,用於對所述疊前地震數據進行疊加速度分析和解釋,生成疊加速度; 反射波波形生成 單元,用於從疊加剖面上識別和拾取反射層,利用所拾取的反射層時間和所述疊加速度提取反射波,得到所述反射波隨炮檢距變化的波形,生成反射波波形;振幅譜生成單元,用於對所述反射波波形進行傅立葉變換,計算反射波振幅譜; 層對數衰減量計算單元,用於根據所述的反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量δ η ; 地層品質因數計算單元,用於根據所述層對數衰減量,利用公式
7.根據權利要求6所述的利用疊前地震道集計算地層品質因數的裝置,其特徵在於,所述預處理單元用於對所述疊前地震數據進行觀測系統加載、道編輯和保幅去噪處理。
8.根據權利要求6所述的利用疊前地震道集計算地層品質因數的裝置,其特徵在於,所述速度分析解釋單元包括: 速度分析單元,用於對所述疊前地震數據進行疊加速度分析,生成疊加速度譜; 速度譜解釋單元,用於對所述疊加速度譜進行解釋,生成疊加速度。
9.根據權利要求6所述的利用疊前地震道集計算地層品質因數的裝置,其特徵在於,所述裝置還包括: 疊加剖面獲取單元,用於對所述疊前地震數據進行動校正和水平疊加,獲得所述疊加剖面。
10.根據權利要求9所述的利用疊前地震道集計算地層品質因數的裝置,其特徵在於,所述層對數衰減量計算單元根據所述的反射波振幅譜,利用所述疊加速度、反射層時間以及濾波器振幅譜,計算得到層對數衰減量Sn,包括: 對於疊前共中心點道集,所述反射波振幅譜為.
【文檔編號】G01V1/28GK103698812SQ201310731132
【公開日】2014年4月2日 申請日期:2013年12月26日 優先權日:2013年12月26日
【發明者】劉洋 申請人:中國石油天然氣集團公司, 中國石油大學(北京)