[發(fā)明專利]一種地震疊前時(shí)間偏移中的射線追蹤方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310739693.1
• 申請(qǐng)日: 2013-12-26
• 公開（公告）號(hào): CN103698813A
• 公開（公告）日: 2014-04-02
• 發(fā)明（設(shè)計(jì)）人: 劉洋
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣集團(tuán)公司;中國(guó)石油大學(xué)(北京)
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/36
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司 11127
• 代理人: 任默聞

1.一種地震疊前時(shí)間偏移中的射線追蹤方法，其特征在于，所述方法包括：

步驟1，獲取M個(gè)反射界面的模型參數(shù)以及N個(gè)炮檢距參數(shù)；

步驟2，根據(jù)K個(gè)界面的模型參數(shù)和炮檢距參數(shù)，計(jì)算來自K個(gè)界面第j個(gè)炮檢距的反射波傳播路徑初值參數(shù)a_i⁽¹⁾及傳播時(shí)間i=1,2,...,M，k=1,2,...,M，j=1,2,...,N；

步驟3，根據(jù)設(shè)定的迭代系數(shù)，利用迭代法計(jì)算來自K個(gè)界面第j個(gè)炮檢距的新的反射波傳播路徑初值參數(shù)a_i⁽²⁾，并計(jì)算來自K個(gè)界面第j個(gè)炮檢距的新的反射波傳播時(shí)間

步驟4，如果ε為設(shè)定的誤差，則即為來自K個(gè)界面第j個(gè)炮檢距的反射波傳播時(shí)間，否則，返回所述步驟3繼續(xù)迭代；

步驟5，如果j<N，則轉(zhuǎn)向步驟2-4繼續(xù)計(jì)算，否則，來自K個(gè)界面反射波傳播時(shí)間計(jì)算完成；

步驟6，如果K<M，則轉(zhuǎn)向步驟2-5繼續(xù)計(jì)算，否則，最終的反射波傳播時(shí)間計(jì)算完成。

2.根據(jù)權(quán)利要求1所述的地震疊前時(shí)間偏移中的射線追蹤方法，其特征在于，在所述步驟1中，獲取的M個(gè)反射界面的模型參數(shù)包括：各層層速度v_k以及各層厚度h_k，其中，k=1,2,...,M。

3.根據(jù)權(quán)利要求2所述的地震疊前時(shí)間偏移中的射線追蹤方法，其特征在于，在所述步驟2中，首先依據(jù)下述公式計(jì)算第k=1個(gè)界面的各個(gè)炮檢距的反射波時(shí)間：

t1,j(2)=(xj2+4h12)1/2/v1,j=1,2,...,N;]]>

其中，x_j為第k=1個(gè)界面上的N個(gè)由小到大排列的炮檢距參數(shù)，v_k為第k=1個(gè)界面上的層速度，h_k為第k=1個(gè)界面的層厚度。