• / 71
  • 下載費用:30 金幣  

疊后地震數據體處理方法及裝置.pdf

摘要
申請專利號:

CN201610058962.1

申請日:

2016.01.28

公開號:

CN106257309A

公開日:

2016.12.28

當前法律狀態:

授權

有效性:

有權

法律詳情: 授權|||實質審查的生效IPC(主分類):G01V 1/30申請日:20160128|||公開
IPC分類號: G01V1/30 主分類號: G01V1/30
申請人: 中國石油天然氣股份有限公司
發明人: 吳吉忠; 穆立華; 趙寶銀; 石文武; 邢文軍; 楊曉利; 郭愛華; 張建坤
地址: 100007 北京市東城區東直門北大街9號中國石油大廈
優先權:
專利代理機構: 北京三高永信知識產權代理有限責任公司 11138 代理人: 周靜
PDF完整版下載: PDF下載
法律狀態
申請(專利)號:

CN201610058962.1

授權公告號:

||||||

法律狀態公告日:

2018.11.16|||2017.01.25|||2016.12.28

法律狀態類型:

授權|||實質審查的生效|||公開

摘要

本發明公開一種疊后地震數據體處理方法及裝置,屬于地震數據處理領域。方法包括:從疊后地震數據體對應的地震測線中選取M條樣本測線,根據M條樣本測線對應的校正Q值場和M條樣本測線對應的加權系數場確定M條樣本測線對應的最終Q值場,根據M條樣本測線對應的樣本速度場和最終Q值場更新李慶忠經驗公式,采用更新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數據體以便于GPU根據地震工區的Q值場對疊后地震數據體進行處理。本發明解決了疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體處理的靈活性和可靠性的效果。本發明用于疊后地震數據體的處理。

權利要求書

1.一種疊后地震數據體處理方法,其特征在于,用于中央處理器CPU,所述方法包括:
從地震工區的疊后地震數據體對應的地震測線中選取M條樣本測線,所述M為大于或者
等于2的整數,且M=0.003N,所述N為所述疊后地震數據體的地震主測線的條數,每條所述
地震測線對應一個地震剖面;
從所述地震工區的疊加速度場中確定所述M條樣本測線中的每條樣本測線對應的速度
場,得到所述M條樣本測線對應的樣本速度場ν;
根據所述M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公式確定所述M條樣本測
線對應的初始Q值場Q0,所述李慶忠經驗公式為Q=α·νβ,α=14,β=2.2;
對所述M條樣本測線對應的初始Q值場Q0進行校正,得到所述M條樣本測線對應的校正Q
值場Q1;
根據所述M條樣本測線對應的校正Q值場Q1和所述M條樣本測線對應的加權系數場
λfinal,確定所述M條樣本測線對應的最終Q值場Q2;
根據所述M條樣本測線對應的樣本速度場ν和所述M條樣本測線對應的最終Q值場Q2,采
用最小二乘法更新所述李慶忠經驗公式,得到更新后的李慶忠經驗公式;
根據所述地震工區的疊加速度場,采用所述更新后的李慶忠經驗公式確定所述地震工
區的Q值場Qall;
向圖形處理器GPU發送所述地震工區的Q值場Qall和所述疊后地震數據體,以便于所述
GPU根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處理,得到處理后的
疊后地震數據體。
2.根據權利要求1所述的方法,其特征在于,在所述根據所述M條樣本測線對應的校正Q
值場Q1和所述M條樣本測線對應的加權系數場λfinal,確定所述M條樣本測線對應的最終Q值
場Q2之前,所述方法還包括:
根據預設的加權系數組對所述M條樣本測線對應的校正Q值場Q1進行處理,得到所述M條
樣本測線對應的加權Q值場組,所述預設的加權系數組中包括P個加權系數λ,所述加權Q值
場組中包括P套加權Q值場,所述P為大于或者等于1的整數;
根據所述加權Q值場組對所述M條樣本測線組成的樣本地震數據體進行反Q濾波處理,
得到反Q濾波數據體組,所述反Q濾波數據體組中包括P套反Q濾波數據體,每套反Q濾波數據
體包括M條反Q樣本測線,每條所述反Q樣本測線對應一個地震剖面;
采用傅里葉變換對所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測
線對應的地震剖面進行頻譜分析;
根據頻譜分析結果進行處理,確定所述M條反Q樣本測線對應的加權系數場;
將所述M條反Q樣本測線對應的加權系數場確定為所述M條樣本測線對應的加權系數場
λfinal。
3.根據權利要求2所述的方法,其特征在于,
所述采用傅里葉變換對所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的地震剖面進行頻譜分析,包括:
在所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地震剖
面上選擇三個矩形時窗,所述三個矩形時窗包括:淺層矩形時窗、中層矩形時窗和深層矩形
時窗;
采用傅里葉變換對所述每條反Q樣本測線對應的地震剖面上位于所述三個矩形時窗中
的每個矩形時窗內的部分進行處理,得到所述每條反Q樣本測線對應的3個反Q樣本頻譜;
根據所述每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻譜,得到所述每套
反Q濾波數據體對應的3M個反Q樣本頻譜;
根據所述反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻譜得到3×M×P個
反Q樣本頻譜,所述3×M×P個反Q樣本頻譜中包括:M×P個淺層矩形時窗對應的反Q樣本頻
譜、M×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層矩形時窗對應的反Q樣本頻譜;
對所述3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。
4.根據權利要求3所述的方法,其特征在于,所述根據頻譜分析結果進行處理,確定所
述M條反Q樣本測線對應的加權系數場,包括:
根據頻譜分析結果確定所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,所述第一目標頻譜為所述每
條反Q樣本測線對應的地震剖面上位于所述淺層矩形時窗內的部分對應的P個反Q樣本頻譜
中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第二目標頻譜為所述每條反Q樣
本測線對應的地震剖面上位于所述中層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫
軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第三目標頻譜為所述每條反Q樣本測線
對應的地震剖面上位于所述深層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成
的封閉區域的面積最大的反Q樣本頻譜;
分別確定所述第一目標頻譜對應的矩形時窗對應的第一目標加權系數、所述第二目標
頻譜對應的矩形時窗對應的第二目標加權系數和所述第三目標頻譜對應的矩形時窗對應
的第三目標加權系數,得到所述每條反Q樣本測線對應的3個目標加權系數,其中,所述第一
目標加權系數、所述第二目標加權系數和所述第三目標加權系數都為所述預設的加權系數
組中的加權系數;
根據所述M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權系數,得到所述M
條反Q樣本測線對應的3M個目標加權系數;
對所述3M個目標加權系數進行空間插值,得到所述地震工區對應的加權系數場;
對所述地震工區對應的加權系數場進行平滑處理,得到處理后的加權系數場;
從所述處理后的加權系數場中抽取出所述M條反Q樣本測線對應的加權系數場。
5.根據權利要求1至4任一所述的方法,其特征在于,
在所述向圖形處理器GPU發送所述地震工區的Q值場Qall和所述疊后地震數據體之前,
所述方法還包括:
按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣,得到采樣Q值場Qs;
對所述疊后地震數據體進行分塊處理,得到D個地震數據塊,所述D為大于或者等于2的
整數;
所述向圖形處理器GPU發送所述地震工區的Q值場Qall和疊后地震數據體,包括:
依次向所述GPU發送所述采樣Q值場Qs和所述D個地震數據塊中的每個地震數據塊。
6.根據權利要求5所述的方法,其特征在于,所述對所述疊后地震數據體進行分塊處
理,得到D個地震數據塊,包括:
計算所述疊后地震數據體中的地震數據的總道數E;
根據所述GPU在一個周期內處理的地震數據的道數Tracemax和所述疊后地震數據體中
的地震數據的總道數E,對所述疊后地震數據體進行分塊處理,得到所述D個地震數據塊;
其中,E=Tracemax×(D-1)+L,所述L為所述D個地震數據塊中的最后一個地震數據塊中
的地震數據的道數。
7.一種疊后地震數據體處理方法,其特征在于,用于圖形處理器GPU,所述方法包括:
接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據體;
根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處理,得到處理后
的疊后地震數據體。
8.根據權利要求7所述的方法,其特征在于,
所述接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據體,包括:
接收所述CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,所述采樣Q值
場Qs是所述CPU按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣得到的,所述D個地
震數據塊是所述CPU對所述疊后地震數據體進行分塊處理得到的;
所述根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處理,得到處
理后的疊后地震數據體,包括:
根據所述采樣Q值場Qs對所述D個地震數據塊中的每個地震數據塊進行反Q濾波處理,得
到D個處理后的地震數據塊;
根據所述D個處理后的地震數據塊,得到所述處理后的疊后地震數據體。
9.根據權利要求8所述的方法,其特征在于,
所述根據所述采樣Q值場Qs對所述D個地震數據塊中的每個地震數據塊進行反Q濾波處
理,得到D個處理后的地震數據塊,包括:
對所述采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ;
根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震數據塊中的所有時
間深度τ處的地震幅值進行反Q濾波處理,得到所述D個處理后的地震數據塊。
10.根據權利要求9所述的方法,其特征在于,所述對所述采樣Q值場Qs進行處理得到每
個時間深度τ處的Q值Qτ,包括:
根據所述采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ;
其中,所述Q值插值計算公式為:
<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mi>=</mi> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>所述τj、所述τ(j-1)和所述τ(j+1)都為時
間深度,所述為時間深度τj處的Q值插值,所述為時間深度τ(j-1)處的Q值插值,所述
為時間深度τ(j+1)處的Q值插值,表示對所述取整。
11.根據權利要求9或10所述的方法,其特征在于,所述根據所有時間深度τ處的Q值Qτ對
所述D個地震數據塊中的每個地震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處
理,得到所述D個處理后的地震數據塊,包括:
根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震數據塊中的每個地
震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理后的地震數據道;
根據所述每個地震數據塊中的所有處理后的地震數據道,得到處理后的地震數據塊;
根據所述D個地震數據塊中的所有處理后的地震數據塊,得到所述D個處理后的地震數
據塊。
12.根據權利要求11所述的方法,其特征在于,
所述根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震數據塊中的每
個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理后的地震數據
道,包括:
分別確定所述GPU中的線程的線程格尺寸和線程塊尺寸;
對所述D個地震數據塊中的每個地震數據塊中的每個地震數據道進行快速傅里葉變
換;
根據快速傅里葉變換結果,確定所述D個地震數據塊中所有地震數據對應的角頻率采
樣間隔△ω;
根據所述角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角頻率值
ωmax,確定所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz;
根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,計算
所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A;
根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,計算所
述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A;
根據所述GPU中的所有線程處理得到的所述每個地震數據道對應的所有時間深度處的
所有頻率點對應的地震幅值,處理得到所述處理后的地震數據道。
13.根據權利要求12所述的方法,其特征在于,
所述GPU中的線程的線程格尺寸為:所述Tracemax為所述GPU
在一個周期內處理的地震數據的道數,b_size=2m,所述m為正整數,且4≤m≤6,所述S為地
震數據道的采樣點數;
所述GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,所述m為正整數,且4
≤m≤6;
所述根據快速傅里葉變換結果,確定所述D個地震數據塊中所有地震數據對應的角頻
率采樣間隔△ω,包括:
根據快速傅里葉變換結果,采用角頻率采樣間隔計算公式確定所述D個地震數據塊中
所有地震數據對應的角頻率采樣間隔△ω;
其中,所述角頻率采樣間隔計算公式為:△ω=2π/(Gdt),所述G為滿足快速傅里葉變
換算法的時間采樣間隔,所述dt為地震數據的時間采樣間隔;
所述根據所述角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角頻率
值ωmax,確定所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,包括:
根據所述目標低角頻率值ωmin和所述角頻率采樣間隔△ω,采用起算角頻率計算公式
計算所述每個地震數據道對應的起算角頻率值ωp;
根據所述目標高角頻率值ωmax和所述角頻率采樣間隔△ω,采用終止角頻率計算公式
計算所述每個地震數據道對應的終止角頻率值ωz;
其中,所述起算角頻率計算公式為:所述終止角頻率計算公式為:
所述△f為頻率增量,所述ωmax+2π△f為通過反Q濾波處理后期望達到的角頻率的上限;
所述根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的
時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,
計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,包括:
根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,采用
地震幅值計算公式計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地
震幅值A;
其中,所述地震幅值計算公式為:所述
A為所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且所述A為復
數,所述k為所述每個頻率點的序號,所述△ω為所述角頻率采樣間隔,所述τj為時間深度,
所述為時間深度τj處的Q值;
所述根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,計
算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A,包括:
根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,采用幅
值和計算公式計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅
值∑A;
其中,所述幅值和計算公式為:
<mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>}</mo> <mo>,</mo> </mrow>
所述A為所述每個地震數據道對應的時間深度τj處的每個
頻率點對應的地震幅值且所述A為復數,所述k為所述每個頻率點的序號,所述△ω為所述
角頻率采樣間隔,所述τj為時間深度,所述為時間深度τj處的Q值,所述Re表示取實部。
14.一種疊后地震數據體處理裝置,其特征在于,用于中央處理器CPU,所述裝置包括:
選取模塊,用于從地震工區的疊后地震數據體對應的地震測線中選取M條樣本測線,所
述M為大于或者等于2的整數,且M=0.003N,所述N為所述疊后地震數據體的地震主測線的
條數,每條所述地震測線對應一個地震剖面;
第一確定模塊,用于從所述地震工區的疊加速度場中確定所述M條樣本測線中的每條
樣本測線對應的速度場,得到所述M條樣本測線對應的樣本速度場ν;
第二確定模塊,用于根據所述M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公式
確定所述M條樣本測線對應的初始Q值場Q0,所述李慶忠經驗公式為Q=α·νβ,α=14,β=
2.2;
校正模塊,用于對所述M條樣本測線對應的初始Q值場Q0進行校正,得到所述M條樣本測
線對應的校正Q值場Q1;
第三確定模塊,用于根據所述M條樣本測線對應的校正Q值場Q1和所述M條樣本測線對應
的加權系數場λfinal,確定所述M條樣本測線對應的最終Q值場Q2;
更新模塊,用于根據所述M條樣本測線對應的樣本速度場ν和所述M條樣本測線對應的
最終Q值場Q2,采用最小二乘法更新所述李慶忠經驗公式,得到更新后的李慶忠經驗公式;
第四確定模塊,用于根據所述地震工區的疊加速度場,采用所述更新后的李慶忠經驗
公式確定所述地震工區的Q值場Qall;
發送模塊,用于向圖形處理器GPU發送所述地震工區的Q值場Qall和所述疊后地震數據
體,以便于所述GPU根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處
理,得到處理后的疊后地震數據體。
15.根據權利要求14所述的裝置,其特征在于,所述裝置還包括:
第一處理模塊,用于根據預設的加權系數組對所述M條樣本測線對應的校正Q值場Q1進
行處理,得到所述M條樣本測線對應的加權Q值場組,所述預設的加權系數組中包括P個加權
系數λ,所述加權Q值場組中包括P套加權Q值場,所述P為大于或者等于1的整數;
第二處理模塊,用于根據所述加權Q值場組對所述M條樣本測線組成的樣本地震數據體
進行反Q濾波處理,得到反Q濾波數據體組,所述反Q濾波數據體組中包括P套反Q濾波數據
體,每套反Q濾波數據體包括M條反Q樣本測線,每條所述反Q樣本測線對應一個地震剖面;
頻譜分析模塊,用于采用傅里葉變換對所述反Q濾波數據體組中的每套反Q濾波數據體
中的每條反Q樣本測線對應的地震剖面進行頻譜分析;
第五確定模塊,用于根據頻譜分析結果進行處理,確定所述M條反Q樣本測線對應的加
權系數場;
第六確定模塊,用于將所述M條反Q樣本測線對應的加權系數場確定為所述M條樣本測
線對應的加權系數場λfinal。
16.根據權利要求15所述的裝置,其特征在于,
所述頻譜分析模塊,具體用于:
在所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地震剖
面上選擇三個矩形時窗,所述三個矩形時窗包括:淺層矩形時窗、中層矩形時窗和深層矩形
時窗;
采用傅里葉變換對所述每條反Q樣本測線對應的地震剖面上位于所述三個矩形時窗中
的每個矩形時窗內的部分進行處理,得到所述每條反Q樣本測線對應的3個反Q樣本頻譜;
根據所述每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻譜,得到所述每套
反Q濾波數據體對應的3M個反Q樣本頻譜;
根據所述反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻譜得到3×M×P個
反Q樣本頻譜,所述3×M×P個反Q樣本頻譜中包括:M×P個淺層矩形時窗對應的反Q樣本頻
譜、M×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層矩形時窗對應的反Q樣本頻譜;
對所述3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。
17.根據權利要求16所述的裝置,其特征在于,
所述第五確定模塊,具體用于:
根據頻譜分析結果確定所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,所述第一目標頻譜為所述每
條反Q樣本測線對應的地震剖面上位于所述淺層矩形時窗內的部分對應的P個反Q樣本頻譜
中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第二目標頻譜為所述每條反Q樣
本測線對應的地震剖面上位于所述中層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫
軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第三目標頻譜為所述每條反Q樣本測線
對應的地震剖面上位于所述深層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成
的封閉區域的面積最大的反Q樣本頻譜;
分別確定所述第一目標頻譜對應的矩形時窗對應的第一目標加權系數、所述第二目標
頻譜對應的矩形時窗對應的第二目標加權系數和所述第三目標頻譜對應的矩形時窗對應
的第三目標加權系數,得到所述每條反Q樣本測線對應的3個目標加權系數,其中,所述第一
目標加權系數、所述第二目標加權系數和所述第三目標加權系數都為所述預設的加權系數
組中的加權系數;
根據所述M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權系數,得到所述M
條反Q樣本測線對應的3M個目標加權系數;
對所述3M個目標加權系數進行空間插值,得到所述地震工區對應的加權系數場;
對所述地震工區對應的加權系數場進行平滑處理,得到處理后的加權系數場;
從所述處理后的加權系數場中抽取出所述M條反Q樣本測線對應的加權系數場。
18.根據權利要求14至17任一所述的裝置,其特征在于,所述裝置還包括:
采樣模塊,用于按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣,得到采樣Q值
場Qs;
分塊模塊,用于對所述疊后地震數據體進行分塊處理,得到D個地震數據塊,所述D為大
于或者等于2的整數;
所述發送模塊,具體用于:
依次向所述GPU發送所述采樣Q值場Qs和所述D個地震數據塊中的每個地震數據塊。
19.根據權利要求18所述的裝置,其特征在于,
所述分塊模塊,具體用于:
計算所述疊后地震數據體中的地震數據的總道數E;
根據所述GPU在一個周期內處理的地震數據的道數Tracemax和所述疊后地震數據體中
的地震數據的總道數E,對所述疊后地震數據體進行分塊處理,得到所述D個地震數據塊;
其中,E=Tracemax×(D-1)+L,所述L為所述D個地震數據塊中的最后一個地震數據塊中
的地震數據的道數。
20.一種疊后地震數據體處理裝置,其特征在于,用于圖形處理器GPU,所述裝置包括:
接收模塊,用于接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據體;
處理模塊,用于根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處
理,得到處理后的疊后地震數據體。
21.根據權利要求20所述的裝置,其特征在于,
所述接收模塊,具體用于:
接收所述CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,所述采樣Q值
場Qs是所述CPU按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣得到的,所述D個地
震數據塊是所述CPU對所述疊后地震數據體進行分塊處理得到的;
所述處理模塊,包括:
處理子模塊,用于根據所述采樣Q值場Qs對所述D個地震數據塊中的每個地震數據塊進
行反Q濾波處理,得到D個處理后的地震數據塊;
得到子模塊,用于根據所述D個處理后的地震數據塊,得到所述處理后的疊后地震數據
體。
22.根據權利要求21所述的裝置,其特征在于,
所述處理子模塊,包括:
第一處理單元,用于對所述采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ;
第二處理單元,用于根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震
數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到所述D個處理后的地震數
據塊。
23.根據權利要求22所述的裝置,其特征在于,
所述第一處理單元,具體用于:
根據所述采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ;
其中,所述Q值插值計算公式為:
<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mi>=</mi> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>所述τj、所述τ(j-1)和所述τ(j+1)都為時間深度,
所述為時間深度τj處的Q值插值,所述為時間深度τ(j-1)處的Q值插值,所述為
時間深度τ(j+1)處的Q值插值,表示對所述取整。
24.根據權利要求22或23所述的裝置,其特征在于,所述第二處理單元,包括:
處理子單元,用于根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震數
據塊中的每個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理后的
地震數據道;
第一得到子單元,用于根據所述每個地震數據塊中的所有處理后的地震數據道,得到
處理后的地震數據塊;
第二得到子單元,用于根據所述D個地震數據塊中的所有處理后的地震數據塊,得到所
述D個處理后的地震數據塊。
25.根據權利要求24所述的裝置,其特征在于,
所述處理子單元,包括:
第一確定子單元,用于分別確定所述GPU中的線程的線程格尺寸和線程塊尺寸;
變換子單元,用于對所述D個地震數據塊中的每個地震數據塊中的每個地震數據道進
行快速傅里葉變換;
第二確定子單元,用于根據快速傅里葉變換結果,確定所述D個地震數據塊中所有地震
數據對應的角頻率采樣間隔△ω;
第三確定子單元,用于根據所述角頻率采樣間隔△ω和預先接收的目標低角頻率值
ωmin、目標高角頻率值ωmax,確定所述每個地震數據道對應的起算角頻率值ωp和終止角頻
率值ωz;
第一計算子單元,用于根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深
度τ處的Q值Qτ中的時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和
終止角頻率值ωz,計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震
幅值A;
第二計算子單元,用于根據所述每個地震數據道對應的時間深度τj處的每個頻率點對
應的地震幅值A,計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震
幅值∑A;
第一處理子單元,用于根據所述GPU中的所有線程處理得到的所述每個地震數據道對
應的所有時間深度處的所有頻率點對應的地震幅值,處理得到所述處理后的地震數據道。
26.根據權利要求25所述的裝置,其特征在于,
所述GPU中的線程的線程格尺寸為:所述Tracemax為所述GPU在
一個周期內處理的地震數據的道數,b_size=2m,所述m為正整數,且4≤m≤6,所述S為地震
數據道的采樣點數;
所述GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,所述m為正整數,且4
≤m≤6;
所述第二確定子單元,具體用于:
根據快速傅里葉變換結果,采用角頻率采樣間隔計算公式確定所述D個地震數據塊中
所有地震數據對應的角頻率采樣間隔△ω;
其中,所述角頻率采樣間隔計算公式為:△ω=2π/(Gdt),所述G為滿足快速傅里葉變
換算法的時間采樣間隔,所述dt為地震數據的時間采樣間隔;
所述第三確定子單元,具體用于:
根據所述目標低角頻率值ωmin和所述角頻率采樣間隔△ω,采用起算角頻率計算公式
計算所述每個地震數據道對應的起算角頻率值ωp;
根據所述目標高角頻率值ωmax和所述角頻率采樣間隔△ω,采用終止角頻率計算公式
計算所述每個地震數據道對應的終止角頻率值ωz;
其中,所述起算角頻率計算公式為:所述終止角頻率計算公式為:
所述△f為頻率增量,所述ωmax+2π△f為通過反Q濾波處理后期望達到
的角頻率的上限;
所述第一計算子單元,具體用于:
根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,采用
地震幅值計算公式計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地
震幅值A;
其中,所述地震幅值計算公式為:所述
A為所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且所述A為復
數,所述k為所述每個頻率點的序號,所述△ω為所述角頻率采樣間隔,所述τj為時間深度,
所述為時間深度τj處的Q值;
所述第二計算子單元,具體用于:
根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,采用幅
值和計算公式計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅
值∑A;
其中,所述幅值和計算公式為:
<mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>}</mo> <mo>,</mo> </mrow>
所述A為所述每個地震數據道對應的時間深度τj處的每個頻
率點對應的地震幅值且所述A為復數,所述k為所述每個頻率點的序號,所述△ω為所述角
頻率采樣間隔,所述τj為時間深度,所述為時間深度τj處的Q值,所述Re表示取實部。

