[發(fā)明專利]一種地質(zhì)條件勘探方法

有效

• 申請?zhí)枺?/td>: 201710953235.6
• 申請日: 2017-10-13
• 公開（公告）號: CN107748388B
• 公開（公告）日: 2019-05-07
• 發(fā)明（設(shè)計）人: 林羿;孫偉家;符力耘;魏偉
• 申請（專利權(quán)）人: 中國科學(xué)院地質(zhì)與地球物理研究所
• 主分類號: G01V1/32
• 分類號: G01V1/32;G01V1/34;G01V1/36
• 代理公司: 北京高沃律師事務(wù)所 11569
• 代理人: 王戈

1.一種地質(zhì)條件勘探方法，其特征在于，所述方法包括：

獲取笛卡爾坐標(biāo)系下地下介質(zhì)傳播速度v(x)、檢波器排列的地震記錄P(x,t)和震源子波的主頻頻帶范圍；

確定一個扁球坐標(biāo)系，所述扁球坐標(biāo)系包括ξ₁、ξ₂、ξ₃三個坐標(biāo)軸，將所述笛卡爾坐標(biāo)系下的地下介質(zhì)傳播速度v(x)和檢波器排列的地震記錄P(x,t)插值到所述扁球坐標(biāo)系中，獲取扁球坐標(biāo)系下的地下介質(zhì)傳播速度模型v(ξ)和檢波器排列的地震記錄P(ξ,t)；

將所述扁球坐標(biāo)系下的所述檢波器排列的地震記錄P(ξ,t)通過傅里葉變換至頻率域的檢波器排列的地震記錄P(ξ,ω)；

在扁球坐標(biāo)系下，對單個頻率的檢波器波場沿ξ₃坐標(biāo)軸方向?qū)⒉▓鱿蛳卵油匾粋€Δξ₃，通過速度時移和有限差分校正獲取ξ₃+Δξ₃處的檢波器波場直至將波場延拓到模型的底部ξ₃′處；

將所述震源子波的主頻頻帶范圍內(nèi)的不同頻率的波場均延拓到模型的底部ξ₃′處，疊加所有頻率的偏移量，得到扁球坐標(biāo)系的疊加后偏移量P(ξ)；

將扁球坐標(biāo)系的疊加后偏移量P(ξ)插值到笛卡爾坐標(biāo)系，得到最終的偏移剖面P(x)；

所述對單個頻率的檢波器波場沿ξ₃坐標(biāo)軸方向通過速度時移，具體包括：所述速度時移通過計算得到，

其中ξ為扁球坐標(biāo)系下一點的空間位置，(ξ₁，ξ₂，ξ₃)為扁球坐標(biāo)系下ξ點的坐標(biāo)，ξ_i(i＝1,2,3)為扁球坐標(biāo)系下ξ點的坐標(biāo)變量，P(ξ,ω)為ξ點的波場，ω為圓頻率，v(ξ)為ξ點的速度，A為控制扁球坐標(biāo)系寬度的常量，i為虛數(shù)單位；

所述有限差分校正通過計算得到，

其中，ξ_i(i＝1,2,3)為扁球坐標(biāo)系下ξ點的坐標(biāo)變量，ξ＝(ξ₁，ξ₂，ξ₃)為扁球坐標(biāo)系下一點的空間位置，P(ξ,ω)為ξ點的波場，ω為圓頻率，v(ξ)為ξ點的速度，A為控制扁球坐標(biāo)系寬度的常量,a_i和b_i為無量綱的常數(shù)，i為虛數(shù)單位。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述笛卡爾坐標(biāo)系下地下介質(zhì)傳播速度v(x)和檢波器排列的地震記錄P(x,t)中的x為笛卡爾坐標(biāo)系下一點的空間位置。

3.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述扁球坐標(biāo)系下的偏移量為P(ξ)＝∑_ωRe[P(ξ,ω)]。

4.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述偏移剖面P(x)為根據(jù)扁球坐標(biāo)系下最終偏移量P(ξ)得到的地震波場，經(jīng)插值變換到原笛卡爾坐標(biāo)系的波場。