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

有效

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

1.一種基于預計算插值的射線走時計算方法，其特征在于，包括以下步驟：

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

步驟2，在每一所述炮點或檢波點的位置發射一組射線，射線到達所述目標層位后終止，并記錄每一所述炮點或檢波點發出的一組射線的射線數據；

步驟3，根據所述炮點或檢波點發出的一組射線的射線數據對所述目標層位的CMP點集進行組內插值計算，獲得所述目標層位的CMP點集對應所述炮點或檢波點的射線數據；

步驟4，在所述炮點或檢波點中選取第一炮點及第一檢波點，并從所述目標層位的CMP點集對應所述炮點或檢波點的射線數據中獲得CMP點集對應第一炮點及第一檢波點的射線數據；

步驟5，根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過迭代插值計算獲得所述第一炮點及第一檢波點在目標層位上的反射點；

步驟6，根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過插值計算獲得所述反射點到述第一炮點及第一檢波點的射線走時之和。

2.根據權利要求1所述的基于預計算插值的射線走時計算方法，其特征在于，所述步驟2中，每一所述炮點或檢波點發出的一組射線的射線數據包括：

每一所述炮點或檢波點發出的一組射線中的每一射線與所述目標層位的交點位置，所述每一射線在交點位置的射線方向，所述每一射線到交點位置的走時。

3.根據權利要求2所述的基于預計算插值的射線走時計算方法，其特征在于，所述步驟3，根據所述炮點或檢波點發出的一組射線的射線數據對所述目標層位的CMP點集進行組內插值計算，獲得所述目標層位的CMP點集對應所述炮點或檢波點的射線數據包括：

步驟31，根據所述每一射線與所述目標層位的交點位置，記錄所述炮點或檢波點發出的一組射線中的相鄰兩條射線與所述目標層位的交點位置之間覆蓋的CMP點；

步驟32，根據所述CMP點的位置以及所述相鄰兩條射線與所述目標層位的交點位置，將所述相鄰兩條射線在交點位置的射線方向以及所述相鄰兩條射線到交點位置的走時插值到所述CMP點，獲得所述CMP點對應所述炮點或檢波點的射線數據。

4.根據權利要求3所述的基于預計算插值的射線走時計算方法，其特征在于，所述步驟5，根據所述CMP點集對應第一炮點及第一檢波點的射線數據，通過迭代插值計算獲得所述第一炮點及第一檢波點在目標層位上的反射點包括：

步驟51，在第一炮點及第一檢波點之間的目標層位上設定一假設點；

步驟52，從所述CMP點集對應第一炮點及第一檢波點的射線數據中，獲取所述假設點兩側相鄰的兩個CMP點對應第一炮點及第一檢波點的射線數據，并通過插值獲得所述假設點對應第一炮點及第一檢波點的射線數據；

步驟53，根據所述假設點對應第一炮點及第一檢波點的射線數據，分別得到第一炮點的射線在所述假設點的射線方向與所述假設點的法向量的夾角、第一檢波點的射線在所述假設點的射線方向與所述假設點的法向量的夾角；

步驟54，比較所述步驟53得到的兩夾角的大小關系，并根據所述大小關系在所述假設點的相應一側以一定步長選取另一假設點；

重復執行所述步驟52至步驟54，當選取到某一假設點，使所述某一假設點的兩夾角之差小于一設定閾值時計算終止，所述某一假設點為第一炮點及第一檢波點在目標層位上的反射點。