說明書

疊后地震數據體處理方法及裝置

技術領域

本發明涉及地震數據處理領域,特別涉及一種疊后地震數據體處理方法及裝置。

背景技術

地震勘探是一種利用地下介質彈性和密度的差異,通過觀測和分析大地對人工激
發地震波的響應,推斷地下巖層的性質和形態的地球物理勘探方法。其中,可以采用采集設
備對地震波進行采集得到地震數據體,地震數據體可以疊加形成疊后地震數據體。在地震
勘探中,為了便于觀測和分析,通常需要對疊后地震數據體進行處理以提高疊后地震數據
體的分辨率。

現有技術中,可以采用反Q濾波法對疊后地震數據體進行處理,以提高疊后地震數
據體的分辨率。具體地,可以從已知的工區疊加速度場中獲取地震數據體對應的速度場v,
該速度場v中包括地震數據體對應的地震波在地層上各點的傳播速度,然后根據該速度場v
采用李慶忠經驗公式Q=α·vβ計算得到地層上各點的Q值,進而根據地層上各點的Q值對疊
后地震數據體進行反Q濾波處理。

在實現本發明的過程中,發明人發現現有技術至少存在以下問題:

現有技術中采用李慶忠經驗公式確定Q值,由于李慶忠經驗公式中的α、β為固定
值,針對不同的地質狀況,采用李慶忠經驗公式確定的Q值與實際Q值相差較大,因此,現有
技術對疊后地震數據體處理的靈活性較低、可靠性較差。

發明內容

為了解決現有技術中的問題,本發明提供一種疊后地震數據體處理方法及裝置。
所述技術方案如下:

第一方面,提供一種疊后地震數據體處理方法,用于中央處理器CPU,所述方法包
括:

從地震工區的疊后地震數據體對應的地震測線中選取M條樣本測線,所述M為大于
或者等于2的整數,且M=0.003N,所述N為所述疊后地震數據體的地震主測線的條數,每條
所述地震測線對應一個地震剖面;

從所述地震工區的疊加速度場中確定所述M條樣本測線中的每條樣本測線對應的
速度場,得到所述M條樣本測線對應的樣本速度場ν;

根據所述M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公式確定所述M條樣
本測線對應的初始Q值場Q0,所述李慶忠經驗公式為Q=α·νβ,α=14,β=2.2;

對所述M條樣本測線對應的初始Q值場Q0進行校正,得到所述M條樣本測線對應的
校正Q值場Q1;

根據所述M條樣本測線對應的校正Q值場Q1和所述M條樣本測線對應的加權系數場
λfinal,確定所述M條樣本測線對應的最終Q值場Q2;

根據所述M條樣本測線對應的樣本速度場ν和所述M條樣本測線對應的最終Q值場
Q2,采用最小二乘法更新所述李慶忠經驗公式,得到更新后的李慶忠經驗公式;

根據所述地震工區的疊加速度場,采用所述更新后的李慶忠經驗公式確定所述地
震工區的Q值場Qall;

向圖形處理器GPU發送所述地震工區的Q值場Qall和所述疊后地震數據體,以便于
所述GPU根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處理,得到處理
后的疊后地震數據體。

進一步地,在所述根據所述M條樣本測線對應的校正Q值場Q1和所述M條樣本測線
對應的加權系數場λfinal,確定所述M條樣本測線對應的最終Q值場Q2之前,所述方法還包括:

根據預設的加權系數組對所述M條樣本測線對應的校正Q值場Q1進行處理,得到所
述M條樣本測線對應的加權Q值場組,所述預設的加權系數組中包括P個加權系數λ,所述加
權Q值場組中包括P套加權Q值場,所述P為大于或者等于1的整數;

根據所述加權Q值場組對所述M條樣本測線組成的樣本地震數據體進行反Q濾波處
理,得到反Q濾波數據體組,所述反Q濾波數據體組中包括P套反Q濾波數據體,每套反Q濾波
數據體包括M條反Q樣本測線,每條所述反Q樣本測線對應一個地震剖面;

采用傅里葉變換對所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的地震剖面進行頻譜分析;

根據頻譜分析結果進行處理,確定所述M條反Q樣本測線對應的加權系數場;

將所述M條反Q樣本測線對應的加權系數場確定為所述M條樣本測線對應的加權系
數場λfinal。

進一步地,所述采用傅里葉變換對所述反Q濾波數據體組中的每套反Q濾波數據體
中的每條反Q樣本測線對應的地震剖面進行頻譜分析,包括:

在所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地
震剖面上選擇三個矩形時窗,所述三個矩形時窗包括:淺層矩形時窗、中層矩形時窗和深層
矩形時窗;

采用傅里葉變換對所述每條反Q樣本測線對應的地震剖面上位于所述三個矩形時
窗中的每個矩形時窗內的部分進行處理,得到所述每條反Q樣本測線對應的3個反Q樣本頻
譜;

根據所述每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻譜,得到所述
每套反Q濾波數據體對應的3M個反Q樣本頻譜;

根據所述反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻譜得到3×M
×P個反Q樣本頻譜,所述3×M×P個反Q樣本頻譜中包括:M×P個淺層矩形時窗對應的反Q樣
本頻譜、M×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層矩形時窗對應的反Q樣本頻
譜;

對所述3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。

進一步地,所述根據頻譜分析結果進行處理,確定所述M條反Q樣本測線對應的加
權系數場,包括:

根據頻譜分析結果確定所述反Q濾波數據體組中的每套反Q濾波數據體中的每條
反Q樣本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,所述第一目標頻譜為所
述每條反Q樣本測線對應的地震剖面上位于所述淺層矩形時窗內的部分對應的P個反Q樣本
頻譜中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第二目標頻譜為所述每條
反Q樣本測線對應的地震剖面上位于所述中層矩形時窗內的部分對應的P個反Q樣本頻譜中
與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第三目標頻譜為所述每條反Q樣本
測線對應的地震剖面上位于所述深層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸
圍成的封閉區域的面積最大的反Q樣本頻譜;

分別確定所述第一目標頻譜對應的矩形時窗對應的第一目標加權系數、所述第二
目標頻譜對應的矩形時窗對應的第二目標加權系數和所述第三目標頻譜對應的矩形時窗
對應的第三目標加權系數,得到所述每條反Q樣本測線對應的3個目標加權系數,其中,所述
第一目標加權系數、所述第二目標加權系數和所述第三目標加權系數都為所述預設的加權
系數組中的加權系數;

根據所述M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權系數,得到所
述M條反Q樣本測線對應的3M個目標加權系數;

對所述3M個目標加權系數進行空間插值,得到所述地震工區對應的加權系數場;

對所述地震工區對應的加權系數場進行平滑處理,得到處理后的加權系數場;

從所述處理后的加權系數場中抽取出所述M條反Q樣本測線對應的加權系數場。

進一步地,在所述向圖形處理器GPU發送所述地震工區的Q值場Qall和所述疊后地
震數據體之前,所述方法還包括:

按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣,得到采樣Q值場Qs;

對所述疊后地震數據體進行分塊處理,得到D個地震數據塊,所述D為大于或者等
于2的整數;

所述向圖形處理器GPU發送所述地震工區的Q值場Qall和疊后地震數據體,包括:

依次向所述GPU發送所述采樣Q值場Qs和所述D個地震數據塊中的每個地震數據
塊。

進一步地,所述對所述疊后地震數據體進行分塊處理,得到D個地震數據塊,包括:

計算所述疊后地震數據體中的地震數據的總道數E;

根據所述GPU在一個周期內處理的地震數據的道數Tracemax和所述疊后地震數據
體中的地震數據的總道數E,對所述疊后地震數據體進行分塊處理,得到所述D個地震數據
塊;

其中,E=Tracemax×(D-1)+L,所述L為所述D個地震數據塊中的最后一個地震數據
塊中的地震數據的道數。

第二方面,提供一種疊后地震數據體處理方法,用于圖形處理器GPU,所述方法包
括:

接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據體;

根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處理,得到處
理后的疊后地震數據體。

進一步地,所述接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據
體,包括:

接收所述CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,所述采
樣Q值場Qs是所述CPU按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣得到的,所述
D個地震數據塊是所述CPU對所述疊后地震數據體進行分塊處理得到的;

所述根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波處理,得
到處理后的疊后地震數據體,包括:

根據所述采樣Q值場Qs對所述D個地震數據塊中的每個地震數據塊進行反Q濾波處
理,得到D個處理后的地震數據塊;

根據所述D個處理后的地震數據塊,得到所述處理后的疊后地震數據體。

進一步地,所述根據所述采樣Q值場Qs對所述D個地震數據塊中的每個地震數據塊
進行反Q濾波處理,得到D個處理后的地震數據塊,包括:

對所述采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ;

根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震數據塊中的所
有時間深度τ處的地震幅值進行反Q濾波處理,得到所述D個處理后的地震數據塊。

進一步地,所述對所述采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ,包括:

根據所述采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ;

其中,所述Q值插值計算公式為:

<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>

<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>所述τj、所述τ(j-1)和所述
τ(j+1)都為時間深度,所述為時間深度τj處的Q值插值,所述為時間深度τ(j-1)處的Q值
插值,所述為時間深度τ(j+1)處的Q值插值,表示對所述取整。

進一步地,所述根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地
震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到所述D個處理后的地震
數據塊,包括:

根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地震數據塊中的每
個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理后的地震數據
道;

根據所述每個地震數據塊中的所有處理后的地震數據道,得到處理后的地震數據
塊;

根據所述D個地震數據塊中的所有處理后的地震數據塊,得到所述D個處理后的地
震數據塊。

進一步地,所述根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個地
震數據塊中的每個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理
后的地震數據道,包括:

分別確定所述GPU中的線程的線程格尺寸和線程塊尺寸;

對所述D個地震數據塊中的每個地震數據塊中的每個地震數據道進行快速傅里葉
變換;

根據快速傅里葉變換結果,確定所述D個地震數據塊中所有地震數據對應的角頻
率采樣間隔△ω;

根據所述角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角頻率
值ωmax,確定所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz;

根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的
時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,
計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A;

根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,
計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A;

根據所述GPU中的所有線程處理得到的所述每個地震數據道對應的所有時間深度
處的所有頻率點對應的地震幅值,處理得到所述處理后的地震數據道。

進一步地,所述GPU中的線程的線程格尺寸為:所述
Tracemax為所述GPU在一個周期內處理的地震數據的道數,b_size=2m,所述m為正整數,且4
≤m≤6,所述S為地震數據道的采樣點數;

所述GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,所述m為正整
數,且4≤m≤6;

所述根據快速傅里葉變換結果,確定所述D個地震數據塊中所有地震數據對應的
角頻率采樣間隔△ω,包括:

根據快速傅里葉變換結果,采用角頻率采樣間隔計算公式確定所述D個地震數據
塊中所有地震數據對應的角頻率采樣間隔△ω;

其中,所述角頻率采樣間隔計算公式為:△ω=2π/(Gdt),所述G為滿足快速傅里
葉變換算法的時間采樣間隔,所述dt為地震數據的時間采樣間隔;

所述根據所述角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角
頻率值ωmax,確定所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,包括:

根據所述目標低角頻率值ωmin和所述角頻率采樣間隔△ω,采用起算角頻率計算
公式計算所述每個地震數據道對應的起算角頻率值ωp;

根據所述目標高角頻率值ωmax和所述角頻率采樣間隔△ω,采用終止角頻率計算
公式計算所述每個地震數據道對應的終止角頻率值ωz;

其中,所述起算角頻率計算公式為:所述終止角頻率計算公式為:
所述△f為頻率增量,所述ωmax+2π△f為通過反Q濾波處理后期望達到
的角頻率的上限;

