[发明专利]多次波发育区地震数据保真速度分析方法

权利要求书

1.一种多次波发育区地震数据保真速度分析方法，其特征在于，包括：

步骤S1、将采集的原始地震数据进行预处理，通过抽道集处理得到叠前共中心点CMP道集；

步骤S2、对步骤S1得到的CMP道集进行随机噪音压制处理，得到高信噪比CMP道集；

步骤S3、对步骤S2得到的高信噪比CMP道集进行粗化网格定义，速度分析网格密度为a*b个面元，采用全频带数据生成用于速度分析的速度谱和超道集；其中，a、b为正整数；

步骤S4、对步骤S3生成的用于速度分析道集进行FK滤波压制相干噪音，之后进行空间预测误差滤波进一步压制随机噪音；

步骤S5、在步骤S4生成的道集中，通过测井曲线、区域地质认识和地震数据综合分析确定该区多次波成因及类型，并根据该区多次波成因及类型选择与其匹配的压制处理方式；

步骤S6、对步骤S5多次波压制后的道集进行时频分析，保留一次反射波优势频带数据进行滤波处理，保留最大偏移距70%偏移距范围内的数据，形成用于速度分析的优势频带/优势偏移距超道集数据；

步骤S7、在步骤S6生成的速度谱和超道集上进行速度分析，得到步骤3定义的网格基础上的一轮速度分析结果；

步骤S8、利用步骤7得到的一轮偏移速度进行与速度分析网格对应的目标线叠前时间偏移，将偏移结果抽CRP道集形成叠加剖面，质控本轮速度分析结果精度；

步骤S9、将步骤8得到的CRP道集进行反动校正，重复3-5轮步骤S3-步骤S8的处理流程，通过3-5轮迭代处理速度分析网格需要加密到a/2*b/2- a/4*b/4个面元，并压制上一轮迭代中道集中残余的多次波；

步骤S10、输出保真偏移速度场： 经过步骤9循环3-5轮迭代后新一轮偏移结果无改善、剩余时差不存在的时候，终止迭代，输出保真偏移速度场。

2.如权利要求1所述的多次波发育区地震数据保真速度分析方法，其特征在于，步骤S1中，所述地震数据准备及预处理包括：将野外采集的原始地震单炮记录通过解编、置观测系统、静校正、去噪、几何扩散补偿与地表一致性振幅补偿，反褶积及反褶积后去噪处理。

3.如权利要求1所述的多次波发育区地震数据保真速度分析方法，其特征在于，所述步骤S5中根据该区多次波成因及类型选择与其匹配的压制处理方式具体包括：若对象工区多次波以自由表面多次波为主，则采用波动方程理论的波场外推SRME方法压制自由表面多次波；若对象工区多次波以层间多次波为主，采用加权最小二乘拉冬变换压制层间多次波。

4.如权利要求1所述的多次波发育区地震数据保真速度分析方法，其特征在于，包括：步骤S9中，重复5轮步骤S3-步骤S8的处理流程，通过5轮处理将速度分析网格加密到a/4*b/4个面元，并压制上一轮迭代中道集中残余的多次波具体包括：当前本轮处理将速度分析网格加密到a/2*b或者a*b/2个面元，并压制上一轮迭代中道集中残余的多次波；当完成第二轮步骤S3-步骤S8的处理后得到第二轮的偏移速度，对第二轮迭代得到的偏移速度对应的偏移结果进行质控，若CRP 道集中仍存在剩余时差，将进行第三步轮重复步骤S3-步骤S8迭代；第三轮处理将速度分析网格需要加密到a/2*b/2个面元；

当完成第三轮的步骤S3-步骤S8处理后得到第三轮的偏移速度，对第二轮迭代得到的偏移速度对应的偏移结果进行质控，若CRP 道集中仍存在剩余时差，将进行下一轮重复步骤S3-步骤S8迭代，进行第四轮速度分析，则第四轮处理将速度分析网格需要加密到a/4*b/2或a/2*b/4个面元，进行第五轮速度分析迭代，则第五轮处理将速度分析网格需要加密到a/4*b/4个面元，其中，a、b为正整数。