[发明专利]一种起伏地表的地下构造成像方法

摘要

本发明公开了一种起伏地表的地下构造成像方法，包括以下步骤：S1，输入观测数据，偏移速度场，同时预先设置反射系数模型的初始值；S2，对观测数据进行基于相关系数的采样或重新采样，得到偏移数据；S3，根据当前的反射系数模型和偏移数据获得数据残差；S4，根据数据残差，确定梯度，基于梯度进行反射系数模型迭代，获得最终的反射系数模型。本发明方法通过引入基于相关系数的采样策略，避免了随机采样，在减小偏移数据量的同时能有效地避免偏移数据随迭代次数产生巨大变化，进而减弱更换偏移数据对梯度方向共轭性的破坏及收敛曲线的抖动现象，提高了迭代时的计算效率和最终的成像质量。

主权项

1.一种起伏地表的地下构造成像方法，其特征在于，包括以下步骤：S1，输入观测数据，偏移速度场，同时预先设置反射系数模型的初始值；S2，对观测数据进行基于相关系数的采样或重新采样，得到偏移数据，具体包括：当迭代次数k为1时，按照固定的采样间隔，选择Ns个沿偏移距方向均匀分布的观测数据作为偏移数据，且采样间隔大于炮间隔；当迭代次数k大于1时，判断迭代次数k除以模型更新方向控制因子o的余数是否为1，若是，对观测数据进行重新采样，否则，不重新采样，得到偏移数据，其中重新采样的具体方法包括：计算剩余数据中每炮数据与上一次迭代时偏移数据的相关系数，选取相关系数最大的Ns个单炮数据作为本次迭代的偏移数据；当剩余数据中单炮数据的数量不足Ns时，不足的部分从已计算过的数据中随机选取，所述相关系数的其计算方法为：其中，d_j(x，t)和为上次迭代时偏移数据中的单炮数据及其均值，d_i(x，t)和为剩余数据中的单炮数据及其均值，N_s为偏移数据中单炮数据的数量，n_x为单炮数据中空间采样点数，n_t为时间采样点数，i，j代表炮号；S3，根据当前的反射系数模型和偏移数据获得数据残差；S4，根据数据残差，确定随机梯度和/或随机共轭梯度，基于随机梯度和/或随机共轭梯度进行反射系数模型迭代，获得最终的反射系数模型。