所述根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ
中的時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值
ωz,計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,包括:

根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的
時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,
采用地震幅值計算公式計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應
的地震幅值A;

其中,所述地震幅值計算公式為:
所述A為所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且所述A
為復數,所述k為所述每個頻率點的序號,所述△ω為所述角頻率采樣間隔,所述τj為時間
深度,所述為時間深度τj處的Q值;

所述根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅
值A,計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A,包
括:

根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,
采用幅值和計算公式計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的
地震幅值∑A;

其中,所述幅值和計算公式為:

<mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> <mo>}</mo> <mo>,</mo> </mrow>
述A為所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且所述A為
復數,所述k為所述每個頻率點的序號,所述△ω為所述角頻率采樣間隔,所述τj為時間深
度,所述為時間深度τj處的Q值,所述Re表示取實部。

第三方面,提供一種疊后地震數據體處理裝置,用于中央處理器CPU,所述裝置包
括:

選取模塊,用于從地震工區的疊后地震數據體對應的地震測線中選取M條樣本測
線,所述M為大于或者等于2的整數,且M=0.003N,所述N為所述疊后地震數據體的地震主測
線的條數,每條所述地震測線對應一個地震剖面;

第一確定模塊,用于從所述地震工區的疊加速度場中確定所述M條樣本測線中的
每條樣本測線對應的速度場,得到所述M條樣本測線對應的樣本速度場ν;

第二確定模塊,用于根據所述M條樣本測線對應的樣本速度場ν,采用李慶忠經驗
公式確定所述M條樣本測線對應的初始Q值場Q0,所述李慶忠經驗公式為Q=α·νβ,α=14,β
=2.2;

校正模塊,用于對所述M條樣本測線對應的初始Q值場Q0進行校正,得到所述M條樣
本測線對應的校正Q值場Q1;

第三確定模塊,用于根據所述M條樣本測線對應的校正Q值場Q1和所述M條樣本測
線對應的加權系數場λfinal,確定所述M條樣本測線對應的最終Q值場Q2;

更新模塊,用于根據所述M條樣本測線對應的樣本速度場ν和所述M條樣本測線對
應的最終Q值場Q2,采用最小二乘法更新所述李慶忠經驗公式,得到更新后的李慶忠經驗公
式;

第四確定模塊,用于根據所述地震工區的疊加速度場,采用所述更新后的李慶忠
經驗公式確定所述地震工區的Q值場Qall;

發送模塊,用于向圖形處理器GPU發送所述地震工區的Q值場Qall和所述疊后地震
數據體,以便于所述GPU根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾波
處理,得到處理后的疊后地震數據體。

進一步地,所述裝置還包括:

第一處理模塊,用于根據預設的加權系數組對所述M條樣本測線對應的校正Q值場
Q1進行處理,得到所述M條樣本測線對應的加權Q值場組,所述預設的加權系數組中包括P個
加權系數λ,所述加權Q值場組中包括P套加權Q值場,所述P為大于或者等于1的整數;

第二處理模塊,用于根據所述加權Q值場組對所述M條樣本測線組成的樣本地震數
據體進行反Q濾波處理,得到反Q濾波數據體組,所述反Q濾波數據體組中包括P套反Q濾波數
據體,每套反Q濾波數據體包括M條反Q樣本測線,每條所述反Q樣本測線對應一個地震剖面;

頻譜分析模塊,用于采用傅里葉變換對所述反Q濾波數據體組中的每套反Q濾波數
據體中的每條反Q樣本測線對應的地震剖面進行頻譜分析;

第五確定模塊,用于根據頻譜分析結果進行處理,確定所述M條反Q樣本測線對應
的加權系數場;

第六確定模塊,用于將所述M條反Q樣本測線對應的加權系數場確定為所述M條樣
本測線對應的加權系數場λfinal。

進一步地,所述頻譜分析模塊,具體用于:

在所述反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地
震剖面上選擇三個矩形時窗,所述三個矩形時窗包括:淺層矩形時窗、中層矩形時窗和深層
矩形時窗;

采用傅里葉變換對所述每條反Q樣本測線對應的地震剖面上位于所述三個矩形時
窗中的每個矩形時窗內的部分進行處理,得到所述每條反Q樣本測線對應的3個反Q樣本頻
譜;

根據所述每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻譜,得到所述
每套反Q濾波數據體對應的3M個反Q樣本頻譜;

根據所述反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻譜得到3×M
×P個反Q樣本頻譜,所述3×M×P個反Q樣本頻譜中包括:M×P個淺層矩形時窗對應的反Q樣
本頻譜、M×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層矩形時窗對應的反Q樣本頻
譜;

對所述3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。

進一步地,所述第五確定模塊,具體用于:

根據頻譜分析結果確定所述反Q濾波數據體組中的每套反Q濾波數據體中的每條
反Q樣本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,所述第一目標頻譜為所
述每條反Q樣本測線對應的地震剖面上位于所述淺層矩形時窗內的部分對應的P個反Q樣本
頻譜中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第二目標頻譜為所述每條
反Q樣本測線對應的地震剖面上位于所述中層矩形時窗內的部分對應的P個反Q樣本頻譜中
與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,所述第三目標頻譜為所述每條反Q樣本
測線對應的地震剖面上位于所述深層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸
圍成的封閉區域的面積最大的反Q樣本頻譜;

分別確定所述第一目標頻譜對應的矩形時窗對應的第一目標加權系數、所述第二
目標頻譜對應的矩形時窗對應的第二目標加權系數和所述第三目標頻譜對應的矩形時窗
對應的第三目標加權系數,得到所述每條反Q樣本測線對應的3個目標加權系數,其中,所述
第一目標加權系數、所述第二目標加權系數和所述第三目標加權系數都為所述預設的加權
系數組中的加權系數;

根據所述M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權系數,得到所
述M條反Q樣本測線對應的3M個目標加權系數;

對所述3M個目標加權系數進行空間插值,得到所述地震工區對應的加權系數場;

對所述地震工區對應的加權系數場進行平滑處理,得到處理后的加權系數場;

從所述處理后的加權系數場中抽取出所述M條反Q樣本測線對應的加權系數場。

進一步地,所述裝置還包括:

采樣模塊,用于按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣,得到采
樣Q值場Qs;

分塊模塊,用于對所述疊后地震數據體進行分塊處理,得到D個地震數據塊,所述D
為大于或者等于2的整數;

所述發送模塊,具體用于:

依次向所述GPU發送所述采樣Q值場Qs和所述D個地震數據塊中的每個地震數據
塊。

進一步地,所述分塊模塊,具體用于:

計算所述疊后地震數據體中的地震數據的總道數E;

根據所述GPU在一個周期內處理的地震數據的道數Tracemax和所述疊后地震數據
體中的地震數據的總道數E,對所述疊后地震數據體進行分塊處理,得到所述D個地震數據
塊;

其中,E=Tracemax×(D-1)+L,所述L為所述D個地震數據塊中的最后一個地震數據
塊中的地震數據的道數。

第四方面,提供一種疊后地震數據體處理裝置,用于圖形處理器GPU,所述裝置包
括:

接收模塊,用于接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據
體;

處理模塊,用于根據所述地震工區的Q值場Qall對所述疊后地震數據體進行反Q濾
波處理,得到處理后的疊后地震數據體。

進一步地,所述接收模塊,具體用于:

接收所述CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,所述采
樣Q值場Qs是所述CPU按照預設采樣間隔對所述地震工區的Q值場Qall進行采樣得到的,所述
D個地震數據塊是所述CPU對所述疊后地震數據體進行分塊處理得到的;

所述處理模塊,包括:

處理子模塊,用于根據所述采樣Q值場Qs對所述D個地震數據塊中的每個地震數據
塊進行反Q濾波處理,得到D個處理后的地震數據塊;

得到子模塊,用于根據所述D個處理后的地震數據塊,得到所述處理后的疊后地震
數據體。

進一步地,所述處理子模塊,包括:

第一處理單元,用于對所述采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ;

第二處理單元,用于根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每
個地震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到所述D個處理后的
地震數據塊。

進一步地,所述第一處理單元,具體用于:

根據所述采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ;

其中,所述Q值插值計算公式為:

<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>

<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>所述τj、所述τ(j-1)和所述τ(j+1)都
為時間深度,所述為時間深度τj處的Q值插值,所述為時間深度τ(j-1)處的Q值插值,
所述為時間深度τ(j+1)處的Q值插值,表示對所述取整。

進一步地,所述第二處理單元,包括:

處理子單元,用于根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個
地震數據塊中的每個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處
理后的地震數據道;

第一得到子單元,用于根據所述每個地震數據塊中的所有處理后的地震數據道,
得到處理后的地震數據塊;

第二得到子單元,用于根據所述D個地震數據塊中的所有處理后的地震數據塊,得
到所述D個處理后的地震數據塊。

進一步地,所述處理子單元,包括:

第一確定子單元,用于分別確定所述GPU中的線程的線程格尺寸和線程塊尺寸;

變換子單元,用于對所述D個地震數據塊中的每個地震數據塊中的每個地震數據
道進行快速傅里葉變換;

第二確定子單元,用于根據快速傅里葉變換結果,確定所述D個地震數據塊中所有
地震數據對應的角頻率采樣間隔△ω;

第三確定子單元,用于根據所述角頻率采樣間隔△ω和預先接收的目標低角頻率
值ωmin、目標高角頻率值ωmax,確定所述每個地震數據道對應的起算角頻率值ωp和終止角
頻率值ωz;

第一計算子單元,用于根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時
間深度τ處的Q值Qτ中的時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值
ωp和終止角頻率值ωz,計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應
的地震幅值A;

第二計算子單元,用于根據所述每個地震數據道對應的時間深度τj處的每個頻率
點對應的地震幅值A,計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的
地震幅值∑A;

第一處理子單元,用于根據所述GPU中的所有線程處理得到的所述每個地震數據
道對應的所有時間深度處的所有頻率點對應的地震幅值,處理得到所述處理后的地震數據
道。

進一步地,所述GPU中的線程的線程格尺寸為:所述
Tracemax為所述GPU在一個周期內處理的地震數據的道數,b_size=2m,所述m為正整數,且4
≤m≤6,所述S為地震數據道的采樣點數;

所述GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,所述m為正整
數,且4≤m≤6;

所述第二確定子單元,具體用于:

根據快速傅里葉變換結果,采用角頻率采樣間隔計算公式確定所述D個地震數據
塊中所有地震數據對應的角頻率采樣間隔△ω;

其中,所述角頻率采樣間隔計算公式為:△ω=2π/(Gdt),所述G為滿足快速傅里
葉變換算法的時間采樣間隔,所述dt為地震數據的時間采樣間隔;

所述第三確定子單元,具體用于:

根據所述目標低角頻率值ωmin和所述角頻率采樣間隔△ω,采用起算角頻率計算
公式計算所述每個地震數據道對應的起算角頻率值ωp;

根據所述目標高角頻率值ωmax和所述角頻率采樣間隔△ω,采用終止角頻率計算
公式計算所述每個地震數據道對應的終止角頻率值ωz;

其中,所述起算角頻率計算公式為:所述終止角頻率計算公式為:
所述△f為頻率增量,所述ωmax+2π△f為通過反Q濾波處理后期望達到
的角頻率的上限;

所述第一計算子單元,具體用于:

根據所述GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的
時間深度τj處的Q值所述每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,
采用地震幅值計算公式計算所述每個地震數據道對應的時間深度τj處的每個頻率點對應
的地震幅值A;

其中,所述地震幅值計算公式為:
所述A為所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且所述A
為復數,所述k為所述每個頻率點的序號,所述△ω為所述角頻率采樣間隔,所述τj為時間
深度,所述為時間深度τj處的Q值;

所述第二計算子單元,具體用于:

根據所述每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,
采用幅值和計算公式計算所述每個地震數據道對應的時間深度τj處的所有頻率點對應的
地震幅值∑A;

其中,所述幅值和計算公式為:

<mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> <mo>}</mo> <mo>,</mo> </mrow>所述A為所述
每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且所述A為復數,所述k
為所述每個頻率點的序號,所述△ω為所述角頻率采樣間隔,所述τj為時間深度,所述為
時間深度τj處的Q值,所述Re表示取實部。

本發明提供的技術方案帶來的有益效果是:

本發明提供的疊后地震數據體處理方法及裝置,方法包括:從地震工區的疊后地
震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣本測
線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣本測
線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條樣本
測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測線對
應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值場;根
據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法更新
李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更新后
的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數據體,
以便于GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理后的疊
后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決了對疊
后地震數據體處理的靈活性較低、可靠性較差的問題。

進一步地,通過CPU對疊后地震數據體進行分塊處理得到地震數據塊,使得GPU根
據地震數據塊對疊后地震數據體進行處理,由于GPU的處理能力較強,本發明通過采用GPU
對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效率。

應當理解的是,以上的一般描述和后文的細節描述僅是示例性的,并不能限制本
發明。

附圖說明

為了更清楚地說明本發明實施例中的技術方案,下面將對實施例描述中所需要使
用的附圖作簡單地介紹,顯而易見地,下面描述中的附圖僅僅是本發明的一些實施例,對于
本領域普通技術人員來講,在不付出創造性勞動的前提下,還可以根據這些附圖獲得其他
的附圖。

圖1是本發明實施例提供的一種疊后地震數據體處理方法的方法流程圖;

圖2是本發明實施例提供的另一種疊后地震數據體處理方法的方法流程圖;

圖3-1是本發明實施例提供的再一種疊后地震數據體處理方法的方法流程圖;

圖3-2是圖3-1所示實施例提供的一種采用傅里葉變換對每條反Q樣本測線對應的
地震剖面進行頻譜分析的方法流程圖;

圖3-3是圖3-1所示實施例提供的一種根據頻譜分析結果進行處理確定M條反Q樣
本測線對應的加權系數場的方法流程圖;

圖3-4是圖3-1所示實施例提供的一種對疊后地震數據體進行分塊處理的方法流
程圖;

圖3-5是圖3-1所示實施例提供的一種根據地震工區的Q值場對疊后地震數據體進
行反Q濾波處理的方法流程圖;

圖3-6是圖3-1所示實施例提供的一種根據采樣Q值場對D個地震數據塊中的每個
地震數據塊進行反Q濾波處理的方法流程圖;

圖3-7是圖3-1所示實施例提供的一種根據所有時間深度處的Q值對每個地震數據
塊中的所有時間深度處的地震幅值進行反Q濾波處理的方法流程圖;

圖3-8是圖3-1所示實施例提供的一種根據所有時間深度處的Q值對每個地震數據
道的所有時間深度處的地震幅值進行反Q濾波處理的方法流程圖;

圖3-9是本發明實施例提供的一種處理前的地震剖面圖;

圖3-10是圖3-9所示的地震剖面圖上的部分區域的放大圖;

圖3-11是本發明實施例提供的對圖3-9所示的地震剖面圖進行處理得到的地震剖
面圖;

圖3-12是圖3-11所示的地震剖面圖上的部分區域的放大圖;

圖3-13是本發明實施例提供的一種處理前的水平時間切片圖;

圖3-14是是本發明實施例提供的對圖3-13所示的水平時間切片圖進行處理得到
的水平時間切片圖;

圖4是本發明實施例提供的一種疊后地震數據體處理裝置的框圖;

圖5是本發明實施例提供的另一種疊后地震數據體處理裝置的框圖;

圖6是本發明實施例提供的再一種疊后地震數據體處理裝置的框圖;

圖7-1是本發明實施例提供的又一種疊后地震數據體處理裝置的框圖;

圖7-2是圖7-1所示實施例提供的一種處理模塊的框圖;

圖7-3是圖7-1所示實施例提供的一種處理子模塊的框圖;

圖7-4是圖7-1所示實施例提供的一種第二處理單元的框圖;

圖7-5是圖7-1所示實施例提供的一種處理子單元的框圖。

此處的附圖被并入說明書中并構成本說明書的一部分,示出了符合本發明的實施
例,并與說明書一起用于解釋本發明的原理。

具體實施方式

為了使本發明的目的、技術方案和優點更加清楚,下面將結合附圖對本發明作進
一步地詳細描述,顯然,所描述的實施例僅僅是本發明一部份實施例,而不是全部的實施
例。基于本發明中的實施例,本領域普通技術人員在沒有做出創造性勞動前提下所獲得的
所有其它實施例,都屬于本發明保護的范圍。

請參考圖1,其示出了本發明實施例提供的一種疊后地震數據體處理方法的方法
流程圖,該疊后地震數據體處理方法可以用于中央處理器(英文:Central Processing
Unit;簡稱:CPU),參見圖1,該方法流程具體包括:

在步驟101中,從地震工區的疊后地震數據體對應的地震測線中選取M條樣本測
線。

其中,M為大于或者等于2的整數,且M=0.003N,N為疊后地震數據體的地震主測線
的條數,每條地震測線對應一個地震剖面。

在步驟102中,從地震工區的疊加速度場中確定M條樣本測線中的每條樣本測線對
應的速度場,得到M條樣本測線對應的樣本速度場ν。

在步驟103中,根據M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公式確定M
條樣本測線對應的初始Q值場Q0。

其中,李慶忠經驗公式為Q=α·νβ,α=14,β=2.2。

在步驟104中,對M條樣本測線對應的初始Q值場Q0進行校正,得到M條樣本測線對
應的校正Q值場Q1。

在步驟105中,根據M條樣本測線對應的校正Q值場Q1和M條樣本測線對應的加權系
數場λfinal,確定M條樣本測線對應的最終Q值場Q2。

在步驟106中,根據M條樣本測線對應的樣本速度場ν和M條樣本測線對應的最終Q
值場Q2,采用最小二乘法更新李慶忠經驗公式,得到更新后的李慶忠經驗公式。

在步驟107中,根據地震工區的疊加速度場,采用更新后的李慶忠經驗公式確定地
震工區的Q值場Qall。

