[发明专利]抗假频的地震数据插值方法和装置

权利要求书

1.一种抗假频的地震数据插值方法，其特征在于，包括：

获取目标区域的时间域空间域的地震数据体；

根据所述时间域空间域的地震数据体，获取多个频率域空间域的单频地震数据体；

利用带有空间平滑算子的插值目标函数，根据所述多个频率域空间域的单频地震数据体中的各个频率域空间域的单频地震数据体，确定对应的插值后的频率域波数域的单频地震数据体；

根据所述插值后的频率域波数域的单频地震数据体，获取插值后的时间域空间域的地震数据体；

其中，利用带有空间平滑算子的插值目标函数，根据所述多个频率域空间域的单频地震数据体中的各个频率域空间域的单频地震数据体，确定对应的插值后的频率域波数域的单频地震数据体，包括：

根据所述多个频率域空间域的单频地震数据体中的各个频率域空间域的单频地震数据体，建立对应频率的带有空间平滑算子的插值目标函数；

分别求解各个对应频率的带有空间平滑算子的插值目标函数，得到对应的插值后的频率域波数域的单频地震数据体；

其中，建立的对应频率的带有空间平滑算子的插值目标函数为：

上式中，d为频率域空间域的单频地震数据体的向量形式，m为插值后的频率域波数域的单频地震数据体的向量形式，G为空间坐标的不等间隔的逆傅里叶变换矩阵，λ_W为第一权重因子，λ_S为第二权重因子，||m||_W²为权值模算子，||m||_S²为空间平滑算子，Minimize为求最小值。

2.根据权利要求1所述的方法，其特征在于，根据所述时间域空间域的地震数据体，获取多个频率域空间域的单频地震数据体，包括：

对所述时间域空间域的地震数据体进行傅里叶变换，得到频率域空间域的地震数据体；

将所述频率域空间域的地震数据体按照频率划分为多个所述频率域空间域的单频地震数据体。

3.根据权利要求1所述的方法，其特征在于，所述空间平滑算子按照以下公式确定：

上式中，||m||_S²为空间平滑算子，m为插值后的频率域波数域的单频地震数据体的向量形式，m^H为m的共轭转置矩阵，D_η为空间差分矩阵，η为空间维度，F^-1为空间坐标的逆傅里叶变换矩阵。

4.根据权利要求1所述的方法，其特征在于，所述第一权重因子的取值范围为大于等于0.9，且小于1。

5.根据权利要求4所述的方法，其特征在于，所述第二权重因子的取值范围为大于0，且小于等于0.1。

6.根据权利要求1所述的方法，其特征在于，分别求解各个对应频率的带有空间平滑算子的插值目标函数，得到对应的插值后的频率域波数域的单频地震数据体，包括：

通过共轭梯度算法，分别求解对应频率的带有空间平滑算子的插值目标函数，得到对应的插值后的频率域波数域的单频地震数据体。

7.根据权利要求1所述的方法，其特征在于，根据所述插值后的频率域波数域的单频地震数据体，获取插值后的时间域空间域的地震数据体，包括：

对所述插值后的频率域波数域的单频数据体进行傅里叶变换，得到插值后的频率域空间域的单频地震数据体；

合并所述插值后的频率域空间域的单频地震数据体，得到插值后的频率域空间域地震数据体；

对所述插值后的频率域空间域地震数据体进行傅里叶变换，得到所述插值后的时间域空间域的地震数据体。

8.根据权利要求1所述的方法，其特征在于，在获取所述插值后的时间域空间域的地震数据体后，所述方法还包括：

根据所述插值后的时间域空间域的地震数据体，对所述目标区域进行地质勘探。

9.一种抗假频的地震数据插值装置，其特征在于，包括：

第一获取模块，用于获取目标区域的时间域空间域的地震数据体；

第二获取模块，用于根据所述时间域空间域的地震数据体，获取多个频率域空间域的单频地震数据体；

确定模块，用于利用带有空间平滑算子的插值目标函数，根据所述多个频率域空间域的单频地震数据体中的各个频率域空间域的单频地震数据体，确定对应的插值后的频率域波数域的单频地震数据体；

第三获取模块，用于根据所述插值后的频率域波数域的单频地震数据体，获取插值后的时间域空间域的地震数据体；

其中，所述确定模块包括：

建立单元，用于根据所述多个频率域空间域的单频地震数据体中的各个频率域空间域的单频地震数据体，建立对应频率的带有空间平滑算子的插值目标函数；

求解单元，用于分别求解各个对应频率的带有空间平滑算子的插值目标函数，得到对应的插值后的频率域波数域的单频地震数据体；

其中，所述建立单元具体用于按照以下公式建立对应频率的带有空间平滑算子的插值目标函数：

上式中，d为频率域空间域的单频地震数据体的向量形式，m为插值后的频率域波数域的单频地震数据体的向量形式，G为空间坐标的不等间隔的逆傅里叶变换矩阵，λ_W为第一权重因子，λ_S为第二权重因子，||m||_W²为权值模算子，||m||_S²为空间平滑算子，Minimize为求最小值。