[發明專利]地震數據處理方法和裝置

有效

• 申請號: 201610204476.6
• 申請日: 2016-04-01
• 公開（公告）號: CN106226818B
• 公開（公告）日: 2018-05-04
• 發明（設計）人: 張光榮;楊躍明;肖富森;張世榮;喻頤;冉崎;廖奇;周肖;馬波;謝冰;梁翰;賴強;陳康;張旋
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/36
• 分類號: G01V1/36;G01V1/30
• 代理公司: 北京三高永信知識產權代理有限責任公司11138
• 代理人: 周靜

技術領域

本發明涉及地質勘測領域，特別涉及一種地震數據處理方法和裝置。

背景技術

地震勘探技術是利用人工激發的地震波在地層內傳播來勘探地下的地質情況。在地面某處激發的地震波向地下傳播時，遇到不同彈性的地層分界面就會產生反射波或折射波返回地面，用專門的儀器可記錄這些波，分析所記錄的波的特點(如波的傳播時間、振動形狀等)，能較準確地測定這些界面的深度和形態，判斷地層的巖性。

由于地震記錄是有噪音的，因此通常會通過不同位置進行激發、接收地下同一點的地震反射波，將這些來自同一點的地震反射波的集合稱為地震道集。通過多次地震信號疊加，可以達到消除噪音，提高信噪比的目的。

由于野外采集的地震資料統一反射點的地震反射時間是不一樣的，其反射時間與偏移距、地震傳播速度有關，因此，需要將不同偏移距地震記錄時間校正到自激自收地震反射時間位置，這個過程叫動校正，其目的是為了便于對地震信號疊加。

但是，因為地震波的傳播速度是未知的，因此動校正的實現過程通常都設計一系列地震傳播速度對地震道集進行動校正處理。

由于野外地震采集數據量非常龐大，人們不可能對每一個地震反射點都進行速度分析。往往都是間隔幾十道才對一個地震反射點進行速度分析，而中間未進行速度分析的道集采用已作速度分析道集的速度進行插值而得，由于地震速度的非均質性，這些插值的速度或多或少都會存在一定的誤差，動校正結果往往很難達到理想效果，很難將所有地震道集同相軸都校正到水平狀態。

發明內容

為了解決現有技術的問題，本發明實施例提供了一種地震數據處理方法和裝置。所述技術方案如下：

一方面，本發明實施例提供了一種地震數據處理方法，所述方法包括：

對動校正處理后的S個地震道集中預設范圍內的地震道集進行疊加得到模型道，S為整數；

計算各個地震道集與所述模型道的相關系數，并選出相關系數最大的第K地震道集；

在預設時間范圍內滑動時窗，計算第K-1地震道集與所述第K地震道集的相關系數，并選取相關系數最大時的時窗作為最優點，將最優點對應的時移量作為第K-1地震道集的剩余動校正量，K＞1且K為整數；

在所述預設時間范圍內滑動時窗，計算第M地震道集與第M+1地震道集的相關系數，并選取相關系數最大時的時窗作為最優點，將最優點對應的時移量作為第M地震道集的剩余動校正量，K-2≥M≥1且M為整數；

在所述預設時間范圍內滑動時窗，計算第K+1地震道集與所述第K地震道集的相關系數，并選取相關系數最大時的時窗作為最優點，將最優點對應的時移量作為第K+1地震道集的剩余動校正量；

在所述預設時間范圍內滑動時窗，計算第N+1地震道集與第N地震道集的相關系數，并選取相關系數最大時的時窗作為最優點，將最優點對應的時移量作為第N+1地震道集的剩余動校正量，S-1≥N≥K+1且N為整數，所述最優點對應的時移量為時窗內的一個地震道集的地震反射時間與偏移距為0的地震道集的地震反射時間的差值；

采用計算出的所述S個地震道集的剩余動校正量對所述S個地震道集進行剩余動校正。

在本發明實施例的一種實現方式中，所述對動校正處理后的S個地震道集中預設范圍內的地震道集進行疊加得到模型道，包括：

采用如下公式對多個地震道集進行疊加得到模型道：

其中，m1、m2為地震道集序號，a_ij為第i道的第j采樣點的振幅值，m1、m2、i和j均為正整數，1≤m1＜m2≤S，i∈[m1，m2]。

在本發明實施例的另一種實現方式中，所述計算各個地震道集與所述模型道的相關系數，包括：

采用以下公式計算模型道與地震道集的相關系數：

其中，為第x道與模型道的相關系數，a_j為模型道第j個采樣點的振幅值，Aa為模型道振幅平均值，b_j為地震道集第j個采樣點的振幅值，Ba為地震道集振幅平均值，n為地震道集中采樣點的總數，n、x和j均為正整數，且1≤x≤S。

在本發明實施例的另一種實現方式中，所述方法還包括：

在對動校正處理后的S個地震道集中預設范圍內的地震道集進行疊加得到模型道之前，對所述各個地震道集進行振幅能量均衡化處理。

在本發明實施例的另一種實現方式中，所述對所述各個地震道集進行振幅能量均衡化處理，包括：

在設定時窗范圍內統計振幅能量：

Q為統計道數，P為采樣點數，a_j為振幅值；