在步驟108中,向圖形處理器GPU發送地震工區的Q值場Qall和疊后地震數據體,以
便于GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理后的疊后
地震數據體。

綜上所述,本發明實施例提供的疊后地震數據體處理方法,通過從地震工區的疊
后地震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣
本測線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣
本測線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條
樣本測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測
線對應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值
場;根據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法
更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更
新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數
據體,以便于GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理
后的疊后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決
了對疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體
處理的靈活性和可靠性的效果。

進一步地,在步驟105之前,該方法還包括:

根據預設的加權系數組對M條樣本測線對應的校正Q值場Q1進行處理,得到M條樣
本測線對應的加權Q值場組,預設的加權系數組中包括P個加權系數λ,加權Q值場組中包括P
套加權Q值場,P為大于或者等于1的整數;

根據加權Q值場組對M條樣本測線組成的樣本地震數據體進行反Q濾波處理,得到
反Q濾波數據體組,反Q濾波數據體組中包括P套反Q濾波數據體,每套反Q濾波數據體包括M
條反Q樣本測線,每條反Q樣本測線對應一個地震剖面;

采用傅里葉變換對反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測
線對應的地震剖面進行頻譜分析;

根據頻譜分析結果進行處理,確定M條反Q樣本測線對應的加權系數場;

將M條反Q樣本測線對應的加權系數場確定為M條樣本測線對應的加權系數場
λfinal。

進一步地,采用傅里葉變換對反Q濾波數據體組中的每套反Q濾波數據體中的每條
反Q樣本測線對應的地震剖面進行頻譜分析,包括:

在反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地震剖
面上選擇三個矩形時窗,三個矩形時窗包括:淺層矩形時窗、中層矩形時窗和深層矩形時
窗;

采用傅里葉變換對每條反Q樣本測線對應的地震剖面上位于三個矩形時窗中的每
個矩形時窗內的部分進行處理,得到每條反Q樣本測線對應的3個反Q樣本頻譜;

根據每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻譜,得到每套反Q
濾波數據體對應的3M個反Q樣本頻譜;

根據反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻譜得到3×M×P個
反Q樣本頻譜,3×M×P個反Q樣本頻譜中包括:M×P個淺層矩形時窗對應的反Q樣本頻譜、M
×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層矩形時窗對應的反Q樣本頻譜;

對3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。

進一步地,根據頻譜分析結果進行處理,確定M條反Q樣本測線對應的加權系數場,
包括:

根據頻譜分析結果確定反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,第一目標頻譜為每條反Q樣本
測線對應的地震剖面上位于淺層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成
的封閉區域的面積最大的反Q樣本頻譜,第二目標頻譜為每條反Q樣本測線對應的地震剖面
上位于中層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成的封閉區域的面積最
大的反Q樣本頻譜,第三目標頻譜為每條反Q樣本測線對應的地震剖面上位于深層矩形時窗
內的部分對應的P個反Q樣本頻譜中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜;

分別確定第一目標頻譜對應的矩形時窗對應的第一目標加權系數、第二目標頻譜
對應的矩形時窗對應的第二目標加權系數和第三目標頻譜對應的矩形時窗對應的第三目
標加權系數,得到每條反Q樣本測線對應的3個目標加權系數,其中,第一目標加權系數、第
二目標加權系數和第三目標加權系數都為預設的加權系數組中的加權系數;

根據M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權系數,得到M條反Q
樣本測線對應的3M個目標加權系數;

對3M個目標加權系數進行空間插值,得到地震工區對應的加權系數場;

對地震工區對應的加權系數場進行平滑處理,得到處理后的加權系數場;

從處理后的加權系數場中抽取出M條反Q樣本測線對應的加權系數場。

進一步地,在步驟108之前,該方法還包括:

按照預設采樣間隔對地震工區的Q值場Qall進行采樣,得到采樣Q值場Qs;

對疊后地震數據體進行分塊處理,得到D個地震數據塊,D為大于或者等于2的整
數;

向圖形處理器GPU發送地震工區的Q值場Qall和疊后地震數據體,包括:

依次向GPU發送采樣Q值場Qs和D個地震數據塊中的每個地震數據塊。

進一步地,對疊后地震數據體進行分塊處理,得到D個地震數據塊,包括:

計算疊后地震數據體中的地震數據的總道數E;

根據GPU在一個周期內處理的地震數據的道數Tracemax和疊后地震數據體中的地
震數據的總道數E,對疊后地震數據體進行分塊處理,得到D個地震數據塊;

其中,E=Tracemax×(D-1)+L,L為D個地震數據塊中的最后一個地震數據塊中的地
震數據的道數。

進一步地,步驟104可以包括:

采用譜比法確定M條樣本測線對應的層Q值場Qi,i為垂直地震剖面VSP測井中的地
層序號;

根據M條樣本測線對應的層Q值場Qi,確定M條樣本測線對應的時間深度的疊加Q值
場Qvsp;

根據M條樣本測線對應的時間深度的疊加Q值場Qvsp,對M條樣本測線對應的初始Q
值場Q0進行校正,得到M條樣本測線對應的校正Q值場Q1。

進一步地,步驟106可以包括:

根據M條樣本測線對應的樣本速度場ν和M條樣本測線對應的最終Q值場Q2,采用最
小二乘法確定李慶忠經驗公式Q=α·νβ中的α值與β值,得到更新α值和更新β值;

將更新α值和更新β值代入李慶忠經驗公式,得到更新后的李慶忠經驗公式。

進一步地,根據M條樣本測線對應的層Q值場Qi,確定M條樣本測線對應的時間深度
的疊加Q值場Qvsp,包括:

根據M條樣本測線對應的層Q值場Qi,利用疊加Q值場公式確定M條樣本測線對應的
時間深度的疊加Q值場Qvsp;

疊加Q值場公式為:

其中,i為垂直地震剖面VSP測井中的地層序號,△Ti為M條樣本測線對
應的地震波在第i層地層中的單程旅行時間,n為地層總數,T為M條樣本測線對應的地震波
從地表垂直傳播至第i層地層的時間。

進一步地,在反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對
應的地震剖面上選擇三個矩形時窗,包括:

在反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地震剖
面上選擇三個矩形時窗,使每個矩形時窗的大小相等,且使每個矩形時窗內的同相軸的清
晰度大于預設清晰度,同相軸的連續度大于預設連續度。

進一步地,在步驟108之前,該方法還包括:

對疊后地震數據體進行去噪處理。

進一步地,步驟101可以包括:

從地震工區的疊后地震數據體對應的地震測線中選取M條存在垂直地震剖面VSP
測井的地震測線作為M條樣本測線;

其中,M條樣本測線包括:地震主測線和聯絡測線中的至少一種。

進一步地,從地震工區的疊后地震數據體對應的地震測線中選取M條存在垂直地
震剖面VSP測井的地震測線作為M條樣本測線,包括:

根據大地坐標,將地震工區的所有VSP測井投射到工區平面圖上;

根據工區平面圖,從地震工區的疊后地震數據體對應的地震測線中選取M條存在
垂直地震剖面VSP測井的地震測線作為M條樣本測線。

進一步地,采用譜比法確定M條樣本測線對應的層Q值場Qi,i為垂直地震剖面VSP
測井中的地層序號,包括:

根據層Q值場公式確定M條樣本測線對應的層Q值場Qi;

層Q值場為: <mrow> <mi>l</mi> <mi>n</mi> <mo>|</mo> <mfrac> <mrow> <msub> <mi>a</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>a</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>|</mo> <mo>=</mo> <mi>c</mi> <mi>o</mi> <mi>n</mi> <mi>s</mi> <mi>t</mi> <mo>+</mo> <mi>&pi;</mi> <mfrac> <mrow> <mi>&delta;</mi> <mi>t</mi> </mrow> <msub> <mi>Q</mi> <mi>i</mi> </msub> </mfrac> <mi>f</mi> <mo>;</mo> </mrow>

其中,f為M條樣本測線中地層的頻譜中的頻率值,ai(f)為M條樣本測線中第i層地
層的頻譜中頻率f對應的地震幅值,ai-1(f)為M條樣本測線中第i-1層地層的頻譜中頻率f對
應的地震幅值,const表示常數,δt為地震波到達第i層地層的時間與地震波到達第i-1層地
層的時間的差值。

進一步地,在步驟108之前,該方法還可以包括:

采用傅里葉變換對M條樣本測線中的每條樣本測線對應的地震剖面進行頻譜分
析,得到M條樣本測線對應的目標低角頻率值ωmin和目標高角頻率值ωmax;

該方法還可以包括:

向GPU發送目標低角頻率值ωmin和目標高角頻率值ωmax。

進一步地,采用傅里葉變換對M條樣本測線中的每條樣本測線對應的地震剖面進
行頻譜分析,得到M條樣本測線對應的目標低角頻率值ωmin和目標高角頻率值ωmax,包括:

在M條樣本測線中的每條樣本測線對應的地震剖面上選擇三個矩形時窗,三個矩
形時窗包括:淺層矩形時窗、中層矩形時窗和深層矩形時窗;

采用傅里葉變換對每條樣本測線對應的地震剖面上位于三個矩形時窗中的每個
矩形時窗內的部分進行處理,得到每條樣本測線對應的3個樣本頻譜;

根據M條樣本測線中的每條樣本測線對應的3個樣本頻譜,得到M條樣本測線對應
的3M個樣本頻譜;

確定3M個樣本頻譜中的每個樣本頻譜對應的最高角頻率值和最低角頻率值,得到
3M個最高角頻率值和3M個最低角頻率值;

將3M個最低角頻率值中最小的角頻率值作為M條樣本測線對應的目標低角頻率值
ωmin;

將3M個最高角頻率值中最大的角頻率值作為M條樣本測線對應的目標高角頻率值
ωmax。

綜上所述,本發明實施例提供的疊后地震數據體處理方法,通過從地震工區的疊
后地震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣
本測線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣
本測線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條
樣本測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測
線對應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值
場;根據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法
更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更
新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數
據體,以便于GPU根據地震工區的Q值場對疊后地震數據體進行反Q濾波處理,得到處理后的
疊后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決了對
疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體處理
的靈活性和可靠性的效果。

進一步地,通過CPU對疊后地震數據體進行分塊處理得到地震數據塊,使得GPU根
據地震數據塊對疊后地震數據體進行處理,由于GPU的處理能力較強,本發明實施例通過采
用GPU對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效
率。

請參考圖2,其示出了本發明實施例提供的另一種疊后地震數據體處理方法的方
法流程圖,該疊后地震數據體處理方法可以用于圖形處理器(英文:Graphic Processing
Unit;簡稱:GPU),參見圖2,該方法流程具體包括:

在步驟201中,接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數據
體。

在步驟202中,根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得
到處理后的疊后地震數據體。

綜上所述,本發明實施例提供的疊后地震數據體處理方法,通過接收CPU發送的地
震工區的Q值場和疊后地震數據體,根據地震工區的Q值場對疊后地震數據體進行反Q濾波
處理,得到處理后的疊后地震數據體,由于地震工區的Q值場是CPU根據更新后的李慶忠經
驗公式計算得到的,且更新后的李慶忠經驗公式是CPU針對疊后地震數據體對李慶忠經驗
公式進行更新得到的,因此,解決了對疊后地震數據體處理的靈活性較低、可靠性較差的問
題,達到了提高疊后地震數據體處理的靈活性和可靠性的效果。

進一步地,步驟201可以包括:

接收CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,采樣Q值場Qs
是CPU按照預設采樣間隔對地震工區的Q值場Qall進行采樣得到的,D個地震數據塊是CPU對
疊后地震數據體進行分塊處理得到的;

步驟202可以包括:

根據采樣Q值場Qs對D個地震數據塊中的每個地震數據塊進行反Q濾波處理,得到D
個處理后的地震數據塊;

根據D個處理后的地震數據塊,得到處理后的疊后地震數據體。

進一步地,根據采樣Q值場Qs對D個地震數據塊中的每個地震數據塊進行反Q濾波
處理,得到D個處理后的地震數據塊,包括:

對采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ;

根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個地震數據塊中的所有時
間深度τ處的地震幅值進行反Q濾波處理,得到D個處理后的地震數據塊。

進一步地,對采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ,包括:

根據采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ;

其中,Q值插值計算公式為:

<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>

<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>τj、τ(j-1)和τ(j+1)都為時間深度,為
時間深度τj處的Q值插值,為時間深度τ(j-1)處的Q值插值,為時間深度τ(j+1)處的Q
值插值,表示對取整。

進一步地,根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個地震數據塊
中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到D個處理后的地震數據塊,包括:

根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個地震數據塊中的每個地
震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理后的地震數據道;

根據每個地震數據塊中的所有處理后的地震數據道,得到處理后的地震數據塊;

根據D個地震數據塊中的所有處理后的地震數據塊,得到D個處理后的地震數據
塊。

進一步地,根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個地震數據塊
中的每個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理后的地震
數據道,包括:

分別確定GPU中的線程的線程格尺寸和線程塊尺寸;

對D個地震數據塊中的每個地震數據塊中的每個地震數據道進行快速傅里葉變
換;

根據快速傅里葉變換結果,確定D個地震數據塊中所有地震數據對應的角頻率采
樣間隔△ω;

根據角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角頻率值
ωmax,確定每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz;

根據GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,計算每個
地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A;

根據每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,計算
每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A;

根據GPU中的所有線程處理得到的每個地震數據道對應的所有時間深度處的所有
頻率點對應的地震幅值,處理得到處理后的地震數據道。

進一步地,GPU中的線程的線程格尺寸為:Tracemax為GPU在
一個周期內處理的地震數據的道數,b_size=2m,m為正整數,且4≤m≤6,S為地震數據道的
采樣點數;

GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,m為正整數,且4≤m
≤6;

根據快速傅里葉變換結果,確定D個地震數據塊中所有地震數據對應的角頻率采
樣間隔△ω,包括:

根據快速傅里葉變換結果,采用角頻率采樣間隔計算公式確定D個地震數據塊中
所有地震數據對應的角頻率采樣間隔△ω;

其中,角頻率采樣間隔計算公式為:△ω=2π/(Gdt),G為滿足快速傅里葉變換算
法的時間采樣間隔,dt為地震數據的時間采樣間隔;

根據角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角頻率值
ωmax,確定每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,包括:

根據目標低角頻率值ωmin和角頻率采樣間隔△ω,采用起算角頻率計算公式計算
每個地震數據道對應的起算角頻率值ωp;

根據目標高角頻率值ωmax和角頻率采樣間隔△ω,采用終止角頻率計算公式計算
每個地震數據道對應的終止角頻率值ωz;

其中,起算角頻率計算公式為:終止角頻率計算公式為:
△f為頻率增量,ωmax+2π△f為通過反Q濾波處理后期望達到的角頻率的上限;

根據GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,計算每個
地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,包括:

根據GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,采用地震
幅值計算公式計算每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A;

其中,地震幅值計算公式為:
A為每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且A為復數,k為每
個頻率點的序號,△ω為角頻率采樣間隔,τj為時間深度,為時間深度τj處的Q值;

根據每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,計算
每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A,包括:

根據每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,采用
幅值和計算公式計算每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值
∑A;

其中,幅值和計算公式為: <mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> <mo>}</mo> <mo>,</mo> </mrow>
A為每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且A為復數,k為每
個頻率點的序號,△ω為角頻率采樣間隔,τj為時間深度,為時間深度τj處的Q值,Re表示
取實部。

進一步地,在所述根據所有時間深度τ處的Q值Qτ對所述D個地震數據塊中的每個
地震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到所述D個處理后的地
震數據塊之前,所述方法還包括:

接收所述CPU發送的目標低角頻率值ωmin和目標高角頻率值ωmax。

進一步地,在步驟202之后,該方法還可以包括:

對處理后的疊后地震數據體進行圖形顯示。

綜上所述,本發明實施例提供的疊后地震數據體處理方法,通過接收CPU發送的地
震工區的Q值場和疊后地震數據體,根據地震工區的Q值場對疊后地震數據體進行反Q濾波
處理,得到處理后的疊后地震數據體,由于地震工區的Q值場是CPU根據更新后的李慶忠經
驗公式計算得到的,且更新后的李慶忠經驗公式是CPU針對疊后地震數據體對李慶忠經驗
公式進行更新得到的,因此,解決了對疊后地震數據體處理的靈活性較低、可靠性較差的問
題,達到了提高疊后地震數據體處理的靈活性和可靠性的效果。

進一步地,通過CPU對疊后地震數據體進行分塊處理得到地震數據塊,GPU根據地
震數據塊對疊后地震數據體進行處理,由于GPU的處理能力較強,本發明實施例通過采用
GPU對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效
率。

請參考圖3-1,其示出了本發明實施例提供的再一種疊后地震數據體處理方法的
方法流程圖,該疊后地震數據體處理方法可以用于CPU-GPU平臺,參見圖3-1,該方法流程具
體包括:

在步驟301中,CPU從地震工區的疊后地震數據體對應的地震測線中選取M條樣本
測線。

其中,M為大于或者等于2的整數,且M=0.003N,N為疊后地震數據體的地震主測線
的條數,每條地震測線對應一個地震剖面。在本發明實施例中,當N大于667時,M=0.003N,
當N小于667時,M等于2。

示例地,在選取樣本測線時,可以從地震工區的疊后地震數據體對應的地震測線
中選取存在垂直地震剖面(英文:Vertical Seismic Profiling;簡稱:VSP)測井的地震測
線作為樣本測線,該樣本測線可以包括:地震主測線和聯絡測線中的至少一種,需要說明的
是,該樣本測線也可以是疊后地震數據體對應的地震測線中除地震主測線和聯絡測線中之
外的其他測線,且M條樣本測線中任意相鄰的兩條樣本測線之間的距離可以相等,也可以不
相等,本發明實施例對此不做限定。

