[发明专利]基于分频动态编码的可控震源混叠数据成像方法

权利要求书

1.基于分频动态编码的可控震源混叠数据成像方法，其特征在于，该基于分频动态编码的可控震源混叠数据成像方法包括：

步骤1，输入震源和炮记录；

步骤2，通过波场对震源和炮记录进行波场延拓后重建编码前的震源波场和接收波场；

步骤3，输入分频编码矩阵；

步骤4，应用编码矩阵对震源波场和接收波场进行编码；

步骤5，对编码后的震源波场和接收波场进行成像；以及

步骤6，输出成像结果。

2.根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在于，在步骤2中，重建的震源波场和接收波场为：

其中：分别代表通过波场对震源和炮记录进行波场延拓后重建的震源波场和接收波场，N_S表示炮记录的总数，W_s、W_R表示物理相关的震源波场和接收波场。

3.根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在于，在步骤3中，使用E_m,n表示编码矩阵，m表示行指数，即炮指数，n是列指数，即编码次数，把矩阵的元素表示成如下的形式：

E_m,n＝A_m,nexp(-jφ_m,n)(3)

A代表第n次编码中第m个炮集的振幅权重，φ_m,n代表第n次编码中第m个炮集的相移，假设编码矩阵为一个频率独立的编码矩阵，这个矩阵为三维的，第三维是频率ω，因此编码矩阵可以写成：

E_m,n(ω)＝A_m,n(ω)exp(-jφ_m,n(ω))(4)

分频编码的编码矩阵表示为：

其中w_i为第i个频率。

4.根据权利要求3所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在于，在步骤3中，分频编码的方式目的是在i≠j时使其中，i、j分别是编码矩阵的纵横坐标，分别代表通过波场对震源和炮记录进行波场延拓后重建的震源波场和接收波场；在对一个超炮道集进行编码时，将频率域的波场进行分频，在对使每个震源波场分别使用不同的频率段，且一个超道集中个震源使用的频率没有重叠，对应的接收波场频率分配情况与震源波场相同。

5.根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在于，在步骤4中，对炮集编码形成超炮道集用炮集与编码矩阵相乘来表示，那么编码后的震源波场和接收波场为：

B_S＝W_SE(5)

B_R＝W_RE(6)

其中：B_s、B_R分别为编码后的震源波场和接受波场，W_s、W_R分别为重建后的震源波场和接收波场，E为编码矩阵。

6.根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在于，在步骤5中，对编码后的震源波场和接收波场进行成像的表达式为：

其中，R_e为编码后的震源波场和接收波场的成像结果，B_s为编码后的震源波场，是编码后接收波场的共轭转置。