[發(fā)明專利]拉東域高角度斷裂提取方法、裝置、存儲介質及設備

有效

• 申請?zhí)枺?/td>: 202110393732.1
• 申請日: 2021-04-13
• 公開（公告）號: CN113126156B
• 公開（公告）日: 2023-02-24
• 發(fā)明（設計）人: 馬淑芳;范廷恩;范洪軍;樊鵬軍;牛濤;周建楠
• 申請（專利權）人: 中海石油（中國）有限公司;中海石油（中國）有限公司北京研究中心
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京紀凱知識產權代理有限公司 11245
• 代理人: 任文娟

1.一種拉東域高角度斷裂提取方法，其特征在于，包括如下步驟：

a測量地震裂縫以獲得原始地震數據體，將原始地震數據體拆分成彼此相鄰的多個子塊體，每個子塊體的大小為N_x×N_y×N_z個采樣點；

b對所述步驟a中的每個子塊體進行處理，具體步驟如下：

b1對所述步驟a中的每個子塊體進行拉東變換；

b2在拉東域，濾除所述步驟b1中每個子塊體進行拉東變換后與高角度斷裂無關的低角度信息；

b3增強所述步驟b1中每個子塊體進行拉東變換后與高角度斷裂有關的高角度信息；

b4對所述步驟b3中的高角度信息進行拉東反變換；

c將所述步驟b處理后的多個子塊體進行拼合，得到高角度斷裂提取結果。

2.根據權利要求1所述的拉東域高角度斷裂提取方法，其特征在于，所述步驟b1包括如下具體步驟：

對一個三維數據塊體f(x,y,z)做如下變換：

其中，z為垂直方向坐標，x,y為兩個水平方向坐標，它們的取值范圍x∈[-x_h,x_h]、y∈[-y_h,y_h]、z∈[-z_h,z_h]為三維數據塊體的坐標范圍，表示沿著空間上一個平面對f(x,y,z)進行積分求和的結果，這個平面的具體參數由ρ,θ,表征，ρ為坐標原點到這個平面的距離，ρ∈[-ρ_h,ρ_h]，θ為平面法向量與z軸正方向的夾角，為平面法向量與x軸正方向的夾角N_x、N_y和N_z分別代表x軸、y軸、z軸三個方向上的采樣點數量；

所述步驟b2包括如下具體步驟：求取使得

其中，θ_max為濾除的低角度信息的最大傾角；

所述步驟b3包括如下具體步驟：求取使得

其中k∈[1,2]，k＝1表示不對高角度信息進行增強，k越大，表示對高角度信息的增強越顯著；

所述步驟b4包括如下具體步驟：

對于一個三維數據塊體做如下變換：

其中，x∈[-x_h,x_h]、y∈[-y_h,y_h]、z∈[-z_h,z_h]，ρ為坐標原點到這個平面的距離，ρ∈[-ρ_h,ρ_h]，θ為平面法向量與z軸正方向的夾角，為平面法向量與x軸正方向的夾角

3.根據權利要求1所述的拉東域高角度斷裂提取方法，其特征在于，所述步驟c包括如下具體步驟：

將所述步驟b處理后的每個子塊體結果f₁(x,y,z)，按照其處理前的每個子塊體f(x,y,z)在原始地震數據體中的位置進行拼合，到高角度斷裂提取結果。

4.一種如權利要求1-3任意一項所述拉東域高角度斷裂提取方法的提取裝置，其特征在于，包括：

第一處理單元a，用于測量地震裂縫以獲得原始地震數據體，將原始地震數據體拆分成彼此相鄰的多個子塊體，每個子塊體的大小為N_x×N_y×N_z個采樣點；

第二處理單元b，用于所述第一處理單元a中的每個子塊體進行處理，具體步驟如下：

b1對所述第一處理單元a中的每個子塊體進行拉東變換；

b2在拉東域，濾除所述步驟b1中每個子塊體進行拉東變換后與高角度斷裂無關的低角度信息；

b3增強所述步驟b1中每個子塊體進行拉東變換后與高角度斷裂有關的高角度信息；

b4對所述步驟b3中的高角度信息進行拉東反變換；

第三處理單元，用于將所述第二處理單元b處理后的多個子塊體進行拼合，得到高角度斷裂提取結果。

5.一種計算機可讀存儲介質，其上存儲有計算機程序，其特征在于，所述計算機程序被處理器執(zhí)行時實現權利要求1-3任意一項所述拉東域高角度斷裂提取方法的步驟。

6.一種計算機設備，包括存儲器、處理器及存儲在存儲器上并可在處理器上運行的計算機程序，其特征在于，所述處理器執(zhí)行所述計算機程序時實現權利要求1-3任意一項所述拉東域高角度斷裂提取方法的步驟。