[發明專利]拾取地震波時距曲線的方法及裝置

有效

• 申請號: 201910217070.5
• 申請日: 2019-03-21
• 公開（公告）號: CN111722274B
• 公開（公告）日: 2023-02-07
• 發明（設計）人: 張文波;李合群;李建峰;聶紅梅;谷逸薇;馬淵明
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 王天堯;任默聞

1.一種拾取地震波時距曲線的方法，其特征在于，所述方法包括：

設定包含時距曲線的時窗與時窗內的時間樣點數，獲取每個時窗內的地震記錄；

將獲取的地震記錄中信噪比最高的地震記錄對應的地震道確定為參考地震道；

確定目標地震道，將每個時窗內獲取的目標地震道的地震記錄與相同時窗內獲取的參考地震道的地震記錄進行同相疊加，得到疊加結果，所述目標地震道為拾取時距曲線的地震道；

根據所述疊加結果，確定目標地震道的地震波時距曲線相對于參考地震道的地震波時距曲線的時移量；

根據參考地震道的地震波時距曲線以及目標地震道的地震波時距曲線相對于參考地震道的地震波時距曲線的時移量，確定目標地震道的地震波時距曲線。

2.根據權利要求1所述的方法，其特征在于，所述將采集的地震記錄中信噪比最高的地震記錄對應的地震道確定為參考地震道，包括：

確定每一道地震記錄分別與所有道地震記錄之間的互相關函數；

構造地震記錄互相關函數矩陣，所述互相關函數矩陣中每一列為一道地震記錄與地震記錄中的每一道地震記錄之間的互相關函數；

確定互相關函數矩陣中每一列互相關函數的平均相似系數；

將所述平均相似系數最高的一列確定為參考列，將所述參考列對應的地震道確定為參考地震道。

3.根據權利要求2所述的方法，其特征在于，所述確定每一道地震記錄分別與所有道地震記錄之間的互相關函數，包括：

利用公式計算每一道地震記錄分別與所有道地震記錄之間的互相關函數；

其中，k和l為地震道的序號，t為時間樣點，Φ_k,l(τ)為第k道與第l道相關的互相關函數，S_k(t)為第k道地震道在時間樣點t采集的地震記錄，S_l(t)為第l道地震道在時間樣點t采集的地震記錄，τ為S_l(t)的延遲時間，k＝1,2,…,n，l＝1,2,…,n，n為設定時窗內的地震道數，t＝1,2,…,m，m為時間樣點的數量。

4.根據權利要求3所述的方法，其特征在于，確定互相關函數矩陣中每一列互相關函數的平均相似系數，包括：

利用公式計算互相關函數矩陣中每一列互相關函數的平均相似系數；

其中，S_k為第k道的平均相似系數，τ₀為使S_k達到最大值時的延遲時間。

5.根據權利要求3或4所述的方法，其特征在于，在將所述平均相似系數最高的一列確定為參考列之后，所述方法還包括：

按照每列對應的平均相似系數由高到低的順序，對互相關函數矩陣中的列進行重新排序，得到排序后的互相關函數矩陣，所述排序后的互相關函數矩陣中的第1列為參考列。

6.根據權利要求5所述的方法，其特征在于，所述將每個時窗內采集的目標地震道的地震記錄與相同時窗內采集的參考地震道的地震記錄進行同相疊加，得到疊加結果，包括：

利用公式計算每列對應的二階互相關函數；

其中，ΦΦ_1,l(δ)為排序后的互相關函數矩陣中，第l道對應的互相關函數與第1列對應行上的互相關函數的二階互相關函數，δ為第l道二階互相關函數的延遲時間，Φ'_1,k(τ)為排序后的互相關函數矩陣中，第k道與第1列對應的地震道相關的互相關函數，Φ'_l,k(τ+δ)為排序后的互相關函數矩陣中，第k道與第l道相關的互相關函數。

7.根據權利要求6所述的方法，其特征在于，根據所述疊加結果，確定目標地震道的時距曲線相對于參考地震道的時距曲線的時移量，包括：

構造函數F(δ)＝ΦΦ_1,l(δ)，確定使函數F(δ)取最大值時δ的值δ₀；

將δ₀確定為第l道地震道的時距曲線相對于參考地震道的時距曲線的時移量。