一种多次波成像方法及系统技术方案

本发明专利技术公开了一种多次波成像方法及系统，方法包括：对原始地震数据进行反射波与多次波分离，获取单独的反射波；提取单阶次的水层多次波与源端微屈多次波；通过离散时间分布函数与概率振幅增益函数计算水层多次波对应的正传水层多次波超道集以及源端微屈多次波对应的反传源端微屈多次波超道集；通过正传水层多次波超道集获取正向延拓的震源波场，通过反传源端微屈多次波超道集获取反向延拓的检波波场；基于改进的激发振幅成像条件计算最终成像剖面。本发明专利技术有效去除了残余同相轴串扰噪声，提升可控阶多次波成像的计算效率。提升可控阶多次波成像的计算效率。提升可控阶多次波成像的计算效率。

【技术实现步骤摘要】
一种多次波成像方法及系统


[0001]本专利技术涉及地震成像
，更具体的说是涉及一种多次波成像方法及系统。

技术介绍

[0002]地球物理勘探是油气资源开发的上游尖端技术。较其他地球物理勘探方法，地震勘探获取地下构造的成像结果信息更加丰富、精度更高，是油气资源发现、监测、增储的主要技术手段，被学术界及工业界广泛地研究并应用。
[0003]地震勘探采集到的数据充斥着大量的多次波，陆上数据以层间多次波为主，海上数据中表面多次波更为发育。多次波同样也是来自地下构造的真实反射，与一次反射波相比，多次波的传播路径、反射角度、覆盖范围等方面含有更为丰富的信息，能够传播到反射波无法照明的构造阴影区。
[0004]标准的逆时偏移成像是基于一次反射波建立的，多次波被视为数据噪声，在偏移成像之前要尽可能的消除。对含有多次波的原始数据直接进行偏移成像，多次波会在地下非构造处形成严重的串扰噪声，掩盖有效的地层，降低成像结果的信噪比，影响随后的地质解释。
[0005]多次波逆时偏移成像最大的难题是串扰噪声的产生与压制。多次波逆时偏移成像正传本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种多次波成像方法，其特征在于，包括：对原始地震数据进行反射波与多次波分离，获取单独的反射波；基于水层模型与波动方程预测得到水层反射波；基于水层反射波频域自乘积，预测0阶至N
‑
1阶水层多次波，其中，N为大于1的整数，0阶水层多次波为水层反射波；将单独的反射波分别与0阶至N
‑
1阶水层多次波进行频域乘积，预测1阶至N阶源端微屈多次波；通过离散时间分布函数与概率振幅增益函数计算0阶至N
‑
1阶水层多次波对应的正传水层多次波超道集以及1阶至N阶源端微屈多次波对应的反传源端微屈多次波超道集；通过正传水层多次波超道集获取正向延拓的震源波场，通过反传源端微屈多次波超道集获取反向延拓的检波波场；基于震源波场和检波波场并通过改进的激发振幅成像条件计算最终成像剖面。2.根据权利要求1所述的一种多次波成像方法，其特征在于，正传水层多次波超道集和反传源端微屈多次波超道集计算公式为：；其中，表示检波器位置，表示第s炮中正传水层多次波超道集，表示第s炮中反传源端微屈多次波超道集；下角标F、B分别表示偏移过程中的正传、反传；为第s炮中j阶多次波的离散时间值，遵循高斯分布；为概率振幅增益函数；N表示多次波阶数最大值；t表示时间；表示时域j阶水层多次波；表示时域j+1阶源端微屈多次波。3.根据权利要求1所述的一种多次波成像方法，其特征在于，震源波场计算公式为：；检波波场计算公式为：；其中，表示拉普拉斯算子；表示速度；表示正向延拓的震源波场；表示反向延拓的检波波场。
4.根据权利要求1所述的一种多次波成像方法，其特征在于，改进的激发振幅成像条件计算公式为：；其中，表示震源波场第l个极大值对应的传播时刻；L表示波场极大值的个数，为大于1的整数；表示总炮数；表示时刻对应的震源波场值，表示时刻对应的检波波场值，表示成像结果。5.根据权利要求1所述的一种多次波成像方法，其特征在于，单独的反射波频域计算公式为：；式中，矩阵分别代表频域原始地震数据、反射波、多次波；矩阵代表频域地表吸收算子。6.根据权利要求1所述的一种多次波成像方法，其特征在于，0阶至N
‑
1阶水层多次波和1阶至N阶源端微屈多次波具体计算过程为：输入水层速度模型，并通过波动方程获得预测水层反射波；将水层反射波进行N
‑
1次频域自乘积，获取预测的0阶至...

【专利技术属性】
技术研发人员：张延保，胡峰，赫曼求，
申请(专利权)人：浙江数智交院科技股份有限公司，
类型：发明
国别省市：