[發明專利]一種基于復波場的方位角度域偏移成像方法及系統

在審

• 申請號: 202011114546.1
• 申請日: 2020-10-16
• 公開（公告）號: CN114442148A
• 公開（公告）日: 2022-05-06
• 發明（設計）人: 許璐;李博;呂慶達
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京知舟專利事務所(普通合伙) 11550
• 代理人: 郭韞

1.一種基于復波場的方位角度域偏移成像方法，其特征在于：所述方法首先進行炮點波場、檢波點波場的復波場的延拓，然后在各個空間窗內實現復波場的局部角度分解，獲得炮點波場、檢波點波場的局部角度分量，最后對不同方位的角度域道集進行疊加，獲得最終的方位角域成像結果。

2.根據權利要求1所述的基于復波場的方位角度域偏移成像方法，其特征在于:所述方法包括：

步驟一：輸入數據：震源、震源HT、炮集、炮集HT；

步驟二：進行炮點復波場外推得到炮點復波場，同時進行檢波點復波場外推得到檢波點復波場；

步驟三：分別將炮點復波場、檢波點復波場分解為多個相同的空間窗，然后獲得每一個空間窗的所有角度分量；

步驟四，利用所有空間窗的角度分量獲得方位角域成像結果。

3.根據權利要求2所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟二的操作包括：

利用下面的公式(1)進行炮點復波場外推得到炮點復波場，利用下面的公式(2)進行檢波點復波場外推得到檢波點復波場：

其中，Us和Ur分別為炮端波場和檢波端波場，v為速度，t為時間，fs為震源子波，fs^h為做完Hilbert變換后的子波，dr和dr^h為炮集和做完Hilbert變換后的炮集，與為炮點復波場與檢波點復波場。

4.根據權利要求3所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟三的操作包括：

(1)：將當前空間窗內的復波場變換到波數域；

(2)：計算波數域內每個斜率上的疊加能量；

(3)：篩選出最大疊加能量角度分量；

(4)：將所述最大疊加能量角度分量反變換到空間域，并求得誤差泛函；

(5)：判斷所述誤差泛函的值是否小于設定的閾值，如果否，則轉入步驟(6)，如果是，則移動到下一個空間窗，然后返回步驟(1)；

(6)：輸出當前角度分量，然后轉入步驟(1)。

5.根據權利要求4所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟(1)的操作包括：

利用下面的公式(5)將當前空間窗內的復波場變換到波數域：

其中K_sx＝(k_sx,k_sz)為空間波數矢量，win為空間窗大小，包括橫向和縱向；t，x分別為步驟一輸入的數據中的時間和空間。

6.根據權利要求5所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟(3)的操作包括：

從步驟(2)計算得到的疊加能量中找到最大的疊加能量，該最大疊加能量對應的斜率即為最大能量角度分量。

7.根據權利要求6所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟(4)中的將所述最大疊加能量角度分量反變換到空間域的操作包括：

利用公式(6)將所述最大疊加能量角度分量進行傅立葉反變換：

其中，其中，代表局部空間波場，

θs表示自變量，θ₀為步驟(3)得到的最大疊加能量角度。

8.根據權利要求7所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟(4)中的求得誤差泛函的操作包括：

利用公式(4)計算得到誤差泛函：

9.根據權利要求8所述的基于復波場的方位角度域偏移成像方法，其特征在于：所述步驟四的操作包括：

利用下面的公式得到方位角域成像結果：

其中，

10.一種基于復波場的方位角度域偏移成像系統，其特征在于：所述系統包括：

輸入單元，用于輸入數據：震源、震源HT、炮集、炮集HT；

波場延拓單元，與所述輸入單元連接，用于進行炮點復波場外推得到炮點復波場，同時進行檢波點復波場外推得到檢波點復波場；

角度分量獲取單元，與所述波場延拓單元連接，用于分別將炮點復波場、檢波點復波場分解為多個相同的空間窗，然后獲得每一個空間窗的所有角度分量；

成像單元，與所述角度分量獲取單元連接，用于利用所有空間窗的角度分量獲得方位角域成像結果。