可選地,在了便于樣本測線的選擇,可以先根據大地坐標,將地震工區的所有VSP
測井投射到工區平面圖上,然后根據工區平面圖,從地震工區的疊后地震數據體對應的地
震測線中選取M條存在垂直地震剖面VSP測井的地震測線作為M條樣本測線。其中,將將地震
工區的所有VSP測井投射到工區平面圖上可以便于從宏觀上了解地震工區,便于樣本測線
的選擇。

在步驟302中,CPU從地震工區的疊加速度場中確定M條樣本測線中的每條樣本測
線對應的速度場,得到M條樣本測線對應的樣本速度場ν。

在本發明實施例中,每條地震測線可以對應一個速度場,該速度場中包括該速度
場對應的地震測線對應的地震剖面上的各個點對應的傳播速度,地震工區的疊加速度場中
包括所有地震測線對應的速度場,CPU可以直接從地震工區的疊加速度場中確定M條樣本測
線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場ν,本發明實施例
對此不做限定。

在步驟303中,CPU根據M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公式確
定M條樣本測線對應的初始Q值場Q0。

其中,李慶忠經驗公式為Q=α·νβ,α=14,β=2.2。

得到M條樣本測線對應的樣本速度場ν后,可以將該M條樣本測線對應的樣本速度
場ν代入李慶忠經驗公式Q=α·νβ,求出M條樣本測線對應的初始Q值場Q0。

在步驟304中,CPU對M條樣本測線對應的初始Q值場Q0進行校正,得到M條樣本測線
對應的校正Q值場Q1。

得到M條樣本測線對應的初始Q值場Q0后,CPU對M條樣本測線對應的初始Q值場Q0進
行校正,得到M條樣本測線對應的校正Q值場Q1。

可選地,在對M條樣本測線對應的初始Q值場Q0進行校正時,CPU可以先采用譜比法
確定M條樣本測線對應的層Q值場Qi,i為垂直地震剖面VSP測井中的地層序號,然后根據M條
樣本測線對應的層Q值場Qi,確定M條樣本測線對應的時間深度的疊加Q值場Qvsp,進而根據M
條樣本測線對應的時間深度的疊加Q值場Qvsp,對M條樣本測線對應的初始Q值場Q0進行校
正,得到M條樣本測線對應的校正Q值場Q1。

在本發明實施例中,CPU采用譜比法確定M條樣本測線對應的層Q值場Qi可以包括:CPU根
據層Q值場公式確定M條樣本測線對應的層Q值場Qi,其中,層Q值場為:f
為M條樣本測線中地層的頻譜中的頻率值,ai(f)為M條樣本測線中第i層地層的頻譜中頻率
f對應的地震幅值,ai-1(f)為M條樣本測線中第i-1層地層的頻譜中頻率f對應的地震幅值,
const表示常數,δt為地震波到達第i層地層的時間與地震波到達第i-1層地層的時間的差
值。其中,f、ai(f)、ai-1(f)和δt都可以預先獲知,且CPU獲知f、ai(f)、ai-1(f)和δt的過程可以
參考現有技術,本發明實施例在此不再贅述。

在本發明實施例中,CPU然后根據M條樣本測線對應的層Q值場Qi,確定M條樣本測
線對應的時間深度的疊加Q值場Qvsp可以包括:CPU根據M條樣本測線對應的層Q值場Qi,利用
疊加Q值場公式確定M條樣本測線對應的時間深度的疊加Q值場Qvsp,其中,疊加Q值場公式
為:i為垂直地震剖面VSP測井中的地層序號,△Ti為M條樣
本測線對應的地震波在第i層地層中的單程旅行時間,n為地層總數,T為M條樣本測線對應
的地震波從地表垂直傳播至第i層地層的時間。在本發明實施例中,i、△Ti、T和n的獲取過
程均可以參考現有技術,本發明實施例在此不再贅述。

在本發明實施例中,CPU根據M條樣本測線對應的時間深度的疊加Q值場Qvsp,對M條
樣本測線對應的初始Q值場Q0進行校正,得到M條樣本測線對應的校正Q值場Q1的過程可以包
括:CPU將M條樣本測線中的每條樣本測線上位于VSP井位置(井口大地坐標為(x,y))的各時
間深度點的初始Q值與該時間深度點對應的Qvsp逐一作相除運算獲得校正系數η(x,y,τ)=
Qvsp(τ)/Q0(τ),其中,Q0(τ)是時間深度τ處的初始Q值,Qvsp(τ)是時間深度τ處由VSP井獲得的
Q值,η(x,y,τ)是大地坐標為(x,y)且時間深度為τ時的時間深度點對應的校正系數,利用校
正系數η(x,y,τ)在其大地坐標所在的樣本測線剖面內進行線性插值并做二維平滑處理,得
到每條樣本測線對應的校正系數場,進而根據M條樣本測線中的所有樣本測線對應的校正
系數場得到M條樣本測線對應的校正系數場,在M條樣本測線對應的地震剖面上逐時間深度
點將校正系數場中的校正系數與初始Q值場Q0中的Q值進行相乘運算,得到校正后的Q值場
Q1。

在步驟305中,CPU根據預設的加權系數組對M條樣本測線對應的校正Q值場Q1進行
處理,得到M條樣本測線對應的加權Q值場組。

其中,預設的加權系數組中包括P個加權系數λ,加權Q值場組中包括P套加權Q值
場,P為大于或者等于1的整數。示例地,在本發明實施例中,P等于5,加權系數組中的加權系
數可以包括:0.8、0.9、1.0、1.1、1.2,本發明實施例對此不做限定。

在本發明實施例中,CPU根據預設的加權系數組對M條樣本測線對應的校正Q值場
Q1進行處理,得到M條樣本測線對應的加權Q值場組可以包括:CPU將加權系數組中的每個加
權系數都與M條樣本測線對應的校正Q值場Q1相乘,得到M條樣本測線對應的加權Q值場組。
示例地,當加權系數可以包括:0.8、0.9、1.0、1.1、1.2時,加權Q值場組中的加權Q值場包括:
0.8Q1、0.9Q1、1.0Q1、1.1Q1和1.2Q1。

在步驟306中,CPU根據加權Q值場組對M條樣本測線組成的樣本地震數據體進行反
Q濾波處理,得到反Q濾波數據體組。

得到M條樣本測線對應的加權Q值場組后,CPU可以根據M條樣本測線對應的加權Q
值場組對M條樣本測線組成的樣本地震數據體進行反Q濾波處理,得到反Q濾波數據體組。其
中,反Q濾波數據體組中包括P套反Q濾波數據體,每套反Q濾波數據體包括M條反Q樣本測線,
每條反Q樣本測線對應一個地震剖面。

在本發明實施例中,CPU根據加權Q值場組對M條樣本測線組成的樣本地震數據體
進行反Q濾波處理的過程可以包括:CPU采用公式對M條樣本測線組成
的樣本地震數據體進行處理,其中,在該公式中,f(t)表示地震幅值,
F(ω)表示地震數據的傅里葉變換結果,ω表示角頻率,τ表示時間深度,Qq表示加權Q值場
組中的Q值,表示傅里葉逆變換,本發明實施例在此不再贅述。

在步驟307中,CPU采用傅里葉變換對反Q濾波數據體組中的每套反Q濾波數據體中
的每條反Q樣本測線對應的地震剖面進行頻譜分析。

CPU得到反Q濾波數據體組后,可以對反Q濾波數據體組中的每套反Q濾波數據體中
的每條反Q樣本測線對應的地震剖面進行頻譜分析。示例地,請參考圖3-2,其示例出的是圖
3-1所示實施例提供的一種采用傅里葉變換對每條反Q樣本測線對應的地震剖面進行頻譜
分析的方法流程圖,參見圖3-2,該方法流程可以包括:

在子步驟3071中,在反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本
測線對應的地震剖面上選擇三個矩形時窗。

在本發明實施例中,反Q濾波數據體組中包括P套反Q濾波數據體,每套反Q濾波數
據體包括M條反Q樣本測線,每條反Q樣本測線對應一個地震剖面。CPU可以在每條飯Q樣本測
線對應的地震剖面上選擇三個矩形時窗,其中,三個矩形時窗包括:淺層矩形時窗、中層矩
形時窗和深層矩形時窗。

需要說明的是,在本發明實施例中,CPU在選擇矩形時窗時,選擇的三個矩形時窗
的大小相等,且三個矩形時窗無重疊區域,每個矩形時窗內的同相軸的清晰度大于預設清
晰度,同相軸的連續度大于預設連續度。預設清晰度用于區分矩形時窗內的同相軸是否清
晰,當矩形時窗內的同相軸的清晰度大于或者等于預設清晰度時,表示矩形時窗內的同相
軸較為清晰,當矩形時窗內的同相軸的清晰度小于預設清晰度時,表示矩形時窗內的同相
軸的不清晰;預設連續度用于分矩形時窗內的同相軸是否連續,當矩形時窗內的同相軸的
連續度大于或者等于預設連續度時,表示矩形時窗內的同相軸連續,當矩形時窗內的同相
軸的連續度小于預設連續度時,表示矩形時窗內的同相軸的不連續,其中,預設清晰度和預
設連續度的具體數值均可以根據實際需要進行設置,本發明實施例對此不做限定。

在子步驟3072中,采用傅里葉變換對每條反Q樣本測線對應的地震剖面上位于三
個矩形時窗中的每個矩形時窗內的部分進行處理,得到每條反Q樣本測線對應的3個反Q樣
本頻譜。

在本發明實施例中,每條反Q樣本測線對應的地震剖面上的三個矩形時窗中的每
個矩形時窗中可以包括每條反Q樣本測線對應的地震剖面的一部分,為了簡化分析過程,
CPU可以采用傅里葉變換對每條反Q樣本測線對應的地震剖面上位于每個矩形時窗內的部
分進行處理,得到每條反Q樣本測線對應的3個反Q樣本頻譜。其中,CPU采用傅里葉變換對每
條反Q樣本測線對應的地震剖面上位于三個矩形時窗中的每個矩形時窗內的部分進行處
理,得到每條反Q樣本測線對應的3個反Q樣本頻譜的過程可以參考現有技術,本發明實施例
在此不再贅述。

在子步驟3073中,根據每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻
譜,得到每套反Q濾波數據體對應的3M個反Q樣本頻譜。

CPU采用傅里葉變換對每條反Q樣本測線對應的地震剖面上位于三個矩形時窗中
的每個矩形時窗內的部分進行處理,得到每條反Q樣本測線對應的3個反Q樣本頻譜后,由于
每套反Q濾波數據體中包括M條反Q樣本測線,因此,CPU可以根據每套反Q濾波數據體中的M
條反Q樣本測線對應的反Q樣本頻譜,得到每套反Q濾波數據體對應的3M個反Q樣本頻譜。示
例地,當M=2時,3M=6,也即是,每套反Q濾波數據體對應的6個反Q樣本頻譜。

在子步驟3074中,根據反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻
譜得到3×M×P個反Q樣本頻譜。

CPU得到每套反Q濾波數據體對應的3M個反Q樣本頻譜后,由于反Q濾波數據體組中
包括P套反Q濾波數據體,因此,CPU可以根據反Q濾波數據體組中的P套反Q濾波數據體對應
的反Q樣本頻譜,得到3×M×P個反Q樣本頻譜。其中,3×M×P個反Q樣本頻譜中包括:M×P個
淺層矩形時窗對應的反Q樣本頻譜、M×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層
矩形時窗對應的反Q樣本頻譜。示例地,當M=2,P=5時,3×M×P=30,也即是,CPU根據反Q
濾波數據體組中的5套反Q濾波數據體對應的反Q樣本頻譜得到30個反Q樣本頻譜,該30個反
Q樣本頻譜中包括:10個淺層矩形時窗對應的反Q樣本頻譜、10個中層矩形時窗對應的反Q樣
本頻譜和10個深層矩形時窗對應的反Q樣本頻譜,本發明實施例在此不再贅述。

在子步驟3075中,對3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。

CPU得到3×M×P個反Q樣本頻譜后,可以對該3×M×P個反Q樣本頻譜中的每個反Q
樣本頻譜進行頻譜分析,示例地,參見子步驟3074可知,3×M×P=30,因此,CPU對該30個反
Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。其中,CPU對每個反Q樣本頻譜進行頻譜分
析的過程可以參考現有技術,本發明實施例在此不再贅述。

在步驟308中,CPU根據頻譜分析結果進行處理,確定M條反Q樣本測線對應的加權
系數場。

CPU采用傅里葉變換對反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的地震剖面進行頻譜分析后,可以根據頻譜分析結果進行處理,確定M條反Q樣
本測線對應的加權系數場。示例地,請參考圖3-3,其示例出的是圖3-1所示實施例提供的一
種根據頻譜分析結果進行處理確定M條反Q樣本測線對應的加權系數場的方法流程圖,參見
圖3-3,該方法流程可以包括:

在子步驟3081中,根據頻譜分析結果確定反Q濾波數據體組中的每套反Q濾波數據
體中的每條反Q樣本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜。

在本發明實施例中,參見子步驟3074可知,每條反Q樣本測線對應的地震剖面上位
于淺層矩形時窗內的部分對應P個反Q樣本頻譜,位于中層矩形時窗內的部分對應P個反Q樣
本頻譜,位于深層矩形時窗內的部分對應P個反Q樣本頻譜。CPU可以根據頻譜分析結果確定
反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的第一目標頻譜、第
二目標頻譜和第三目標頻譜。其中,第一目標頻譜為每條反Q樣本測線對應的地震剖面上位
于淺層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成的封閉區域的面積最大的
反Q樣本頻譜,第二目標頻譜為每條反Q樣本測線對應的地震剖面上位于中層矩形時窗內的
部分對應的P個反Q樣本頻譜中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,第三目
標頻譜為每條反Q樣本測線對應的地震剖面上位于深層矩形時窗內的部分對應的P個反Q樣
本頻譜中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜,該第一目標頻譜、第二目標頻
譜和第三目標頻譜都可以稱為恰當的補償頻譜,本發明實施例對此不做限定。

在子步驟3082中,分別確定第一目標頻譜對應的矩形時窗對應的第一目標加權系
數、第二目標頻譜對應的矩形時窗對應的第二目標加權系數和第三目標頻譜對應的矩形時
窗對應的第三目標加權系數,得到每條反Q樣本測線對應的3個目標加權系數。

CPU確定每條反Q樣本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜
后,可以分別確定第一目標頻譜對應的矩形時窗對應的第一目標加權系數、第二目標頻譜
對應的矩形時窗對應的第二目標加權系數和第三目標頻譜對應的矩形時窗對應的第三目
標加權系數,得到每條反Q樣本測線對應的3個目標加權系數。其中,第一目標加權系數、第
二目標加權系數和第三目標加權系數都為預設的加權系數組中的加權系數。

需要說明的是,由于第一目標頻譜、第二目標頻譜和第三目標頻譜都是反Q濾波數
據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的反Q樣本頻譜,而反Q濾波數據
體組是CPU根據M條樣本測線對應的加權Q值場組對M條樣本測線組成的樣本地震數據體進
行反Q濾波處理得到的,因此,CPU可以反推出第一目標頻譜對應的矩形時窗對應的第一目
標加權系數、第二目標頻譜對應的矩形時窗對應的第二目標加權系數和第三目標頻譜對應
的矩形時窗對應的第三目標加權系數,本發明實施例在此不再贅述。

在子步驟3083中,根據M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權
系數,得到M條反Q樣本測線對應的3M個目標加權系數。

由于每條反Q樣本測線對應的3個目標加權系數,因此,M條樣本測線可以對應3M個
目標加權系數,CPU可以直接根據M條反Q樣本測線中的每條煩Q樣本測線對應的3個目標加
權系數得到3M個目標加權系數,本發明實施在此不再贅述。

在子步驟3084中,對3M個目標加權系數進行空間插值,得到地震工區對應的加權
系數場。

得到M條反Q樣本測線對應的3M個目標加權系數后,CPU可以對M條反Q樣本測線對
應的3M個目標加權系數進行空間插值,得到地震工區對應的加權系數場。

在本發明實施例中,CPU對3M個目標加權系數進行空間插值的過程可以包括:對于
每條反Q樣本測線而言,假設該每條反Q樣本測線對應的地震剖面上淺層矩形時窗的中心點
對應的大地坐標是(x1,y1),在時間深度τ處,該淺層矩形時窗的中心點對應的目標加權系
數為λ1,假設該每條反Q樣本測線對應的地震剖面上中層矩形時窗的中心點對應的大地坐
標是(x2,y2),在時間深度τ處,該淺層矩形時窗的中心點對應的目標加權系數為λ2,假設該
每條反Q樣本測線對應的地震剖面上深層矩形時窗的中心點對應的大地坐標是(x3,y3),在
時間深度τ處,該淺層矩形時窗的中心點對應的目標加權系數為λ3,可以得到3個空間加權系
數點,也即,每條反Q樣本測線對應3個空間加權系數點,因此,M條反Q樣本測線可以對應3M
個空間加權系數點,該3M個空間加權系數點都為已知的空間加權系數點,可以采用空間三
維插值方法對該3M個已知的空間加權系數點進行空間插值,得到地震工區對應的加權系數
場。其中,采用空間三維插值方法對該3M個已知的空間加權系數點進行空間插值的過程可
以參考現有技術,本發明實施例在此不再贅述。

在子步驟3085中,對地震工區對應的加權系數場進行平滑處理,得到處理后的加
權系數場。

由于在子步驟3084中,CPU對3M個目標加權系數進行空間插值得到地震工區對應
的加權系數場是離散的加權系數場,該離散的加權系數場中包括的是多個離散的空間加權
系數點,因此,CPU可以對該離散的加權系數場進行平滑處理,得到處理后的加權系數場。示
例地,CPU對地震工區對應的加權系數場做空間三維平滑處理,得到處理后的加權系數場,
本發明實施例在此不再贅述。

在子步驟3086中,從處理后的加權系數場中抽取出M條反Q樣本測線對應的加權系
數場。

