[发明专利]一种波阻抗谱反演方法及系统

说明书

本发明提供了一种波阻抗谱反演方法及系统，该方法包括：获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；获取地震记录并据此得到地震记录频率域振幅对数谱；根据井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；将地震记录与算子结合，形成相对波阻抗。该系统包括：第一处理单元，用于获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；第二处理单元，用于获取地震记录并据此得到地震记录频率域振幅对数谱；算子获取单元，用于根据井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；波阻抗形成单元，将地震记录与算子结合，形成相对波阻抗。本发明可以获得较为准确的波阻抗反演结果。

技术领域

本发明属于地震信号处理领域，具体涉及一种波阻抗谱反演方法及系统。

背景技术

地震解释通常需要对地下地层的弹性、物性信息进行定性甚至定量的获取，以便达到更为精细的储层刻画描述。常规反演在井比较少的工区多用稀疏脉冲反演。稀疏脉冲反演非常依赖于子波的提取，而要达到全区的优化，子波必须是时空变，目前精确的子波提取和子波的时空变都是尚未完全解决的难题。一般的方法是通过多井分析，得到一个全区的综合子波，但这会影响到储层预测的精度和反演的分辨率。而地震解释人员通常希望获得的地震资料不仅仅限于高分辨率地震剖面，相应直观准确的波阻抗数据体可作为地层属性用于定量地震解释。

如何通过接收到的地震数据获取波阻抗信息也一直在为了克服各种问题而持续得到深入的研究。

发明内容

本发明的特征和优点在下文的描述中部分地陈述，或者可从该描述显而易见，或者可通过实践本发明而学习。

为克服现有技术的问题，本发明提供可选地一种波阻抗谱反演方法，包括：

S1、获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；

S2、获取地震记录并据此得到地震记录频率域振幅对数谱；

S3、根据所述井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；

S4、将所述地震记录与算子结合，形成相对波阻抗。

可选地，所述步骤S1包括：先将井曲线通过傅立叶变换变换到频率域得到井曲线频率域振幅谱，再将井曲线频率域振幅谱取对数得到井曲线频率域振幅对数谱，随后计算井曲线频率域振幅对数谱的各频率的平均值曲线，并对平均值曲线进行拟合，得到井波阻抗频率域振幅对数谱。

可选地，所述步骤S2包括：先将地震道曲线通过傅立叶变换变换到频率域得到频率域振幅谱，再将频率域振幅谱取对数得到地震记录频率域振幅对数谱，计算地震记录频率域振幅对数谱的各频率的平均值曲线，对平均值曲线进行光滑，得到地震记录频率域振幅对数谱。

可选地，所述步骤S3包括：先对井波阻抗频率域振幅对数谱进行带通滤波，再结合地震记录频率域振幅对数谱计算算子的频率域振幅对数谱，将频率域振幅对数谱取幂得到算子频率域振幅谱后，并通过反傅立叶变换变换到时间域得到时间域算子，对所述时间域算子进行截短，得到算子。

可选地，所述步骤S4包括：将算子与地震记录进行褶积得到相对波阻抗，其后加载波阻抗低频分量，得到绝对波阻抗。

本发明提供一种波阻抗谱反演系统，包括：

第一处理单元，用于获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；

第二处理单元，用于获取地震记录并据此得到地震记录频率域振幅对数谱；

算子获取单元，用于根据所述井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；

波阻抗形成单元，将所述地震记录与算子结合，形成相对波阻抗。