[发明专利]基于插值与去混叠同时处理的地震数据分离方法及系统

权利要求书

1.一种基于插值与去混叠同时处理的地震数据分离方法，其特征在于，该方法包括如下步骤：

(1)基于观测地震数据提取伪分离地震数据；

(2)基于Curvelet变换对伪分离地震数据进行插值与去混叠同时处理迭代更新直至迭代结束，输出主震源、副震源数据，完成地震数据分离；

步骤(1)具体为：

(11)确定混叠算子Γ；

(12)获取伪分离地震数据：F_t^HΓ^HF_td，其中，F_t、F_t^H分别为时间方向的一维正、反傅里叶变换，Γ^H为混叠算子Γ的共轭算子，d为观测地震数据；

步骤(2)具体为：

(21)赋值n＝1，其中，为主震源第n次迭代值，为副震源第n次迭代值；

(22)确定去混叠项和插值项并更新主震源和副震源数据：

其中，和为去混叠项，和为插值项，为主震源第n+1次迭代值，为副震源第n+1次迭代值，C、C^H为Curvelet正、反变换，T_λ为硬阈值函数，R为采样算子，且R满足可交换性，即RF_t＝F_tR，RΓ＝ΓR；

(23)判断n是否小于N，若是则赋值n＝n+1并返回步骤(22)，否则输出主震源和副震源数据，其中，N为最大迭代次数。

2.根据权利要求1所述的一种基于插值与去混叠同时处理的地震数据分离方法，其特征在于，所述的混叠算子Γ具体确定为：

其中，e为自然常数，i为虚数单位，ω为圆频率，t_k为副震源中第k炮的时间延迟，k＝1,...,M，M为副震源中炮的数量。

3.一种基于插值与去混叠同时处理的地震数据分离系统，该系统包括：

地震数据提取模块：该模块基于观测地震数据提取伪分离地震数据；

插值与去混叠分离模块：该模块对伪分离地震数据进行插值与去混叠并迭代更新，得到分离的主震源、副震源数据；

输出模块：该模块输出分离的主震源、副震源数据；

所述的地震数据提取模块获取的伪分离地震数据为F_t^HΓ^HF_td，其中，F_t、F_t^H分别为时间方向的一维正、反傅里叶变换，Γ^H为混叠算子Γ的共轭算子，d为观测地震数据；

所述的插值与去混叠分离模块执行以下步骤完成主震源、副震源数据分离：

(a)赋值n＝1，其中，为主震源第n次迭代值，为副震源第n次迭代值；

(b)确定去混叠项和插值项并更新主震源和副震源数据：

其中，和为去混叠项，和为插值项，为主震源第n+1次迭代值，为副震源第n+1次迭代值，C、C^H为Curvelet正、反变换，T_λ为硬阈值函数，R为采样算子，且R满足可交换性，即RF_t＝F_tR，RΓ＝ΓR；

(c)判断n是否小于N，若是则赋值n＝n+1并返回步骤(b)，否则输出主震源和副震源数据，其中，N为最大迭代次数。

4.根据权利要求3所述的一种基于插值与去混叠同时处理的地震数据分离系统，其特征在于，所述的混叠算子Γ具体确定为：

其中，e为自然常数，i为虚数单位，ω为圆频率，t_k为副震源中第k炮的时间延迟，k＝1,...,M，M为副震源中炮的数量。