[发明专利]沿起伏地表的地震数据重构方法、装置、介质及设备

权利要求书

1.一种沿起伏地表的地震数据重构方法，其特征在于，包括：

获取地震数据并确定待插值点空间坐标，在所有待插值点补齐地震道，并使补齐后的所述地震数据的道数与完整的待重建地震数据的地震道数相同；

利用时-深转换关系将起伏地表曲面深度表示成时间的形式，将所述地震数据区域延拓成二维平面数据体或三维立体数据体；

对补齐后的所述地震数据，利用曲线、曲面傅里叶变换进行傅里叶变换，根据需要选择适当的滤波函数，实现在频率-波数域中的稀疏地震数据的滤波，将在频率-波数域中滤波后的稀疏地震数据利用曲线、曲面傅里叶变换做逆傅里叶变换并进行迭代计算，得到变换数据；

将所述变换数据恢复成分布于起伏地表曲面的地震记录表示，得到二维和三维地段地表曲面的地震数据。

2.根据权利要求1所述的沿起伏地表的地震数据重构方法，其特征在于，所述方法还包括：根据所述地震数据确定采用二维插值重构或三维插值重构方法。

3.根据权利要求1所述的沿起伏地表的地震数据重构方法，其特征在于，

基于地震数据在频率-波数域中稀疏性特点，建立二维曲线和三维曲面傅里叶变换了的沿起伏地表地震数据重构迭代公式。

4.根据权利要求1至3任一所述的沿起伏地表的地震数据重构方法，其特征在于，确定先验信息的目标函数，所述先验信息的目标函数表示为

J_min=||d^obs-Sd||₂²⁺μ||d^obs-Sd||₁

其中，d是完整的待重建地震数据，S为采样矩阵，d^obs观测到的已有地震数据，||·||₂和||·||₁分别表示L-2及L-1范数,μ为调节前后两项权重的参数。

5.根据权利要求1至3任一所述的沿起伏地表的地震数据重构方法，其特征在于，对于含缺失地震道的观测地震数据d^obs (x,y,t)，分别利用二维快速傅里叶变换和三维快速傅里叶变换的正变换得到其F-k谱D^obs (f,k_x,k_y)；

对在频率-波数域中的稀疏数据根据需要进行适当的滤波处理，设滤波函数为W(f,k_x,k_y)，最小化目标函数公式则得到频率-波数域的最优化反演迭代公式为

d^k＝(I-S)F^-1WFd^k-1+d^obs，

k是迭代次数，d^k是第k次迭代后得到的插值重构结果，I是单位矩阵，S是采样矩阵，F及F^-1分别是正反曲线曲面傅里叶变换算子,d^obs为观测地震数据。

6.根据权利要求1至3任一所述的沿起伏地表的地震数据重构方法，其特征在于，由曲面方程z=F_S1(x,y)计算各网格点的坐标，利用曲面S₁、S₂方程的坐标关系，截取两曲线S₁～S₂之间的数据作为沿曲面S₁的三维傅里叶变换数据，以在沿曲面S₁上任意一点的正反傅里叶变换，可用于实现沿曲面S₁实际地震记录的数据重构。