[發(fā)明專利]基于傾角約束的多道反褶積處理方法

在審

• 申請?zhí)枺?/td>: 202011045519.3
• 申請日: 2020-09-29
• 公開（公告）號: CN112162318A
• 公開（公告）日: 2021-01-01
• 發(fā)明（設(shè)計）人: 王磊;劉春強;馬雄;黃保;趙妮
• 申請（專利權(quán)）人: 地球脈動（無錫）科技有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 無錫市匯誠永信專利代理事務(wù)所(普通合伙) 32260
• 代理人: 王晨光

1.基于傾角約束的多道反褶積處理方法，其特征在于：包括以下步驟：

S1：根據(jù)原始地震數(shù)據(jù)d(t,x)計算地層傾角θ(t,x)；

S2：根據(jù)地層傾角θ(t,x)計算反射系數(shù)r(t,x)沿傾角方向的導(dǎo)數(shù)

S3：對反射系數(shù)r(t,x)沿傾角方向的導(dǎo)數(shù)進行L2范數(shù)約束，構(gòu)建具有傾角約束的多道反褶積目標(biāo)函數(shù)J(r)；

S4：求解反褶積目標(biāo)函數(shù)J(r)，得到反射系數(shù)迭代求解公式。

2.根據(jù)權(quán)利要求1所述的基于傾角約束的多道反褶積處理方法，其特征在于：步驟S1具體如下：

S10：計算原始地震數(shù)據(jù)d(t,x)沿時間方向的偏導(dǎo)數(shù)

S11:計算原始地震數(shù)據(jù)d(t,x)沿空間方向的偏導(dǎo)數(shù)

S12：按照公式：計算地層傾角θ(t,x)。

3.根據(jù)權(quán)利要求1所述的基于傾角約束的多道反褶積處理方法，其特征在于：步驟S2具體如下：

S20:根據(jù)地層傾角θ(t,x)計算反射系數(shù)r(t,x)沿時間方向的偏導(dǎo)數(shù)在傾角方向的投影

S21：根據(jù)地層傾角θ(t,x)計算反射系數(shù)r(t,x)沿空間方向的偏導(dǎo)數(shù)在傾角方向的投影

S22：根據(jù)反射系數(shù)r(t,x)沿時間方向的偏導(dǎo)數(shù)在傾角方向的投影和沿空間方向的偏導(dǎo)數(shù)在傾角方向的投影計算反射系數(shù)r(t,x)沿傾角方向的導(dǎo)數(shù)

4.根據(jù)權(quán)利要求3所述的基于傾角約束的多道反褶積處理方法，其特征在于：步驟S20中，反射系數(shù)r(t,x)沿時間方向的偏導(dǎo)數(shù)在傾角方向的投影步驟S21中，反射系數(shù)r(t,x)沿空間方向的偏導(dǎo)數(shù)在傾角方向的投影步驟S22中，反射系數(shù)r(t,x)沿傾角方向的導(dǎo)數(shù)

5.根據(jù)權(quán)利要求1所述的基于傾角約束的多道反褶積處理方法，其特征在于：所述反褶積目標(biāo)函數(shù)其中，λ為縱向的調(diào)節(jié)因子，μ為傾角方向的調(diào)節(jié)因子。

6.根據(jù)權(quán)利要求5所述的基于傾角約束的多道反褶積處理方法，其特征在于：λ的取值范圍為0.01～0.5，μ的取值范圍為0.01～1。

7.根據(jù)權(quán)利要求5所述的基于傾角約束的多道反褶積處理方法，其特征在于：對反褶積目標(biāo)函數(shù)的求解過程如下：先將反射系數(shù)r(t,x)的初始值設(shè)置為0，然后利用迭代重加權(quán)算法求解目標(biāo)函數(shù)J(r),得到反射系數(shù)迭代求解公式：其中,r^k+1為經(jīng)過第(k+1)步迭代后的離散反射系數(shù)序列，W為多道地震子波矩陣，為迭代更新算子，D_θ為傾角約束算子。