[發明專利]一種基于預計算插值的射線走時計算方法及系統

有效

• 申請號: 201410363926.7
• 申請日: 2014-07-28
• 公開（公告）號: CN104133238B
• 公開（公告）日: 2016-10-26
• 發明（設計）人: 梁兼棟;郝曉光
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 湯在彥

技術領域

本發明是關于石油勘探技術，主要涉及地震數據處理技術，尤指一種基于預計算插值的射線走時計算方法。

背景技術

在地震數據處理技術中，介質速度的計算可以利用相干反演技術，根據給定的若干個掃描速度，對其中每一個掃描速度生成一個相應的介質模型，在此模型上計算射線走時并與道集數據比較得出在當前速度模型下的相似度，最終獲得介質速度。因此，射線走時的計算很大程度影響到最終結果的精確程度和計算效率。現有的射線走時計算方法一般是采用射線追蹤方法，這種方法具有很高的精度，但是計算效率低下，而且在大規模數據處理中往往存在嚴重的重復計算問題。

發明內容

本發明為了克服現有的射線走時計算方法的不足，提出了一種基于預計算插值的射線走時計算方法及系統，以達到高效率的走時計算并能保證計算精度。

為達到上述目的，本發明提出的基于預計算插值的射線走時計算方法包括以下步驟：步驟1，采集實驗地區數據建立所述實驗地區的速度模型，在所述速度模型上設定炮點、檢波點，并設定所述速度模型的目標層位；步驟2，在每一所述炮點或檢波點的位置發射一組射線，射線到達所述目標層位后終止，并記錄每一所述炮點或檢波點發出的一組射線的射線數據；步驟3，根據所述炮點或檢波點發出的一組射線的射線數據對所述目標層位的CMP點集進行組內插值計算，獲得所述目標層位的CMP點集對應所述炮點或檢波點的射線數據；步驟4，在所述炮點或檢波點中選取第一炮點及第一檢波點，并從所述目標層位的CMP點集對應所述炮點或檢波點的射線數據中獲得CMP點集對應第一炮點及第一檢波點的射線數據；步驟5，根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過迭代插值計算獲得所述第一炮點及第一檢波點在目標層位上的反射點；步驟6，根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過插值計算獲得所述反射點到述第一炮點及第一檢波點的射線走時之和。

為達到上述目的，本發明還提出了一種基于預計算插值的射線走時計算系統，包括：模型建立模塊，用于根據實驗地區的數據建立所述實驗地區的速度模型，在所述速度模型上設定炮點、檢波點，并設定所述速度模型的目標層位；射線數據采集模塊，用于在每一所述炮點或檢波點的位置發射一組射線，射線到達所述目標層位后終止，并記錄每一所述炮點或檢波點發出的一組射線的射線數據；插值計算模塊，用于根據所述炮點或檢波點發出的一組射線的射線數據對所述目標層位的CMP點集進行組內插值計算，獲得所述目標層位的CMP點集對應所述炮點或檢波點的射線數據；射線數據選取模塊，用于在所述炮點或檢波點中選取第一炮點及第一檢波點，并從所述目標層位的CMP點集對應所述炮點或檢波點的射線數據中獲得CMP點集對應第一炮點及第一檢波點的射線數據；反射點獲取模塊，用于根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過迭代插值計算獲得所述第一炮點及第一檢波點在目標層位上的反射點；射線走時計算模塊，用于根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過插值計算獲得所述反射點到述第一炮點及第一檢波點的射線走時之和。

本發明的基于預計算插值的射線走時計算方法及系統能夠高效地計算相干反演算法中所需要的射線走時信息，并且也能達到較高的精度，同時該算法可以并行運行，提高了地震勘探數據處理的效率。

附圖說明

此處所說明的附圖用來提供對本發明的進一步理解，構成本申請的一部分，并不構成對本發明的限定。在附圖中：

圖1為本發明一實施例的基于預計算插值的射線走時計算方法流程圖。

圖2為本發明一實施例的基于預計算插值的射線走時計算方法中步驟3的詳細流程圖。

圖3為本發明一實施例的基于預計算插值的射線走時計算方法中步驟5的詳細流程圖。

圖4為本發明一具體實施例的射線方向與目標層位的法向量的夾角示意圖。

圖5為本發明一實施例的基于預計算插值的射線走時計算系統的結構示意圖。

圖6為本發明一實施例的基于預計算插值的射線走時計算系統中插值計算模塊的結構示意圖。

圖7為本發明一實施例的基于預計算插值的射線走時計算系統中反射點獲取模塊的結構示意圖。

具體實施方式

以下配合圖式及本發明的較佳實施例，進一步闡述本發明為達成預定發明目的所采取的技術手段。

圖1為本發明一實施例的基于預計算插值的射線走時計算方法流程圖。如圖1所示，該方法包括以下步驟：

步驟1，采集實驗地區數據建立實驗地區的速度模型，在速度模型上設定炮點、檢波點，并設定速度模型的目標層位。