[发明专利]频率波数域高角度断裂提取方法、装置、存储介质及设备

权利要求书

1.一种频率波数域高角度断裂提取方法，其特征在于，包括如下步骤：

a测量地震裂缝以获得原始地震数据体，将原始地震数据体拆分成彼此相邻的多个子块体，每个子块体的大小为N_x×N_y×N_z个采样点；

b对所述步骤a中的每个子块体进行f-k变换和f-k反变换处理，所述步骤b包括如下具体步骤：

b1对所述步骤a中的每个子块体进行f-k变换；

b2在频率波数域，滤除所述步骤b1中每个子块体进行f-k变换后与高角度断裂无关的低角度信息，求取G₁(k_x,k_y,ω)；

b3对所述步骤b2中的G₁(k_x,k_y,ω)进行f-k反变换；

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]为三维数据块体的坐标范围，G(k_x,k_y,ω)表示对f(x,y,z)进行f-k变换的结果，k_x∈[-k_xh,k_xh]、k_y∈[-k_yh,k_yh]、ω∈[-ω_h,ω_h]，k_xh、k_yh、ω_h的取值要求是包含G(k_x,k_y,ω)的有效信息，为虚数单位；N_x、N_y和N_z分别代表x轴、y轴、z轴三个方向上的采样点数量；

所述步骤b2包括如下具体步骤：求取G₁(k_x,k_y,ω)，使得

其中，θ_max为滤除的低角度信息的最大倾角；

所述步骤b3包括如下具体步骤：

对于一个三维数据块体G₁(k_x,k_y,ω)做如下变换：

其中，x∈[-x_h,x_h]、y∈[-y_h,y_h]、z∈[-z_h,z_h]。

3.根据权利要求2所述的频率波数域高角度断裂提取方法，其特征在于，所述步骤c包括如下具体步骤：

将所述步骤b处理后的每个子块体结果f₁(x,y,z)，按照其处理前的每个子块体f(x,y,z)在原始地震数据体中的位置进行拼合，到高角度断裂提取结果。