[發明專利]基于反射結構約束的多道吸收補償處理方法

在審

• 申請號: 202110078488.X
• 申請日: 2021-01-20
• 公開（公告）號: CN112904416A
• 公開（公告）日: 2021-06-04
• 發明（設計）人: 趙東鳳;李國發;李皓;馬雄;張思超;高洋;楊娜霞;丁禹行;梁伽福;張平民;彭建偉;李佼佼
• 申請（專利權）人: 中國石油大學（北京）
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京紀凱知識產權代理有限公司 11245
• 代理人: 李曉紅

1.一種基于反射結構約束的多道吸收補償處理方法，其特征在于，包括以下步驟：

S1：根據原始地震數據d(t,x)計算反射結構表征算子r_a,b；

S2：根據反射結構表征算子r_a,b構造反射結構約束正則化項R(t,x)；

S3：將反射結構約束正則化項R(t,x)加入到吸收補償目標函數中，構建具有反射結構約束的多道吸收補償目標函數J(s)；

式中，s(t,x)為無衰減的地震記錄；a(τ,t)為地層吸收衰減函數；λ為縱向的調節因子；μ為反射結構約束的調節因子，表示褶積運算符；

S4：采用共軛梯度法對基于反射結構約束的多道吸收補償目標函數進行求解，得到吸收補償后的地震數據；

s＝(A^TA+λI+μR^TR)^-1A^Td，

式中，A為吸收衰減褶積矩陣，是由a(τ,t)離散得到的矩陣；I為單位矩陣；R是反射結構約束褶積算子，是由r_a,b離散構造得到的矩陣；d為含噪聲衰減地震數據向量；是由d(t,x)離散得到；s為吸收補償后的地震數據，即無衰減的地震記錄。

2.如權利要求1所述的基于反射結構約束的多道吸收補償處理方法，其特征在于，在上述步驟S1中，反射結構表征算子r_a,b按照如下公式計算：

式中，d_a,b(t,x)表示從原始地震數據d(t,x)中沿時間方向平移a、空間方向平移b截取的子數據；M和N分別為反射結構表征算子r_a,b在水平方向和垂直方向的大小。

3.如權利要求1所述的基于發射結構約束的多道吸收補償處理方法，其特征在于，在上述步驟S2中，反射結構約束正則化項R(t,x)按照下述公式來構建：

4.如權利要求1所述的基于反射結構約束的多道吸收補償處理方法，其特征在于，在上述步驟S3中，吸收補償目標函數如下：

5.如權利要求1所述的基于反射結構約束的多道吸收補償處理方法，其特征在于：在上述步驟S4中，A為由a(τ,t)通過時間離散采樣得到的矩陣，R為由r_a,b通過時間離散采樣構造得到的矩陣；d是由d(t,x)通過時間離散采樣得到。

6.如權利要求1所述的基于反射結構約束的多道吸收補償處理方法，其特征在于，λ的取值范圍為0.02～1；μ的取值范圍為0.01～2。

7.如權利要求2所述的基于反射結構約束的多道吸收補償處理方法，其特征在于，M取5～9，N取5～7。