[發明專利]基于Shearlet變換的多分量地震數據Corssline方向波場重建方法

在審

• 申請號: 201710310807.9
• 申請日: 2017-05-05
• 公開（公告）號: CN107121701A
• 公開（公告）日: 2017-09-01
• 發明（設計）人: 劉成明;王德利
• 申請（專利權）人: 吉林大學
• 主分類號: G01V1/32
• 分類號: G01V1/32;G01V1/34
• 代理公司: 長春吉大專利代理有限責任公司22201
• 代理人: 王立文

1.一種基于Shearlet變換的多分量地震數據Corssline方向波場重建方法，其特征在于，是將Shearlet變換與閾值迭代法相結合，地震數據插值歸結為l₁范數最優化方程：

上式中M為觀測矩陣，M中元素為0表示數據缺失，元素為1表示數據存在，S^-1為Shearlet逆變換；x為一組Shearlet系數；

將地震數據插值轉換到Shearlet域重構出地震數據，反演出一組l₁范數最小的Shearlet系數x，公式(1)將地震數據采集系統作為觀測系統，Shearlet變換作為壓縮感知理論實施的稀疏變換，構成了多分量地震數據插值的基礎，在多分量地震數據插值中，引入了Vy分量，由于

P和分量作為相互的約束條件，為彼此提供額外的插值信息，將P和同時作為插值輸入，壓力波場和加速度波場同時重建，將公式(1)與公式(2)相加得到：

觀測到的壓力和加速度分量作為已知輸入，x_p,x_v分別為重構出壓力波場和加速度波場的shearlet系數，每次迭代過程中P和相互更新，利用P的導數求取積分求取P。

2.一種基于Shearlet變換的多分量地震數據Corssline方向波場重建方法，其特征在于，包括以下步驟：

A、初始化迭代參數：迭代次數N,閾值因子L_max參數λ；

B、初始化壓力波場以及梯度波場P＝0，

C、計算初始迭代參數δ_p＝λ_p·L_max，δ_v＝λ_v·L_max；

D、主循環For i＝1 to N do；

E、利用基于Shearlet變換的地震數據重構方法重建壓力場和梯度場；

F、更新閾值參數δ_p andδ_v；

G、Ifλ_p≤λ_min&λ_v≤λ_min結束.Else返回D；

H、輸出:重構的壓力波場P以及速度波場

3.按照權利要求2所述的基于Shearlet變換的多分量地震數據Corssline方向波場重建方法，其特征在于，E步驟所述的利用基于Shearlet變換的地震數據重構方法重建壓力場，包括以下步驟：

①、通過更新P_i：

②、計算P波場的殘差resp_i＝M(P_M-P_i)

③、通過閾值算法，求解壓力波場的Shearlet系數

xi+1=Tλp(xi+S(Pi+respi)).]]>

④、通過Shearlet逆變換更新壓力波場P_i+1＝Sx_i+1。

4.按照權利要求2所述的基于Shearlet變換的多分量地震數據Corssline方向波場重建方法，其特征在于，E步驟所述的利用基于Shearlet變換的地震數據重構方法重建梯度場，包括以下步驟：

I、通過P_i更新

Ⅱ、計算加速度波場的殘差

Ⅲ、通過閾值算法，求解加速度波場的Shearlet系數

xi+1=Tλv(xi+S(Vgi+resvi)).]]>

Ⅳ、通過Shearlet逆變換更新加速度波場