CPU對地震工區對應的加權系數場進行平滑處理得到處理后的加權系數場后,可
以從處理后的加權系數場中抽取出M條反Q樣本測線對應的加權系數場,示例地,CPU根據M
條樣本測線對應的空間坐標值,從處理后的加權系數場中抽取出M條反Q樣本測線對應的加
權系數場。其中,CPU根據M條樣本測線對應的空間坐標值從處理后的加權系數場中抽取出M
條反Q樣本測線對應的加權系數場的過程可以參考現有技術,本發明實施例在此不再贅述。

在步驟309中,CPU將M條反Q樣本測線對應的加權系數場確定為M條樣本測線對應
的加權系數場λfinal。

抽取出M條反Q樣本測線對應的加權系數場后,CPU可以將M條反Q樣本測線對應的
加權系數場確定為M條樣本測線對應的加權系數場λfinal,本發明實施例在此不再贅述。

在步驟310中,CPU根據M條樣本測線對應的校正Q值場Q1和M條樣本測線對應的加
權系數場λfinal,確定M條樣本測線對應的最終Q值場Q2。

由于在步驟304中,CPU得到了M條樣本測線對應的校正Q值場Q1,在步驟309中CPU
得到了M條樣本測線對應的加權系數場λfinal,因此,CPU可以根據步驟304中得到的M條樣本
測線對應的校正Q值場Q1和步驟309中得到的M條樣本測線對應的加權系數場λfinal確定M條
樣本測線對應的最終Q值場Q2。

示例地,在本發明實施例中,CPU可以將M條樣本測線對應的校正Q值場Q1和M條樣
本測線對應的加權系數場λfinal相乘,將M條樣本測線對應的校正Q值場Q1與M條樣本測線對
應的加權系數場λfinal的乘積作為M條樣本測線對應的最終Q值場Q2,本發明實施例對此不做
限定。

在步驟311中,CPU根據M條樣本測線對應的樣本速度場ν和M條樣本測線對應的最
終Q值場Q2,采用最小二乘法更新李慶忠經驗公式,得到更新后的李慶忠經驗公式。

由于在步驟302中CPU得到了M條樣本測線對應的樣本速度場ν,在步驟310中CPU得
到了M條樣本測線對應的最終Q值場Q2,因此,CPU可以根據M條樣本測線對應的樣本速度場ν
和M條樣本測線對應的最終Q值場Q2采用最小二乘法更新李慶忠經驗公式,得到更新后的李
慶忠經驗公式。

其中,CPU根據M條樣本測線對應的樣本速度場ν和M條樣本測線對應的最終Q值場
Q2采用最小二乘法更新李慶忠經驗公式得到更新后的李慶忠經驗公式可以包括:CPU根據M
條樣本測線對應的樣本速度場ν和M條樣本測線對應的最終Q值場Q2,采用最小二乘法確定
李慶忠經驗公式Q=α·νβ中的α值與β值,得到更新α值和更新β值,然后將更新α值和更新β
值代入李慶忠經驗公式,得到更新后的李慶忠經驗公式。具體地,CPU將李慶忠經驗公式Q=
α·νβ中的α與β都看做未知數,將M條樣本測線對應的樣本速度場ν中的速度值和M條樣本測
線對應的最終Q值場Q2中的Q值代入李慶忠經驗公式Q=α·νβ,計算得到更新α值和更新β值,
然后將Q和ν看做未知數,將更新α值和更新β值代入李慶忠經驗公式,得到更新后的李慶忠
經驗公式,本發明實施例在此不再贅述。

在步驟312中,CPU根據地震工區的疊加速度場,采用更新后的李慶忠經驗公式確
定地震工區的Q值場Qall。

得到更新后的李慶忠經驗公式后,CPU可以根據地震工區的疊加速度場,采用更新
后的李慶忠經驗公式確定地震工區的Q值場Qall,示例地,CPU將地震工區的疊加速度場代入
更新后的李慶忠經驗公式,得到地震工區的Q值場Qall,本發明實施例在此不再贅述。

在步驟313中,按照預設采樣間隔對地震工區的Q值場Qall進行采樣,得到采樣Q值
場Qs。

CPU得到地震工區的Q值場Qall后,可以按照預設采樣間隔對地震工區的Q值場Qall
進行采樣,得到采樣Q值場Qs,其中,預設采樣間隔可以根據實際需要進行設置,本發明實施
例對此不做限定,示例地,預設采樣間隔為20ms(中文:毫秒),CPU對地震工區的Q值場Qall進
行采樣的過程可以參考現有技術,本發明實施例在此不再贅述。

在步驟314中,對疊后地震數據體進行分塊處理,得到D個地震數據塊,D為大于或
者等于2的整數。

由于疊后地震數據體中的數據量較大,因此,CPU可以對疊后地震數據體進行分塊
處理,得到D個地震數據塊,其中,D為大于或者等于2的整數,且D的具體數值可以根據實際
情況確定,本發明實施例對此不做限定。

示例地,請參考圖3-4,其示例出的是圖3-1所示實施例提供的一種對疊后地震數
據體進行分塊處理的方法流程圖,參見圖3-4,該方法流程可以包括:

在子步驟3141中,計算疊后地震數據體中的地震數據的總道數E。

其中,CPU計算疊后地震數據體中的地震數據的總道數E具體實現過程可以參考現
有技術,本發明實施例在此不再贅述。

在子步驟3142中,根據GPU在一個周期內處理的地震數據的道數Tracemax和疊后地
震數據體中的地震數據的總道數E,對疊后地震數據體進行分塊處理,得到D個地震數據塊。

得到疊后地震數據體中的地震數據的總道數E后,CPU可以根據GPU在一個周期內
處理的地震數據的道數Tracemax和疊后地震數據體中的地震數據的總道數E,對疊后地震數
據體進行分塊處理,得到D個地震數據塊,其中,E=Tracemax×(D-1)+L,L為D個地震數據塊
中的最后一個地震數據塊中的地震數據的道數,其中,Tracemax的取值為2u,256≤2u≤
0.001E,且u為正整數,本發明實施例對此不做限定。

在步驟315中,CPU向GPU發送地震工區的Q值場Qall和疊后地震數據體。

CPU得到地震工區的Q值場Qall后,可以向GPU發送地震工區的Q值場Qall和疊后地震
數據體,以便于GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處
理后的疊后地震數據體。

需要說明的是,由于在步驟314中CPU對疊后地震數據體進行分塊處理得到了D個
地震數據塊,且在步驟313中CPU按照預設采樣間隔對地震工區的Q值場Qall進行采樣得到采
樣Q值場Qs,因此,CPU向GPU發送地震工區的Q值場Qall和疊后地震數據體也即是CPU依次向
GPU發送采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,本發明實施例對此不作限定。

還需要說明的是,為了便于對疊后地震數據體進行處理,CPU在向GPU發送疊后地
震數據體之前,可以對疊后地震數據體進行去噪處理,示例地,CPU對疊后地震數據體進行
保幅去噪處理,本發明實施例對此不做限定。

還需要說明的是,在向GPU發送地震工區的Q值場Qall和疊后地震數據體之前,CPU
可以采用傅里葉變換對M條樣本測線中的每條樣本測線對應的地震剖面進行頻譜分析,得
到M條樣本測線對應的目標低角頻率值ωmin和目標高角頻率值ωmax,CPU還可以向GPU發送
該目標低角頻率值ωmin和目標高角頻率值ωmax。其中,CPU采用傅里葉變換對M條樣本測線
中的每條樣本測線對應的地震剖面進行頻譜分析,得到M條樣本測線對應的目標低角頻率
值ωmin和目標高角頻率值ωmax的過程可以包括:CPU在疊后地震數據體對應的M條樣本測線
中的每條樣本測線對應的地震剖面上選擇三個矩形時窗,然后采用傅里葉變換對每條樣本
測線對應的地震剖面上位于三個矩形時窗中的每個矩形時窗內的部分進行處理,得到每條
樣本測線對應的3個樣本頻譜,根據M條樣本測線中的每條樣本測線對應的3個樣本頻譜,得
到M條樣本測線對應的3M個樣本頻譜,進而確定3M個樣本頻譜中的每個樣本頻譜對應的最
高角頻率值和最低角頻率值,得到3M個最高角頻率值和3M個最低角頻率值,然后將3M個最
低角頻率值中最小的角頻率值作為M條樣本測線對應的目標低角頻率值ωmin,將3M個最高
角頻率值中最大的角頻率值作為M條樣本測線對應的目標高角頻率值ωmax。其中,該三個矩
形時窗包括:淺層矩形時窗、中層矩形時窗和深層矩形時窗,在本發明實施例中,CPU在選擇
矩形時窗時,選擇的三個矩形時窗的大小相等,且三個矩形時窗無重疊區域,每個矩形時窗
內的同相軸的清晰度大于預設清晰度,同相軸的連續度大于預設連續度。預設清晰度用于
區分矩形時窗內的同相軸是否清晰,當矩形時窗內的同相軸的清晰度大于或者等于預設清
晰度時,表示矩形時窗內的同相軸較為清晰,當矩形時窗內的同相軸的清晰度小于預設清
晰度時,表示矩形時窗內的同相軸的不清晰;預設連續度用于分矩形時窗內的同相軸是否
連續,當矩形時窗內的同相軸的連續度大于或者等于預設連續度時,表示矩形時窗內的同
相軸連續,當矩形時窗內的同相軸的連續度小于預設連續度時,表示矩形時窗內的同相軸
的不連續,其中,預設清晰度和預設連續度的具體數值均可以根據實際需要進行設置,本發
明實施例對此不做限定。

在步驟316中,GPU接收CPU發送的地震工區的Q值場Qall和疊后地震數據體。

CPU向GPU發送地震工區的Q值場Qall和疊后地震數據體時,GPU可以接收CPU發送的
地震工區的Q值場Qall和疊后地震數據體。在本發明實施例中,由于CPU向GPU發送地震工區
的Q值場Qall和疊后地震數據體可以包括CPU依次向GPU發送采樣Q值場Qs和D個地震數據塊
中的每個地震數據塊。因此,GPU接收CPU發送的地震工區的Q值場Qall和疊后地震數據體可
以包括GPU接收CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,其中,采樣Q
值場Qs是CPU按照預設采樣間隔對地震工區的Q值場Qall進行采樣得到的,D個地震數據塊是
CPU對疊后地震數據體進行分塊處理得到的。

在步驟317中,GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,
得到處理后的疊后地震數據體。

接收到CPU發送的地震工區的Q值場Qall和疊后地震數據體后,GPU可以根據地震工
區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理后的疊后地震數據體。示例
地,請參考圖3-5,其示例出的是圖3-1所示實施例提供的一種根據地震工區的Q值場對疊后
地震數據體進行反Q濾波處理的方法流程圖,參見圖3-5,該方法流程可以包括:

在子步驟3171中,根據采樣Q值場Qs對D個地震數據塊中的每個地震數據塊進行反
Q濾波處理,得到D個處理后的地震數據塊。

由于GPU接收到的實際上是地震數據塊和采樣Q值場Qs,因此,GPU根據地震工區的
Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理后的疊后地震數據體也即是GPU根
據采樣Q值場Qs對D個地震數據塊中的每個地震數據塊進行反Q濾波處理,得到D個處理后的
地震數據塊,進而根據D個處理后的地震數據塊處理得到處理后的疊后地震數據體。

示例地,請參考圖3-6,其示例出的是圖3-1所示實施例提供的一種根據采樣Q值場
對D個地震數據塊中的每個地震數據塊進行反Q濾波處理的方法流程圖,參見圖3-6,該方法
流程可以包括:

在子步驟31711中,對采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ。

由于采樣Q值場Qs是CPU按照預設采樣間隔對地震工區的Q值場Qall進行采樣得到
的,因此,采樣Q值場Qs中的Q值較少,因此,GPU可以對采樣Q值場Qs進行處理得到每個時間深
度τ處的Q值Qτ。

其中,GPU對采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ可以包括:GPU根
據采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ。其中,Q值插值計算公式
為:

<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>

<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>τj、τ(j-1)和τ(j+1)都為時間深度,
為時間深度τj處的Q值插值,為時間深度τ(j-1)處的Q值插值,為時間深度τ(j+1)
處的Q值插值,表示對取整。

在子步驟31712中,根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個地震
數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到D個處理后的地震數據塊。

得到每個時間深度τ處的Q值Qτ后,GPU可以根據所有時間深度τ處的Q值Qτ對D個地
震數據塊中的每個地震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到D
個處理后的地震數據塊。示例地,請參考圖3-7,其示例出的是圖3-1所示實施例提供的一種
根據所有時間深度處的Q值對D個地震數據塊中的每個地震數據塊中的所有時間深度處的
地震幅值進行反Q濾波處理的方法流程圖,參見圖3-7,該方法流程可以包括:

在子步驟317121中,根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個地
震數據塊中的每個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處理
后的地震數據道。

在本發明實施例中,每個地震數據道包括多個時間深度τ,每個時間深度τ對應一
個Q值Qτ和一個地震幅值,因此,GPU可以根據D個地震數據塊中的每個地震數據塊中的每個
地震數據道的每個時間深度τ處的Q值Qτ對該每個時間深度τ處的地震幅值進行反Q濾波處
理,進而根據每個地震數據道的所有時間深度τ處的Q值Qτ對該每個地震數據道的所有時間
深度τ處的地震幅值進行反Q濾波處理。

示例地,請參考圖3-8,其示例出的是圖3-1所示實施例提供的一種根據所有時間
深度處的Q值對D個地震數據塊中的每個地震數據塊中的每個地震數據道的所有時間深度
處的地震幅值進行反Q濾波處理的方法流程圖,參見圖3-8,該方法流程可以包括:

在子步驟3171211中,分別確定GPU中的線程的線程格尺寸和線程塊尺寸。

GPU中可以包括多個線程,GPU可以確定GPU中的線程的線程格尺寸和線程塊尺寸。
其中,GPU中的線程的線程格尺寸為:Tracemax為GPU在一個周期
內處理的地震數據的道數,b_size=2m,m為正整數,且4≤m≤6,S為地震數據道的采樣點
數,GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,m為正整數,且4≤m≤6,
本發明實施例對此不做限定。

在子步驟3171212中,對D個地震數據塊中的每個地震數據塊中的每個地震數據道
進行快速傅里葉變換。

其中,GPU對D個地震數據塊中的每個地震數據塊中的每個地震數據道進行快速傅
里葉變換的過程可以參考現有技術,本發明實施例在此不再贅述。

在子步驟3171213中,根據快速傅里葉變換結果,確定D個地震數據塊中所有地震
數據對應的角頻率采樣間隔△ω。

對D個地震數據塊中的每個地震數據塊中的每個地震數據道進行快速傅里葉變換
后,GPU可以根據快速傅里葉變換結果,確定D個地震數據塊中所有地震數據對應的角頻率
采樣間隔△ω。

其中,GPU根據快速傅里葉變換結果,確定D個地震數據塊中所有地震數據對應的
角頻率采樣間隔△ω可以包括:GPU根據快速傅里葉變換結果,采用角頻率采樣間隔計算公
式確定D個地震數據塊中所有地震數據對應的角頻率采樣間隔△ω。其中,角頻率采樣間隔
計算公式為:△ω=2π/(Gdt),G為滿足快速傅里葉變換算法的時間采樣間隔,dt為地震數
據的時間采樣間隔,G和dt的求取過程可以參考現有技術,本發明實施例在此不再贅述。

在子步驟3171214中,根據角頻率采樣間隔△ω和預先接收的目標低角頻率值
ωmin、目標高角頻率值ωmax,確定每個地震數據道對應的起算角頻率值ωp和終止角頻率值
ωz。

得到D個地震數據塊中所有地震數據對應的角頻率采樣間隔△ω后,GPU可以根據
角頻率采樣間隔△ω和預先接收的目標低角頻率值ωmin、目標高角頻率值ωmax,確定每個
地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,其中,GPU根據角頻率采樣間隔△
ω和預先接收的目標低角頻率值ωmin、目標高角頻率值ωmax,確定每個地震數據道對應的
起算角頻率值ωp和終止角頻率值ωz的過程可以包括:GPU根據目標低角頻率值ωmin和角頻
率采樣間隔△ω,采用起算角頻率計算公式計算每個地震數據道對應的起算角頻率值ωp,
根據目標高角頻率值ωmax和角頻率采樣間隔△ω,采用終止角頻率計算公式計算每個地震
數據道對應的終止角頻率值ωz,其中,起算角頻率計算公式為:終止角頻
率計算公式為:△f為頻率增量,ωmax+2π△f為通過反Q濾波處理
后期望達到的角頻率的上限,△f和ωmax+2π△f的求取過程可以參考現有技術,本發明實施
例在此不再贅述。

需要說明的是,GPU還可以計算xindex=blockDim.x*blockIdx.x+threadIdx.x,
yindex=blockDim.y*blockIdx.y+threadIdx.y,其中xindex與yindex分別為線程在線程
格x方向與y方向的索引值,blockDim.x與blockDim.y分別是線程在線程塊x方向與y方向的
尺寸值,blockIdx.x與blockIdx.y分別是線程塊在線程格中x方向與y方向的索引值,
threadIdx.x與threadIdx.y分別是線程在線程塊中x方向與y方向的索引值,本發明實施例
在此不再贅述。

在子步驟3171215中,根據GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度
τ處的Q值Qτ中的時間深度τj處的Q值每個地震數據道對應的起算角頻率值ωp和終止角
頻率值ωz,計算每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A。

