[0047] Εη，η(ω) = An，n(co)exp(-j φη η(ω)) (4)
[0048] 分频编码的编码矩阵可以表示为：
[0049]
[0050] 其中Wi为第i个频率。
[0051] 也可以理解为分频编码的方式目的是在i辛j时使
其中，i、j 分别是编码矩阵的纵横坐标，W,1分别代表通过波场对震源和炮记录进行波场延拓 后重建的震源波场和接收波场。为了达到这个目的在对一个超炮道集进行编码时，将频率 域的波场进行分频，在对使每个震源波场分别使用不同的频率段，且一个超道集中个震源 使用的频率没有重叠，对应的接收波场频率分配情况与震源波场相同。
[0052] 在步骤104,应用编码矩阵对震源波场和接收波场进行编码。
[0053] 对炮集编码形成超炮道集可以用炮集与编码矩阵相乘来表示，那么编码后的震源 波场和接收波场为：
[0054] Bs= WsE (5)
[0055] Br= WrE (6)
[0056] 其中：BS、BR分别为编码后的震源波场和接受波场，W S、WR分别为重建后的震源波场 和接收波场，E为编码矩阵。
[0057] 在步骤105,对编码后的震源波场和接收波场进行成像；
[0059] 其中，艮为编码后的震源波场和接收波场的成像结果，B 3为编码后的震源波场， /，是编码后接收波场的共辄转置。
[0060] 在步骤106,输出成像结果。流程结束。
[0061] 为了研究分频范围与成像结果的关系，我们分别对不同频率范围进行分频编码比 较。震源的主频为30赫兹，围绕这一主频我们分别对0-60hz、5-55hz、10-50hz、15-45hz、 20-40hz以及25-35hz六个频率范围进行分频编码测试，分析分频范围对成像结果的影响。 图3到图8是对不同频带进行分频编码成像的结果，从图中我们可以看出随着频率范围的 变小，成像结果的质量在逐渐变差，但是当频率范围变小的范围不大时并不明显，而且频率 范围的变小可以减小计算量从而减少计算时间。图9是不同频带进行测试的耗时与频带 范围的关系，可见随着频带长度变窄计算时间也是在逐渐缩短的。因此，我们可以在保证成 像质量的前提下，缩短分频范围，对我们提高计算效率是有帮助的。
[0062] 可控震源高效采集方式，如滑动扫描采集、Vl采集等，大大提高了勘探效率，降低 了可控震源地震勘探的成本。因此，对可控震源的工作原理、激发参数、处理技巧等方面进 行研究，提高采集、处理可控震源数据的能力对适应国内地震勘探的发展趋势是非常必要 的。为此本发明提供一种分频动态编码混叠数据成像技术，开发基于交替扫描、滑动扫描、 独立同步扫描（ISS)、远距离同步扫描（DSSS)等高效采集技术的混叠数据成像技术，为高 效采集方式采集的数据进行后续的解释工作提供成像基础。
[0063] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，应此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 基于分频动态编码的可控震源混叠数据成像方法，其特征在于，该基于分频动态编 码的可控震源混叠数据成像方法包括： 步骤1，输入震源和炮记录； 步骤2,通过波场对震源和炮记录进行波场延拓后重建编码前的震源波场和接收波 场； 步骤3,输入分频编码矩阵； 步骤4,应用编码矩阵对震源波场和接收波场进行编码； 步骤5,对编码后的震源波场和接收波场进行成像；W及 步骤6,输出成像结果。2. 根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在 于，在步骤2中，重建的震源波场和接收波场为：其中：W,、分别代表通过波场对震源和炮记录进行波场延拓后重建的震源波场 和接收波场，Ns表示炮记录的总数，Wg、WK表示物理相关的震源波场和接收波场。3. 根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在 于，在步骤3中，使用Em,。表示编码矩阵，m表示行指数，即炮指数，η是列指数，即编码次数， 把矩阵的元素表示成如下的形式： Em,n=Am,neXP(-j4m,n)做 A代表第η次编码中第m个炮集的振幅权重，dv。代表第η次编码中第m个炮集的相 移，假设编码矩阵为一个频率独立的编码矩阵，运个矩阵为Ξ维的，第Ξ维是频率ω，因此 编码矩阵可W写成： Em,n(") = Am,n(")eXP(-j Φη,η(")) (4) 分频编码的编码矩阵表示为：其中为第i个频率。4. 根据权利要求3所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在 于，在步骤3中，分频编码的方式目的是在i声j时使3 4，其中，i、j分别是 编码矩阵的纵横坐标，W：、wi分别代表通过波场对震源和炮记录进行波场延拓后重建的 震源波场和接收波场；在对一个超炮道集进行编码时，将频率域的波场进行分频，在对使每 个震源波场分别使用不同的频率段，且一个超道集中个震源使用的频率没有重叠，对应的 接收波场频率分配情况与震源波场相同。5. 根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在 于，在步骤4中，对炮集编码形成超炮道集用炮集与编码矩阵相乘来表示，那么编码后的震 源波场和接收波场为： Bs=WsE(5) Br=WrE(6) 其中化、B汾别为编码后的震源波场和接受波场，Wg、Wu分别为重建后的震源波场和接 收波场，E为编码矩阵。6.根据权利要求1所述的基于分频动态编码的可控震源混叠数据成像方法，其特征在 于，在步骤5中，对编码后的震源波场和接收波场进行成像的表达式为：(7) 其中，R。为编码后的震源波场和接收波场的成像结果，Bg为编码后的震源波场，^<是 编码后接收波场的共辆转置。
【专利摘要】本发明提供一种基于分频动态编码的可控震源混叠数据成像方法，该基于分频动态编码的可控震源混叠数据成像方法包括：步骤1，输入震源和炮记录；步骤2，通过波场对震源和炮记录进行波场延拓后重建编码前的震源波场和接收波场；步骤3，输入分频编码矩阵；步骤4，应用编码矩阵对震源波场和接收波场进行编码；步骤5，对编码后的震源波场和接收波场进行成像；以及步骤6，输出成像结果。该基于分频动态编码的可控震源混叠数据成像方法能够减少混叠数据偏移中产生的串扰噪声并且减少偏移的计算量提高效率。
【IPC分类】G01V1/28
【公开号】CN105319583
【申请号】CN201510272528
【发明人】于富文, 李振春, 魏褔吉, 胡立新, 段卫星, 付大志, 卢湘鹏, 张剑
【申请人】中石化石油工程地球物理有限公司胜利分公司
【公开日】2016年2月10日
【申请日】2015年5月26日