杨育臣，方金伟，王宁，李松龄

1 Geology and Geophysics Program，Missouri University of Science and Technology，Rolla，MO 65409，USA 2 深部岩土力学与地下工程国家重点实验室，中国矿业大学力学与土木工程学院，徐州 221116 3 东北石油大学地球科学学院，大庆 163318

0 引言

近年来，发展了宽频带、宽方位、高密度的地震数据采集技术，通过地震成像技术实现地下参数的高信噪比、高分辨率和高保真度的建模.其中，具有最高成像精度的反演成像以全波形反演(FWI)为代表，实现构造信息和岩性油气藏的高精度参数建模.全波形反演通过建立观测和合成数据的某种距离度量的目标泛函，利用优化理论计算出参数下降方向和更新步长，不断修改模型直至合成数据接近观测数据，从而获得准确的模型参数.全波形反演能够充分利用地震数据中的振幅和相位信息，获得高分辨率高精度的速度模型.

计算效率低下是制约全波形反演进行大规模计算以及实际资料处理的主要因素，因此提升全波形反演的效率就显得十分重要(邢贞贞等，2019).除了使用高性能计算平台加速之外(Yang et al.，2015；Gokhberg and Fichtner，2016；张猛等，2014；桂生等，2017)，通常使用算法类的效率优化方案来加速反演.基于算法类的加速方法通常使用震源编码策略提升全波形反演效.震源编码的思想是通过某种震源组合方式将很多单炮形成一个超级震源(编码炮)，实现几十炮甚至上百炮的同时模拟，从而减少全波形反演中波场模拟的次数，极大地减少全波形反演的计算量.然而，相比于传统的同时源模拟(震源编码)方法，比如随机源编码(Krebs et al.，2009)、随机时移(Dai and Schuster，2009)、平面波编码(Zhang et al.，2003；……