得到GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ、每個地
震數據道對應的起算角頻率值ωp和終止角頻率值ωz后,GPU可以根據GPU中的線程的線程
格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間深度τj處的Q值每個地震數據
道對應的起算角頻率值ωp和終止角頻率值ωz,計算每個地震數據道對應的時間深度τj處
的每個頻率點對應的地震幅值A。示例地,GPU可以根據GPU中的線程的線程格尺寸、線程塊
尺寸、所有時間深度τ處的Q值Qτ中的時間深度τj處的Q值每個地震數據道對應的起算角
頻率值ωp和終止角頻率值ωz,采用地震幅值計算公式計算每個地震數據道對應的時間深
度τj處的每個頻率點對應的地震幅值A;其中,地震幅值計算公式為:
A為每個地震數據道對應的時
間深度τj處的每個頻率點對應的地震幅值且A為復數,k為每個頻率點的序號,△ω為角頻
率采樣間隔,τj為時間深度,為時間深度τj處的Q值。

在子步驟3171216中,根據每個地震數據道對應的時間深度τj處的每個頻率點對
應的地震幅值A,計算每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值
∑A。

得到每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A后,GPU
可以根據每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,計算每個
地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A。示例地,GPU可以根據
每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,采用幅值和計算公
式計算每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A;

其中,幅值和計算公式為: <mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> <mo>}</mo> <mo>,</mo> </mrow>
A為每個地震數據道對應的時間深度τj處的每個頻率點對應
的地震幅值且A為復數,k為每個頻率點的序號,△ω為角頻率采樣間隔,τj為時間深度,
為時間深度τj處的Q值,Re表示取實部。

在子步驟3171217中,根據GPU中的所有線程處理得到的每個地震數據道對應的所
有時間深度處的所有頻率點對應的地震幅值,處理得到處理后的地震數據道。

得到每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值∑A后,
GPU可以根據GPU中的所有線程處理得到的每個地震數據道對應的所有時間深度處的所有
頻率點對應的地震幅值,處理得到處理后的地震數據道,本發明實施例在此不再贅述。

在子步驟317122中,根據每個地震數據塊中的所有處理后的地震數據道,得到處
理后的地震數據塊。

得到每個地震數據塊中的處理后的地震數據道后,GPU可以根據每個地震數據塊
中的所有處理后的地震數據道,得到處理后的地震數據塊。

在子步驟317123中,根據D個地震數據塊中的所有處理后的地震數據塊,得到D個
處理后的地震數據塊。

得到D個地震數據塊中的每個處理后的地震數據塊后,GPU可以根據D個地震數據
塊中的所有處理后的地震數據塊,得到D個處理后的地震數據塊。

在子步驟3172中,根據D個處理后的地震數據塊,得到處理后的疊后地震數據體。

得到D個處理后的地震數據塊后,GPU可以根據D個處理后的地震數據塊,得到處理
后的疊后地震數據體。

需要說明的是,在本發明實施例中,GPU在根據所有時間深度τ處的Q值Qτ對D個地
震數據塊中的每個地震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到D
個處理后的地震數據塊之前,可以先接收CPU發送的目標低角頻率值ωmin和目標高角頻率
值ωmax;且GPU在根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理
后的疊后地震數據體之后,還可以對處理后的疊后地震數據體進行圖形顯示,其中,GPU顯
示的圖形可以指示地下構造的地層沉積樣式與微小斷裂,可以用于構造解釋與砂體追蹤,
本發明實施例在此不再贅述。

綜上所述,本發明實施例提供的疊后地震數據體處理方法,通過從地震工區的疊
后地震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣
本測線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣
本測線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條
樣本測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測
線對應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值
場;根據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法
更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更
新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數
據體,以便于GPU根據地震工區的Q值場對疊后地震數據體進行反Q濾波處理,得到處理后的
疊后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決了對
疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體處理
的靈活性和可靠性的效果。

進一步地,現有技術中,通常在CPU中對疊后地震數據體進行處理,由于CPU在一個
計算周期內只能對一個地震數據道進行處理,且其為單道輸入輸出,讀寫頻繁,處理效率較
低。本發明實施例提供的疊后地震數據體處理方法,通過CPU對疊后地震數據體進行分塊處
理得到地震數據塊,GPU根據地震數據塊對疊后地震數據體進行處理,由于GPU在一個計算
周期內可以對多個地震數據道進行處理,GPU的處理能力較強,本發明實施例通過采用GPU
對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效率。

本發明實施例提供的疊后地震數據體處理方法可以提高疊后地震數據的分辨率,
特別適用于處理大規模的三維偏移疊加地震數據體,對提高疊后地震數據體的分辨率具有
重要應用價值。

下面以一個具體的例子對本發明實施例提供的疊后地震數據體處理方法進行簡
單說明。具體地,以冀東油田柏各莊區塊為例,CPU型號為Xeon E5-2687,GPU型號為Quadro
K5000,具體的處理過程可以包括以下步驟:

(1)CPU從地震工區的疊后地震數據體對應的地震測線中選取2條樣本測線。

具體包括:CPU統計地震工區內VSP測井數為3,根據大地坐標,將地震工區的所有
VSP測井投射到工區平面圖上,統計地震工區內地震主測線條數N=500,從500條地震主測
線條中選取M=2條地震主測線作為樣本測線,且使每條樣本測線上分布一口VSP測井。

在選取樣本測線之前,CPU對疊后地震數據體進行去噪處理,其中,疊后地震數據
體的采樣間隔為2ms,記錄時長為2500ms,三維偏移疊后地震數據體中共包含300000個成像
道(地震數據道),水平地震測線之間的道間距為15m(中文:米),垂直地震測線的道間距為
40m,疊后地震數據體的覆蓋面積是180平方千米,示例地,如圖3-9所示,其示出的是垂直坐
標為8000m處的沿地震測線方向切片的三維偏移疊加剖面。其中,在該圖3-9中,橫軸表示垂
直距離LS,單位為km(中文:千米),縱軸表示時間深度τ,單位為s(中文:秒)。如圖3-10所示,
其示出的是圖3-9所示的區域390部分的放大圖。

(2)CPU從地震工區的疊加速度場中抽取出2條樣本測線中的每條樣本測線對應的
速度場,得到該2條樣本測線對應的樣本速度場ν,該樣本速度場ν內的速度值按照樣本測線
的測線號的順序依次排列。

(3)CPU將樣本速度場ν內的速度值代入李慶忠經驗公式Q=α·νβ(α=14,β=2.2)
得到該2條樣本測線對應的初始的Q值場Q0。

(4)CPU對該2條樣本測線對應的初始Q值場Q0進行校正,得到M條樣本測線對應的
校正Q值場Q1。

具體包括:采用譜比法確定該2條樣本測線對應的層Q值場Qi,i為垂直地震剖面
VSP測井中的地層序號,然后將層Q值場Qi轉換為該2條樣本測線對應的時間深度的疊加Q值
場Qvsp,進而根據該2條樣本測線對應的時間深度的疊加Q值場Qvsp對該2條樣本測線對應的
初始Q值場Q0進行校正,得到M條樣本測線對應的校正Q值場Q1。

(5)CPU根據預設的加權系數組對2條樣本測線對應的校正Q值場Q1進行處理,得到
2條樣本測線對應的加權Q值場組。

具體包括:CPU將校正Q值場Q1分別與加權系數組中的加權系數λ進行相乘(λ取值
為0.8,0.9,1.0,1.1,1.2),得到加權Q值場組,該加權Q值場組中包括5套加權Q值場。

(6)CPU根據加權Q值場組對2條樣本測線組成的樣本地震數據體進行反Q濾波處
理,得到反Q濾波數據體組,該反Q濾波數據體組中包括5套反Q濾波數據體,每套反Q濾波數
據體中包括2條反Q樣本測線。

(7)CPU采用傅里葉變換對反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q
樣本測線對應的地震剖面進行頻譜分析。

(8)CPU根據頻譜分析結果進行處理,確定2條反Q樣本測線對應的加權系數場。

具體包括:,根據頻譜分析結果確定反Q濾波數據體組中的每套反Q濾波數據體中
的每條反Q樣本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,分別確定第一目
標頻譜對應的矩形時窗對應的第一目標加權系數、第二目標頻譜對應的矩形時窗對應的第
二目標加權系數和第三目標頻譜對應的矩形時窗對應的第三目標加權系數,得到每條反Q
樣本測線對應的3個目標加權系數,根據2條反Q樣本測線中的每條反Q樣本測線對應的3個
目標加權系數,得到M條反Q樣本測線對應的6個目標加權系數,對6個目標加權系數進行空
間插值,得到地震工區對應的加權系數場,對地震工區對應的加權系數場進行平滑處理,得
到處理后的加權系數場,從處理后的加權系數場中抽取出2條反Q樣本測線對應的加權系數
場。

(9)CPU將2條反Q樣本測線對應的加權系數場確定為2條樣本測線對應的加權系數
場λfinal。

(10)CPU將步驟(9)得到的加權系數場λfinal與步驟(4)得到的校正Q值場Q1進行相
乘運算,得到2條樣本測線對應的最終Q值場Q2。

(11)CPU根據2條樣本測線對應的樣本速度場ν和2條樣本測線對應的最終Q值場
Q2,采用最小二乘法更新李慶忠經驗公式,得到更新后的李慶忠經驗公式。

(12)CPU根據地震工區的疊加速度場,采用更新后的李慶忠經驗公式確定地震工
區的Q值場Qall。

(13)CPU按照預設采樣間隔對地震工區的Q值場Qall進行采樣,得到采樣Q值場Qs。

(14)CPU對疊后地震數據體進行分塊處理,得到D個地震數據塊。

具體包括:CPU計算疊后數據體的總道數為E=300000,設GPU在一個周期內處理的
地震數據的道數Trace_max=256,將總道數為E=300000的疊后數據體按照每塊含有
Trace_max=256道進行分塊操作,共產生D=1171塊內含256道地震數據的數據塊,會剩余
不足256道的尾道部分,該尾道部分的道數為224,將該部分數據作為最后一塊數據體。

(15)CPU將得到的地震工區的Q值場Qall和地震數據塊輸入到GPU。

具體地,CPU將分塊得到的D個地震數據塊逐塊輸入GPU內。

(16)GPU根據地震工區的Q值場Qall對D個地震數據塊中的每個地震數據愉塊進行
反Q濾波處理,得到處理后的疊后地震數據體。

GPU處理完后,可以將處理結果輸出到CPU端,CPU與GPU相互協同,共同完成所有分
塊數據的輸入、反Q濾波處理與輸出操作。本發明應用CPU-GPU平臺對地震工區的所有地震
數據處理完用時618秒,而現有技術中僅在CPU上執行相同數據的反Q濾波處理用時13600
秒,通過對比可以看出,應用CPU-GPU平臺后,計算效率顯著提升,加速比達到22倍。

(17)GPU通過顯示軟件對處理后的疊后地震數據體進行圖像顯示。

GPU顯示的處理后的疊后地震數據體的圖像可以如圖3-11所示,其中,該圖3-11是
對圖3-9所示的地震剖面處理后的圖,在該圖3-11中,橫軸表示垂直距離LS,單位為km,縱軸
表示時間深度τ,單位為s。圖3-12是圖3-11中的區域3110的放大圖。對比圖3-9與圖3-11,以
及圖3-10與圖3-12可見,應用本發明后,地震數據的分辨率明顯提高,原先疊合的同相軸
380被更好的分離開來,同相軸380橫向連續性變好。

圖3-13是處理前的疊后地震數據體在1200ms時刻的水平時間切片,圖3-14是應用
本發明實施例提供的疊后地震數據體處理方法對圖3-13所示的水平時間切片處理后的時
間切片圖,參見圖3-13和圖3-14,其中,橫坐標Line表示地震主測線的序號,縱坐標表示共
深度點(英文:Common Depth Point;簡稱:CDP),通過對比圖3-13和圖3-14可以看出,應用
本發明實施例提供的疊后地震數據體處理方法對圖3-13所示的疊后地震數據體在1200ms
時刻的水平時間切片處理后,河道空間展布更加清楚,細節信息明顯增多。

綜上所述,本發明實施例提供的疊后地震數據體處理方法,通過從地震工區的疊
后地震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣
本測線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣
本測線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條
樣本測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測
線對應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值
場;根據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法
更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更
新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數
據體,以便于GPU根據地震工區的Q值場對疊后地震數據體進行反Q濾波處理,得到處理后的
疊后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決了對
疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體處理
的靈活性和可靠性的效果。

進一步地,現有技術中,通常在CPU中對疊后地震數據體進行處理,由于CPU在一個
計算周期內只能對一個地震數據道進行處理,且其為單道輸入輸出,讀寫頻繁,處理效率較
低。本發明實施例提供的疊后地震數據體處理方法,通過CPU對疊后地震數據體進行分塊處
理得到地震數據塊,GPU根據地震數據塊對疊后地震數據體進行處理,由于GPU在一個計算
周期內可以對多個地震數據道進行處理,GPU的處理能力較強,本發明實施例通過采用GPU
對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效率。

本發明實施例提供的疊后地震數據體處理方法可以提高疊后地震數據的分辨率,
特別適用于處理大規模的三維偏移疊加地震數據體,對提高疊后地震數據體的分辨率具有
重要應用價值。

下述為本發明裝置實施例,可以用于執行本發明方法實施例。對于本發明裝置實
施例中未披露的細節,請參照本發明方法實施例。

請參考圖4,其示出了本發明實施例提供的一種疊后地震數據體處理裝置400的框
圖,該疊后地震數據體處理裝置400可以用于執行圖1所示實施例提供的疊后地震數據體處
理方法或者圖3-1所示實施例提供的疊后地震數據體處理方法中的部分方法,參見圖4,該
疊后地震數據體處理裝置400可以包括:

選取模塊401,用于從地震工區的疊后地震數據體對應的地震測線中選取M條樣本
測線,M為大于或者等于2的整數,且M=0.003N,N為疊后地震數據體的地震主測線的條數,
每條地震測線對應一個地震剖面;

第一確定模塊402,用于從地震工區的疊加速度場中確定M條樣本測線中的每條樣
本測線對應的速度場,得到M條樣本測線對應的樣本速度場ν;

第二確定模塊403,用于根據M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公
式確定M條樣本測線對應的初始Q值場Q0,李慶忠經驗公式為Q=α·νβ,α=14,β=2.2;

校正模塊404,用于對M條樣本測線對應的初始Q值場Q0進行校正,得到M條樣本測
線對應的校正Q值場Q1;

第三確定模塊405,用于根據M條樣本測線對應的校正Q值場Q1和M條樣本測線對應
的加權系數場λfinal,確定M條樣本測線對應的最終Q值場Q2;

更新模塊406,用于根據M條樣本測線對應的樣本速度場ν和M條樣本測線對應的最
終Q值場Q2,采用最小二乘法更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;

第四確定模塊407,用于根據地震工區的疊加速度場,采用更新后的李慶忠經驗公
式確定地震工區的Q值場Qall;

發送模塊408,用于向圖形處理器GPU發送地震工區的Q值場Qall和疊后地震數據
體,以便于GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理后
的疊后地震數據體。

綜上所述,本發明實施例提供的疊后地震數據體處理裝置,通過從地震工區的疊
后地震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣
本測線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣
本測線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條
樣本測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測
線對應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值
場;根據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法
更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更
新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數
據體,以便于GPU根據地震工區的Q值場對疊后地震數據體進行反Q濾波處理,得到處理后的
疊后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決了對
疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體處理
的靈活性和可靠性的效果。

請參考圖5,其示出了本發明實施例提供的另一種疊后地震數據體處理裝置400的
框圖,該疊后地震數據體處理裝置500可以用于執行圖1所示實施例提供的疊后地震數據體
處理方法或者圖3-1所示實施例提供的疊后地震數據體處理方法中的部分方法,參見圖5,
該疊后地震數據體處理裝置500可以包括但不限于:

選取模塊501,用于從地震工區的疊后地震數據體對應的地震測線中選取M條樣本
測線,M為大于或者等于2的整數,且M=0.003N,N為疊后地震數據體的地震主測線的條數,
每條地震測線對應一個地震剖面;

第一確定模塊502,用于從地震工區的疊加速度場中確定M條樣本測線中的每條樣
本測線對應的速度場,得到M條樣本測線對應的樣本速度場ν;

第二確定模塊503,用于根據M條樣本測線對應的樣本速度場ν,采用李慶忠經驗公
式確定M條樣本測線對應的初始Q值場Q0,李慶忠經驗公式為Q=α·νβ,α=14,β=2.2;

校正模塊504,用于對M條樣本測線對應的初始Q值場Q0進行校正,得到M條樣本測
線對應的校正Q值場Q1;

第三確定模塊505,用于根據M條樣本測線對應的校正Q值場Q1和M條樣本測線對應
的加權系數場λfinal,確定M條樣本測線對應的最終Q值場Q2;

更新模塊506,用于根據M條樣本測線對應的樣本速度場ν和M條樣本測線對應的最
終Q值場Q2,采用最小二乘法更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;

第四確定模塊507,用于根據地震工區的疊加速度場,采用更新后的李慶忠經驗公
式確定地震工區的Q值場Qall;

發送模塊508,用于向圖形處理器GPU發送地震工區的Q值場Qall和疊后地震數據
體,以便于GPU根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處理,得到處理后
的疊后地震數據體。

進一步地,請繼續參考圖5,該疊后地震數據體處理裝置500還可以包括:

第一處理模塊509,用于根據預設的加權系數組對M條樣本測線對應的校正Q值場
Q1進行處理,得到M條樣本測線對應的加權Q值場組,預設的加權系數組中包括P個加權系數
λ,加權Q值場組中包括P套加權Q值場,P為大于或者等于1的整數;

第二處理模塊510,用于根據加權Q值場組對M條樣本測線組成的樣本地震數據體
進行反Q濾波處理,得到反Q濾波數據體組,反Q濾波數據體組中包括P套反Q濾波數據體,每
套反Q濾波數據體包括M條反Q樣本測線,每條反Q樣本測線對應一個地震剖面;

頻譜分析模塊511,用于采用傅里葉變換對反Q濾波數據體組中的每套反Q濾波數
據體中的每條反Q樣本測線對應的地震剖面進行頻譜分析;

