[發(fā)明專利]基于共散射角道集的散射波分離方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201710512085.5
• 申請(qǐng)日: 2017-06-28
• 公開(kāi)（公告）號(hào): CN109143362B
• 公開(kāi)（公告）日: 2020-06-19
• 發(fā)明（設(shè)計(jì)）人: 白英哲
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號(hào): G01V1/34
• 分類號(hào): G01V1/34
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于共散射角道集的散射波分離方法，其特征在于，該方法包括：

通過(guò)共散射角偏移將數(shù)據(jù)偏移到散射角域，生成共散射角道集；

根據(jù)反射波和散射波在共散射角道集中的不同響應(yīng)特征分離散射波和反射波；

其中，通過(guò)以下公式生成共散射角道集：

其中，G表示共散射角道集，是空間散射入射角，是空間散射反射角，是空間坐標(biāo)，θ_s和φ_s分別為散射入射角和散射入射方位角，θ_r和φ_r分別為散射反射角和散射反射方位角，t是時(shí)間，t_max是最大傳播時(shí)間，S_min和S_max分別是最小和最大炮號(hào)，σ是小的實(shí)數(shù)，P_D和P_U分別是正傳波場(chǎng)和反傳波場(chǎng)；通過(guò)波動(dòng)方程形成波場(chǎng)傳播系統(tǒng)，計(jì)算正傳波場(chǎng)和反傳波場(chǎng)：

其中，P_D是正傳波場(chǎng)，P_U是反傳波場(chǎng)，S是震源函數(shù)，D是觀測(cè)數(shù)據(jù)，c是傳播速度，z是成像深度；根據(jù)正傳波場(chǎng)和反傳波場(chǎng)求取正傳波場(chǎng)和反傳波場(chǎng)的坡印廷矢量，

其中，和分別是正傳波場(chǎng)和反傳波場(chǎng)的坡印廷矢量；

通過(guò)以下公式求取散射角：

其中，θ_s和φ_s分別為散射入射角和散射入射方位角，θ_r和φ_r分別為散射反射角和散射反射方位角，(k_Sx,k_Sy,k_Sz)和(k_Rx,k_Ry,k_Rz)分別表示正傳波場(chǎng)和反傳波場(chǎng)的波數(shù)。

2.如權(quán)利要求1所述的散射波分離方法，其特征在于，通過(guò)傅里葉變換法求取正傳波場(chǎng)和反傳波場(chǎng)的波數(shù)，傅里葉變換為：

其中，F(xiàn)T[·]表示傅里葉變換，和分別表示正傳波場(chǎng)和反傳波場(chǎng)的傅里葉變換，(x,y,z)表示 空間坐標(biāo)，(k_Sx,k_Sy,k_Sz)和(k_Rx,k_Ry,k_Rz)分別表示正傳波場(chǎng)和反傳波場(chǎng)的波數(shù)。

3.如權(quán)利要求1所述的散射波分離方法，其特征在于，將傅里葉變換后的波場(chǎng)和重排到散射角域，得到和將重排后的波場(chǎng)變換回空間域得到帶有角度信息的空間波場(chǎng)：FT^-1[·]表示傅里葉反變換。

4.如權(quán)利要求1所述的散射波分離方法，其特征在于，利用反褶積成像條件作為保幅成像條件進(jìn)行擴(kuò)展，引入帶有角度的波場(chǎng)進(jìn)行成像，生成所述共散射角道集，反褶積成像條件表達(dá)式為：

其中，I是像，P_D和P_U分別是正傳波場(chǎng)和反傳波場(chǎng)。

5.如權(quán)利要求1所述的散射波分離方法，其特征在于，通過(guò)入射反射線共面，分布于法線兩側(cè)，并且入射角等于反射角三個(gè)判別準(zhǔn)則，獲得切除角度，通過(guò)沿所述切除角度進(jìn)行切除，實(shí)現(xiàn)消除反射波，保留散射波。

6.如權(quán)利要求5所述的散射波分離方法，其特征在于，利用三個(gè)判別準(zhǔn)則，獲得切除角度包括：

1)入射反射線和法線共面，有：

2)入射角等于反射角，有：

cosθ_scosθ_d+sinθ_ssinθ_dcos(φ_s-φ_d)

＝cosθ_rcosθ_d+sinθ_rsinθ_dcos(φ_r-φ_d)

3)入射角和反射角位于法線兩側(cè)，有：

其中，θ_s和φ_s分別為散射入射角和散射入射方位角，θ_r和φ_r分別為散射反射角和散射反射方位角，θ_d和φ_d分別是界面傾角和傾角方位角。