[發(fā)明專利]一種海上地震資料中側(cè)反射波的壓制方法和裝置

有效

• 申請?zhí)枺?/td>: 201410014290.5
• 申請日: 2014-01-13
• 公開（公告）號: CN104777015B
• 公開（公告）日: 2018-06-01
• 發(fā)明（設(shè)計）人: 王童奎;路鵬飛;郭愛華;李瑩;孫明;何青城;李浩宇
• 申請（專利權(quán)）人: 中國石油天然氣股份有限公司
• 主分類號: G01N1/28
• 分類號: G01N1/28;G01N1/30
• 代理公司: 北京三高永信知識產(chǎn)權(quán)代理有限責(zé)任公司 11138
• 代理人: 劉映東

1.一種海上地震資料中側(cè)反射波的壓制方法，其特征在于，所述方法為：

解析海上地震資料，得到海上地震總反射波信息，所述海上地震總反射波信息至少包括：所述側(cè)反射波的旅行時、所述側(cè)反射波在海水中的傳播速度、所述側(cè)反射波的入射角、偏移距、有效反射波在海水中的傳播速度和炮源的地理坐標(biāo)；

根據(jù)所述海上地震總反射波信息，計算所述海上地震資料中側(cè)反射源的地理坐標(biāo)；

根據(jù)所述海上地震總反射波信息，通過繞射波的波場方程模擬所述地理坐標(biāo)處的側(cè)反射波的波場；

在所述海上地震總反射波的波場上，運用匹配相減法去除模擬的側(cè)反射波的波場，得到處理后的海上地震反射波；

所述海上地震資料通過使用炮源-側(cè)反射源-檢波器探測系統(tǒng)采集得到；

所述繞射波的波場方程為：

u 1 ( r S , t ) = ∫ V 0 v 1 ( r ) v 0 2 ( r ) u 0 2 ( r , t ) d V ]]>

其中，u1(rS,t)為所述側(cè)反射波的波場、rS為所述側(cè)反射源的地理坐標(biāo)、u0(r,t)為有效反射波的波場、r為有效反射波的波場中任一點的地理坐標(biāo)、t為側(cè)反射波的旅行時、v1(r)為側(cè)反射波在海水中的傳播速度、v0(r)為有效反射波在海水中的傳播速度、V為波場；

所述根據(jù)所述海上地震總反射波信息，計算所述海上地震資料中側(cè)反射源的地理坐標(biāo)包括：

根據(jù)所述側(cè)反射波的旅行時、所述側(cè)反射波在海水中的傳播速度、所述側(cè)反射波的入射角和偏移距，通過側(cè)反射波的時距曲線方程確定炮源到所述側(cè)反射源的距離；

根據(jù)所述炮源到所述該側(cè)反射源的距離和所述炮源的地理坐標(biāo)，通過側(cè)反射源地理坐標(biāo)的求解方程，確定所述側(cè)反射源的地理坐標(biāo)；

所述側(cè)反射波的時距曲線方程為：

t = O ′ O v + ( O ′ O ) 2 + x 2 - 2 O ′ O × x × c o s ( α ) v ]]>

其中，O'O為所述炮源到所述側(cè)反射源的距離、t為所述側(cè)反射波的旅行時、x為所述偏移距、v為所述側(cè)反射波在海水中的傳播速度、α為所述側(cè)反射波的入射角；

所述側(cè)反射源地理坐標(biāo)的求解方程為：

P 1 2 x + x 2 - 2 P 1 x × x + P 1 2 y + y 2 - 2 y × P 1 y = O ′ O 1 2 ]]>

P 2 2 x + x 2 - 2 x × P 2 x + P 2 2 y + y 2 - 2 y × P 2 y = O ′ O 2 2 ]]>

其中，x、y分別為該側(cè)反射源地理坐標(biāo)的橫、縱坐標(biāo)點；O'O1、O'O2為不同炮源到該側(cè)反射源的距離；P1x、P2x為不同炮源的地理坐標(biāo)的橫坐標(biāo)；P1y、P2y為不同炮源的地理坐標(biāo)的縱坐標(biāo)。

2.一種用于權(quán)利要求1所述海上地震資料中側(cè)反射波的壓制方法的海上地震資料中側(cè)反射波的壓制裝置，其特征在于，所述裝置包括：

解析模塊，用于解析海上地震資料，得到海上地震總反射波信息，所述海上地震總反射波信息至少包括：所述側(cè)反射波的旅行時、所述側(cè)反射波在海水中的傳播速度、所述側(cè)反射波的入射角、偏移距、有效反射波在海水中的傳播速度和炮源的地理坐標(biāo)；

計算模塊，用于根據(jù)所述海上地震總反射波信息，計算所述海上地震資料中側(cè)反射源的地理坐標(biāo)；

模擬模塊，用于根據(jù)所述海上地震總反射波信息，通過繞射波的波場方程模擬所述地理坐標(biāo)處的側(cè)反射波的波場；

處理模塊，用于在所述海上地震總反射波的波場上，運用匹配相減法去除模擬的側(cè)反射波的波場，得到處理后的海上地震反射波；

所述海上地震資料通過使用炮源-側(cè)反射源-檢波器探測系統(tǒng)采集得到；

所述繞射波的波場方程為：

u 1 ( r S , t ) = ∫ V 0 v 1 ( r ) v 0 2 ( r ) u 0 2 ( r , t ) d V ]]>

其中，u1(rS,t)為所述側(cè)反射波的波場、rS為所述側(cè)反射源的地理坐標(biāo)、u0(r,t)為有效反射波的波場、r為有效反射波的波場中任一點的地理坐標(biāo)、t為側(cè)反射波的旅行時、v1(r)為側(cè)反射波在海水中的傳播速度、v0(r)為有效反射波在海水中的傳播速度、V為波場；

所述計算模塊用于根據(jù)所述側(cè)反射波的旅行時、所述側(cè)反射波在海水中的傳播速度、所述側(cè)反射波的入射角和該偏移距，通過側(cè)反射波的時距曲線方程確定炮源到所述側(cè)反射源的距離；根據(jù)所述炮源到所述該側(cè)反射源的距離和所述炮源的地理坐標(biāo)，通過側(cè)反射源地理坐標(biāo)的求解方程，確定所述側(cè)反射源的地理坐標(biāo)；

所述側(cè)反射波的時距曲線方程為：

t = O ′ O v + ( O ′ O ) 2 + x 2 - 2 O ′ O × x × c o s ( α ) v ]]>

其中，O'O為所述炮源到所述側(cè)反射源的距離、t為所述側(cè)反射波的旅行時、x為所述偏移距、v為所述側(cè)反射波在海水中的傳播速度、α為所述側(cè)反射波的入射角；

側(cè)反射源地理坐標(biāo)的求解方程為：

P 1 2 x + x 2 - 2 P 1 x × x + P 1 2 y + y 2 - 2 y × P 1 y = O ′ O 1 2 ]]>

P 2 2 x + x 2 - 2 x × P 2 x + P 2 2 y + y 2 - 2 y × P 2 y = O ′ O 2 2 ]]>

其中，x、y分別為該側(cè)反射源地理坐標(biāo)的橫、縱坐標(biāo)點；O'O1、O'O2為不同炮源到該側(cè)反射源的距離；P1x、P2x為不同炮源的地理坐標(biāo)的橫坐標(biāo)；P1y、P2y為不同炮源的地理坐標(biāo)的縱坐標(biāo)。