第五確定模塊512,用于根據頻譜分析結果進行處理,確定M條反Q樣本測線對應的
加權系數場;

第六確定模塊513,用于將M條反Q樣本測線對應的加權系數場確定為M條樣本測線
對應的加權系數場λfinal。

進一步地,頻譜分析模塊511,具體用于:

在反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣本測線對應的地震剖
面上選擇三個矩形時窗,三個矩形時窗包括:淺層矩形時窗、中層矩形時窗和深層矩形時
窗;

采用傅里葉變換對每條反Q樣本測線對應的地震剖面上位于三個矩形時窗中的每
個矩形時窗內的部分進行處理,得到每條反Q樣本測線對應的3個反Q樣本頻譜;

根據每套反Q濾波數據體中的M條反Q樣本測線對應的反Q樣本頻譜,得到每套反Q
濾波數據體對應的3M個反Q樣本頻譜;

根據反Q濾波數據體組中的P套反Q濾波數據體對應的反Q樣本頻譜得到3×M×P個
反Q樣本頻譜,3×M×P個反Q樣本頻譜中包括:M×P個淺層矩形時窗對應的反Q樣本頻譜、M
×P個中層矩形時窗對應的反Q樣本頻譜和M×P個深層矩形時窗對應的反Q樣本頻譜;

對3×M×P個反Q樣本頻譜中的每個反Q樣本頻譜進行頻譜分析。

進一步地,第五確定模塊512,具體用于:

根據頻譜分析結果確定反Q濾波數據體組中的每套反Q濾波數據體中的每條反Q樣
本測線對應的第一目標頻譜、第二目標頻譜和第三目標頻譜,第一目標頻譜為每條反Q樣本
測線對應的地震剖面上位于淺層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成
的封閉區域的面積最大的反Q樣本頻譜,第二目標頻譜為每條反Q樣本測線對應的地震剖面
上位于中層矩形時窗內的部分對應的P個反Q樣本頻譜中與橫軸圍成的封閉區域的面積最
大的反Q樣本頻譜,第三目標頻譜為每條反Q樣本測線對應的地震剖面上位于深層矩形時窗
內的部分對應的P個反Q樣本頻譜中與橫軸圍成的封閉區域的面積最大的反Q樣本頻譜;

分別確定第一目標頻譜對應的矩形時窗對應的第一目標加權系數、第二目標頻譜
對應的矩形時窗對應的第二目標加權系數和第三目標頻譜對應的矩形時窗對應的第三目
標加權系數,得到每條反Q樣本測線對應的3個目標加權系數,其中,第一目標加權系數、第
二目標加權系數和第三目標加權系數都為預設的加權系數組中的加權系數;

根據M條反Q樣本測線中的每條反Q樣本測線對應的3個目標加權系數,得到M條反Q
樣本測線對應的3M個目標加權系數;

對3M個目標加權系數進行空間插值,得到地震工區對應的加權系數場;

對地震工區對應的加權系數場進行平滑處理,得到處理后的加權系數場;

從處理后的加權系數場中抽取出M條反Q樣本測線對應的加權系數場。

進一步地,請繼續參考圖5,該疊后地震數據體處理裝置500還可以包括:

采樣模塊514,用于按照預設采樣間隔對地震工區的Q值場Qall進行采樣,得到采樣
Q值場Qs;

分塊模塊515,用于對疊后地震數據體進行分塊處理,得到D個地震數據塊,D為大
于或者等于2的整數;

發送模塊508,具體用于:

依次向GPU發送采樣Q值場Qs和D個地震數據塊中的每個地震數據塊。

進一步地,分塊模塊515,具體用于:

計算疊后地震數據體中的地震數據的總道數E;

根據GPU在一個周期內處理的地震數據的道數Tracemax和疊后地震數據體中的地
震數據的總道數E,對疊后地震數據體進行分塊處理,得到D個地震數據塊;

其中,E=Tracemax×(D-1)+L,L為D個地震數據塊中的最后一個地震數據塊中的地
震數據的道數。

綜上所述,本發明實施例提供的疊后地震數據體處理裝置,通過從地震工區的疊
后地震數據體對應的地震測線中選取M條樣本測線;從地震工區的疊加速度場中確定M條樣
本測線中的每條樣本測線對應的速度場,得到M條樣本測線對應的樣本速度場;根據M條樣
本測線對應的樣本速度場,采用李慶忠經驗公式確定M條樣本測線對應的初始Q值場;對M條
樣本測線對應的初始Q值場進行校正,得到M條樣本測線對應的校正Q值場;根據M條樣本測
線對應的校正Q值場和M條樣本測線對應的加權系數場,確定M條樣本測線對應的最終Q值
場;根據M條樣本測線對應的樣本速度場和M條樣本測線對應的最終Q值場,采用最小二乘法
更新李慶忠經驗公式,得到更新后的李慶忠經驗公式;根據地震工區的疊加速度場,采用更
新后的李慶忠經驗公式確定地震工區的Q值場;向GPU發送地震工區的Q值場和疊后地震數
據體,以便于GPU根據地震工區的Q值場對疊后地震數據體進行反Q濾波處理,得到處理后的
疊后地震數據體。由于針對疊后地震數據體對李慶忠經驗公式進行了更新,因此,解決了對
疊后地震數據體處理的靈活性較低、可靠性較差的問題,達到了提高疊后地震數據體處理
的靈活性和可靠性的效果。

進一步地,現有技術中,通常在CPU中對疊后地震數據體進行處理,由于CPU在一個
計算周期內只能對一個地震數據道進行處理,且其為單道輸入輸出,讀寫頻繁,處理效率較
低。本發明實施例提供的疊后地震數據體處理裝置,通過CPU對疊后地震數據體進行分塊處
理得到地震數據塊,GPU根據地震數據塊對疊后地震數據體進行處理,由于GPU在一個計算
周期內可以對多個地震數據道進行處理,GPU的處理能力較強,本發明實施例通過采用GPU
對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效率。

本發明實施例提供的疊后地震數據體處理裝置可以提高疊后地震數據的分辨率,
特別適用于處理大規模的三維偏移疊加地震數據體,對提高疊后地震數據體的分辨率具有
重要應用價值。

請參考圖6,其示出了本發明實施例提供的再一種疊后地震數據體處理裝置400的
框圖,該疊后地震數據體處理裝置600可以用于執行圖2所示實施例提供的疊后地震數據體
處理方法或者圖3-1所示實施例提供的疊后地震數據體處理方法中的部分方法,參見圖6,
該疊后地震數據體處理裝置600可以包括:

接收模塊610,用于接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數
據體;

處理模塊620,用于根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處
理,得到處理后的疊后地震數據體。

綜上所述,本發明實施例提供的疊后地震數據體處理裝置,通過接收CPU發送的地
震工區的Q值場和疊后地震數據體,根據地震工區的Q值場對疊后地震數據體進行反Q濾波
處理,得到處理后的疊后地震數據體,由于地震工區的Q值場是CPU根據更新后的李慶忠經
驗公式計算得到的,且更新后的李慶忠經驗公式是CPU針對疊后地震數據體對李慶忠經驗
公式進行更新得到的,因此,解決了對疊后地震數據體處理的靈活性較低、可靠性較差的問
題,達到了提高疊后地震數據體處理的靈活性和可靠性的效果。

請參考圖7-1,其示出了本發明實施例提供的又一種疊后地震數據體處理裝置400
的框圖,該疊后地震數據體處理裝置700可以用于執行圖2所示實施例提供的疊后地震數據
體處理方法或者圖3-1所示實施例提供的疊后地震數據體處理方法中的部分方法,參見圖
7-1,該疊后地震數據體處理裝置700可以包括:

接收模塊710,用于接收中央處理器CPU發送的地震工區的Q值場Qall和疊后地震數
據體;

處理模塊720,用于根據地震工區的Q值場Qall對疊后地震數據體進行反Q濾波處
理,得到處理后的疊后地震數據體。

進一步地,接收模塊710,具體用于:

接收CPU發送的采樣Q值場Qs和D個地震數據塊中的每個地震數據塊,采樣Q值場Qs
是CPU按照預設采樣間隔對地震工區的Q值場Qall進行采樣得到的,D個地震數據塊是CPU對
疊后地震數據體進行分塊處理得到的;

請參考圖7-2,其示出的是圖7-1所示實施例提供的一種處理模塊720的框圖,參見
圖7-2,該處理模塊720可以包括:

處理子模塊721,用于根據采樣Q值場Qs對D個地震數據塊中的每個地震數據塊進
行反Q濾波處理,得到D個處理后的地震數據塊;

得到子模塊722,用于根據D個處理后的地震數據塊,得到處理后的疊后地震數據
體。

進一步地,請參考圖7-3,其示出的是圖7-1所示實施例提供的一種處理子模塊721
的框圖,參見圖7-3,該處理子模塊721可以包括:

第一處理單元7211,用于對采樣Q值場Qs進行處理得到每個時間深度τ處的Q值Qτ;

第二處理單元7212,用于根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每
個地震數據塊中的所有時間深度τ處的地震幅值進行反Q濾波處理,得到D個處理后的地震
數據塊。

進一步地,第一處理單元7211,具體用于:

根據采樣Q值場Qs和Q值插值計算公式確定每個時間深度τ處的Q值Qτ;

其中,Q值插值計算公式為:

<mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>

<mrow> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> <msub> <mi>&tau;</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>int</mi> <mo>(</mo> <mfrac> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mn>0.02</mn> </mfrac> <mo>)</mo> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <mn>0.02</mn> <mo>,</mo> </mrow>τj、τ(j-1)和τ(j+1)都為時間深
度,為時間深度τj處的Q值插值,為時間深度τ(j-1)處的Q值插值,為時間深度
τ(j+1)處的Q值插值,表示對取整。

進一步地,請參考圖7-4,其示出的是圖7-1所示實施例提供的一種第二處理單元
7212的框圖,參見圖7-4,該第二處理單元7212可以包括:

處理子單元72121,用于根據所有時間深度τ處的Q值Qτ對D個地震數據塊中的每個
地震數據塊中的每個地震數據道的所有時間深度τ處的地震幅值進行反Q濾波處理,得到處
理后的地震數據道;

第一得到子單元72122,用于根據每個地震數據塊中的所有處理后的地震數據道,
得到處理后的地震數據塊;

第二得到子單元72123,用于根據D個地震數據塊中的所有處理后的地震數據塊,
得到D個處理后的地震數據塊。

進一步地,請參考圖7-5,其示出的是圖7-1所示實施例提供的一種處理子單元
72121的框圖,參見圖7-5,該處理子單元72121可以包括:

第一確定子單元721211,用于分別確定GPU中的線程的線程格尺寸和線程塊尺寸;

變換子單元721212,用于對D個地震數據塊中的每個地震數據塊中的每個地震數
據道進行快速傅里葉變換;

第二確定子單元721213,用于根據快速傅里葉變換結果,確定D個地震數據塊中所
有地震數據對應的角頻率采樣間隔△ω;

第三確定子單元721214,用于根據角頻率采樣間隔△ω和預先接收的目標低角頻
率值ωmin、目標高角頻率值ωmax,確定每個地震數據道對應的起算角頻率值ωp和終止角頻
率值ωz;

第一計算子單元721215,用于根據GPU中的線程的線程格尺寸、線程塊尺寸、所有
時間深度τ處的Q值Qτ中的時間深度τj處的Q值每個地震數據道對應的起算角頻率值ωp
和終止角頻率值ωz,計算每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅
值A;

第二計算子單元721216,用于根據每個地震數據道對應的時間深度τj處的每個頻
率點對應的地震幅值A,計算每個地震數據道對應的時間深度τj處的所有頻率點對應的地
震幅值∑A;

第一處理子單元72127,用于根據GPU中的所有線程處理得到的每個地震數據道對
應的所有時間深度處的所有頻率點對應的地震幅值,處理得到處理后的地震數據道。

進一步地,GPU中的線程的線程格尺寸為:Tracemax為
GPU在一個周期內處理的地震數據的道數,b_size=2m,m為正整數,且4≤m≤6,S為地震數
據道的采樣點數;

GPU中的線程的線程塊尺寸為:(b_size,b_size),b_size=2m,m為正整數,且4≤m
≤6;

第二確定子單元721213,具體用于:

根據快速傅里葉變換結果,采用角頻率采樣間隔計算公式確定D個地震數據塊中
所有地震數據對應的角頻率采樣間隔△ω;

其中,角頻率采樣間隔計算公式為:△ω=2π/(Gdt),G為滿足快速傅里葉變換算
法的時間采樣間隔,dt為地震數據的時間采樣間隔;

第三確定子單元721214,具體用于:

根據目標低角頻率值ωmin和角頻率采樣間隔△ω,采用起算角頻率計算公式計算
每個地震數據道對應的起算角頻率值ωp;

根據目標高角頻率值ωmax和角頻率采樣間隔△ω,采用終止角頻率計算公式計算
每個地震數據道對應的終止角頻率值ωz;

其中,起算角頻率計算公式為:終止角頻率計算公式為:
△f為頻率增量,ωmax+2π△f為通過反Q濾波處理后期望達到的角頻
率的上限;

第一計算子單元721215,具體用于:

根據GPU中的線程的線程格尺寸、線程塊尺寸、所有時間深度τ處的Q值Qτ中的時間
深度τj處的Q值每個地震數據道對應的起算角頻率值ωp和終止角頻率值ωz,采用地震
幅值計算公式計算每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A;

其中,地震幅值計算公式為:
A為每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值且A為復數,k為每
個頻率點的序號,△ω為角頻率采樣間隔,τj為時間深度,為時間深度τj處的Q值;

第二計算子單元721216,具體用于:

根據每個地震數據道對應的時間深度τj處的每個頻率點對應的地震幅值A,采用
幅值和計算公式計算每個地震數據道對應的時間深度τj處的所有頻率點對應的地震幅值
∑A;

其中,幅值和計算公式為:

<mrow> <mo>&Sigma;</mo> <mi>A</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&omega;</mi> <mi>min</mi> </msub> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>int</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&omega;</mi> <mi>max</mi> </msub> <mo>+</mo> <mn>2</mn> <mi>&pi;</mi> <mi>&Delta;</mi> <mi>f</mi> </mrow> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </munderover> <mi>Re</mi> <mo>{</mo> <mi>F</mi> <mrow> <mo>(</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>k&Delta;&omega;&tau;</mi> <mi>j</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>Q</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </msub> </mrow> </mfrac> <mo>)</mo> <mo>}</mo> <mo>,</mo> </mrow>

A為每個地震數據道對應的時間深度τj處的每個頻
率點對應的地震幅值且A為復數,k為每個頻率點的序號,△ω為角頻率采樣間隔,τj為時間
深度,為時間深度τj處的Q值,Re表示取實部。

綜上所述,本發明實施例提供的疊后地震數據體處理裝置,通過接收CPU發送的地
震工區的Q值場和疊后地震數據體,根據地震工區的Q值場對疊后地震數據體進行反Q濾波
處理,得到處理后的疊后地震數據體,由于地震工區的Q值場是CPU根據更新后的李慶忠經
驗公式計算得到的,且更新后的李慶忠經驗公式是CPU針對疊后地震數據體對李慶忠經驗
公式進行更新得到的,因此,解決了對疊后地震數據體處理的靈活性較低、可靠性較差的問
題,達到了提高疊后地震數據體處理的靈活性和可靠性的效果。

進一步地,現有技術中,通常在CPU中對疊后地震數據體進行處理,由于CPU在一個
計算周期內只能對一個地震數據道進行處理,且其為單道輸入輸出,讀寫頻繁,處理效率較
低。本發明實施例提供的疊后地震數據體處理裝置,通過CPU對疊后地震數據體進行分塊處
理得到地震數據塊,GPU根據地震數據塊對疊后地震數據體進行處理,由于GPU在一個計算
周期內可以對多個地震數據道進行處理,GPU的處理能力較強,本發明實施例通過采用GPU
對疊后地震數據體進行反Q濾波處理,提高了對疊后地震數據體的處理速度和處理效率。

本發明實施例提供的疊后地震數據體處理裝置可以提高疊后地震數據的分辨率,
特別適用于處理大規模的三維偏移疊加地震數據體,對提高疊后地震數據體的分辨率具有
重要應用價值。

需要說明的是:上述實施例提供的疊后地震數據體處理裝置在處理疊后地震數據
體時,僅以上述各功能模塊的劃分進行舉例說明,實際應用中,可以根據需要而將上述功能
分配由不同的功能模塊完成,即將設備的內部結構劃分成不同的功能模塊,以完成以上描
述的全部或者部分功能。另外,上述實施例提供的疊后地震數據體處理裝置與疊后地震數
據體處理方法實施例屬于同一構思,其具體實現過程詳見方法實施例,這里不再贅述。

本領域普通技術人員可以理解實現上述實施例的全部或部分步驟可以通過硬件
來完成,也可以通過程序來指令相關的硬件完成,所述的程序可以存儲于一種計算機可讀
存儲介質中,上述提到的存儲介質可以是只讀存儲器,磁盤或光盤等。

以上所述僅為本發明的較佳實施例,并不用以限制本發明,凡在本發明的精神和
原則之內,所作的任何修改、等同替換、改進等,均應包含在本發明的保護范圍之內。

關 鍵 詞:
地震 數據 處理 方法 裝置
  專利查詢網所有資源均是用戶自行上傳分享,僅供網友學習交流,未經上傳用戶書面授權,請勿作他用。
關于本文
本文標題:疊后地震數據體處理方法及裝置.pdf
鏈接地址:http://www.rgyfuv.icu/p-6100667.html
關于我們 - 網站聲明 - 網站地圖 - 資源地圖 - 友情鏈接 - 網站客服客服 - 聯系我們

[email protected] 2017-2018 zhuanlichaxun.net網站版權所有
經營許可證編號:粵ICP備17046363號-1 
 


收起
展開
山东11选5中奖结果走势图