[發明專利]基于上下行波場分解的逆時偏移方法及系統

在審

• 申請號: 201710633034.8
• 申請日: 2017-07-28
• 公開（公告）號: CN109307890A
• 公開（公告）日: 2019-02-05
• 發明（設計）人: 許璐;李博;王小青
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種基于上下行波場分解的逆時偏移方法，其特征在于，該方法包括：

構建震源與地震記錄的解析波場；

對解析波場進行波場的正向與反向延拓；

在成像網格沿深度方向對每個傳播時刻的波場進行希爾伯特變換，得到變換后的變換波場；

由原始波場與所述變換波場求解上下行波場；

選擇特定方向的波場利用互相關成像條件進行成像。

2.根據權利要求1所述的基于上下行波場分解的逆時偏移方法，其特征在于，由震源與地震記錄原始波場及其希爾伯特變換結果，構建解析波場：

式中，其中為解析波場，u為原始波場，u^h為u的希爾伯特變換結果。

3.根據權利要求1所述的基于上下行波場分解的逆時偏移方法，其特征在于，基于以下波動方程組實現對解析波場的正向與反向波場延拓：

式中，為震源或地震記錄，為震源或地震記錄的希爾伯特變換，為傳播過程中原始波場及其希爾伯特變換波場，為地震波速度。

4.根據權利要求3所述的基于上下行波場分解的逆時偏移方法，其特征在于，對解析波場進行波場的正向與反向延拓，得到每個傳播時刻的原始波場與沿時間方向做完希爾伯特變換后的變換波場。

5.根據權利要求1所述的基于上下行波場分解的逆時偏移方法，其特征在于，在波場傳播的每個時刻，沿成像網格深度方向對u^h進行z方向的希爾伯特變換，得到每個傳播時刻做完時間與深度希爾伯特變換后的變換波場：

式中，u^h為u的希爾伯特變換結果，H_tH_z(u(x,t))為時間方向與深度方向的希爾伯特變換波場。

6.根據權利要求1所述的基于上下行波場分解的逆時偏移方法，其特征在于，由原始波場與變換波場求解上下行波場，其表達式為：

u_down(x,t)＝[u(x,t)+H_tH_z(u(x,t))]/2

u_up(x,t)＝[u(x,t)-H_tH_z(u(x,t))]/2

式中，H_t與H_z分別表示在時間方向與深度方向上的希爾伯特變換，u(x,t)為原始傳播波場，u_down(x,t)、u_up(x,t)為下行波場與上行波場。

7.根據權利要求6所述的基于上下行波場分解的逆時偏移方法，其特征在于，基于上下行波場與原始波場以及變換波場間的關系式，完成炮點與檢波點的上下行波場分離。

8.根據權利要求1所述的基于上下行波場分解的逆時偏移方法，其特征在于，利用傳播方向相反波場進行互相關成像，成像條件為：

式中，S_u(t,m)、S_d(t,m)為震源上行與下行波場；R_u(t,m)、R_d(t,m)為檢波端上行與下行波場。

9.一種基于上下行波場分解的逆時偏移系統，其特征在于，該系統包括：

存儲器，存儲有計算機可執行指令；

處理器，所述處理器運行所述存儲器中的計算機可執行指令，執行以下步驟：

構建震源與地震記錄的解析波場；

對解析波場進行波場的正向與反向延拓；

在成像網格沿深度方向對每個傳播時刻的波場進行希爾伯特變換，得到變換后的變換波場；

由原始波場與所述變換波場求解上下行波場；

選擇特定方向的波場利用互相關成像條件進行成像。

10.根據權利要求9所述的基于上下行波場分解的逆時偏移系統，其特征在于，對炮點波場與檢波點波場進行上下行波分離后，只對傳播方向相反的波場進行相關成像。