一种波阻抗谱反演方法及系统技术方案

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

【技术实现步骤摘要】
一种波阻抗谱反演方法及系统


[0001]本专利技术属于地震信号处理领域，具体涉及一种波阻抗谱反演方法及系统。

技术介绍

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

技术实现思路

[0004]本专利技术的特征和优点在下文的描述中部分地陈述，或者可从该描述显而易见，或者可通过实践本专利技术而学习。
[0005]为克服现有技术的问题，本专利技术提供可选地一种波阻抗谱反演方法，包括：
[0006]S1、获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；
[0007]S2、获取地震记录并据此得到地震记录频率域振幅对数谱；
[0008]S3、根据所述井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；
[0009]S4、将所述地震记录与算子结合，形成相对波阻抗。
[0010]可选地，所述步骤S1包括：先将井曲线通过傅立叶变换变换到频率域得到井曲线频率域振幅谱，再将井曲线频率域振幅谱取对数得到井曲线频率域振幅对数谱，随后计算井曲线频率域振幅对数谱的各频率的平均值曲线，并对平均值曲线进行拟合，得到井波阻抗频率域振幅对数谱。
[0011]可选地，所述步骤S2包括：先将地震道曲线通过傅立叶变换变换到频率域得到频率域振幅谱，再将频率域振幅谱取对数得到地震记录频率域振幅对数谱，计算地震记录频率域振幅对数谱的各频率的平均值曲线，对平均值曲线进行光滑，得到地震记录频率域振幅对数谱。
[0012]可选地，所述步骤S3包括：先对井波阻抗频率域振幅对数谱进行带通滤波，再结合地震记录频率域振幅对数谱计算算子的频率域振幅对数谱，将频率域振幅对数谱取幂得到算子频率域振幅谱后，并通过反傅立叶变换变换到时间域得到时间域算子，对所述时间域算子进行截短，得到算子。
[0013]可选地，所述步骤S4包括：将算子与地震记录进行褶积得到相对波阻抗，其后加载
波阻抗低频分量，得到绝对波阻抗。
[0014]本专利技术提供一种波阻抗谱反演系统，包括：
[0015]第一处理单元，用于获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；
[0016]第二处理单元，用于获取地震记录并据此得到地震记录频率域振幅对数谱；
[0017]算子获取单元，用于根据所述井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；
[0018]波阻抗形成单元，将所述地震记录与算子结合，形成相对波阻抗。
[0019]可选地，所述第一处理单元具体用于：先将井曲线通过傅立叶变换变换到频率域得到井曲线频率域振幅谱，再将井曲线频率域振幅谱取对数得到井曲线频率域振幅对数谱，随后计算井曲线频率域振幅对数谱的各频率的平均值曲线，并对平均值曲线进行拟合，得到井波阻抗频率域振幅对数谱。
[0020]可选地，所述算子获取单元具体用于：先对井波阻抗频率域振幅对数谱进行带通滤波，再结合地震记录频率域振幅对数谱计算算子的频率域振幅对数谱，将频率域振幅对数谱取幂得到算子频率域振幅谱后，并通过反傅立叶变换变换到时间域得到时间域算子，对所述时间域算子进行截短，得到算子。
[0021]可选地，所述波阻抗形成单元具体用于：将算子与地震记录进行褶积得到相对波阻抗，其后加载波阻抗低频分量，得到绝对波阻抗。
[0022]本专利技术提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机可执行的至少一个程序，所述至少一个程序被所述计算机执行时使所述计算机执行本专利技术任一实施例提供的方法中的步骤。
[0023]本专利技术提供的波阻抗谱反演方法及系统可以获得较为准确的波阻抗反演结果。
附图说明
[0024]图1是本专利技术实施例的波阻抗谱反演方法的流程图；
[0025]图2是本专利技术实施例的波阻抗谱反演系统的结构示意图；
[0026]图3A是处理前的频谱示意图；
[0027]图3B是处理后的频谱示意图；
[0028]图3C是波阻抗低频模型剖面(10-15Hz)；
[0029]图3D是采用本专利技术一实施例获得的绝对波阻抗反演结果；
[0030]图3E是平滑处理后的绝对波阻抗反演结果；
具体实施方式
[0031]下面结合附图对本专利技术作进一步详细描述：
[0032]波阻抗谱反演算法是一种对井的依赖程度比较小的反演算法。它实际上是用地震的频谱和井的波阻抗的频谱相匹配来完成反演，所以这种反演方法没有明显的子波提取过程，也不需要初始模型，纵向分辨率比稀疏脉冲高，由于其客观、简洁特别适合于储层非均质很强的岩性圈闭预测。
[0033]如图1所示，本专利技术提供波阻抗谱反演方法，包括以下步骤：
[0034]S10、获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；
[0035]在具体实施时，先将井曲线通过傅立叶变换变换到频率域得到井曲线频率域振幅谱，再将井曲线频率域振幅谱取对数得到井曲线频率域振幅对数谱，随后计算井曲线频率域振幅对数谱的各频率的平均值，最后通过曲线拟合，得到光滑的井波阻抗频率域振幅对数谱。
[0036]S20、获取地震记录并据此得到地震记录频率域振幅对数谱；
[0037]地震记录是由地震道组成的，将地震道用曲线绘制出来即可得到地震道曲线。在具体实施时，先将地震道曲线通过傅立叶变换变换到频率域得到频率域振幅谱，再将频率域振幅谱取对数得到地震记录频率域振幅对数谱，计算地震记录频率域振幅对数谱的各频率的平均值曲线，对平均值曲线进行光滑后，得到光滑的地震记录频率域振幅对数谱。
[0038]S30、根据所述井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；
[0039]在具体实施时，先对井波阻抗频率域振幅对数谱进行带通滤波，再结合地震记录频率域振幅对数谱计算算子的频率域振幅对数谱，将频率域振幅对数谱取幂得到算子频率域振幅谱后，并通过反傅立叶变换变换到时间域得到时间域算子，对所述时间域算子进行截短，得到算子。
[0040]S40、将所述地震记录与算子结合，形成相对波阻抗。
[0041]在具体实施时，将算子与地震记录进行褶积得到相对波阻抗，其后加载波阻抗低频分量，得到绝对波阻抗。
[0042]步骤S10中，设有M条井曲线p
m
(t),m＝1,2L M，此时，井波阻抗处理具体包本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种波阻抗谱反演方法，其特征在于，包括：S1、获取井曲线，并进行井波阻抗处理，得到井波阻抗频率域振幅对数谱；S2、获取地震记录并据此得到地震记录频率域振幅对数谱；S3、根据所述井波阻抗频率域振幅对数谱和地震记录频率域振幅对数谱获取算子；S4、将所述地震记录与算子结合，形成相对波阻抗。2.根据权利要求1所述波阻抗谱反演方法，其特征在于，所述步骤S1包括：先将井曲线通过傅立叶变换变换到频率域得到井曲线频率域振幅谱，再将井曲线频率域振幅谱取对数得到井曲线频率域振幅对数谱，随后计算井曲线频率域振幅对数谱的各频率的平均值曲线，并对平均值曲线进行拟合，得到井波阻抗频率域振幅对数谱。3.根据权利要求1所述波阻抗谱反演方法，其特征在于，所述步骤S2包括：先将地震道曲线通过傅立叶变换变换到频率域得到频率域振幅谱，再将频率域振幅谱取对数得到地震记录频率域振幅对数谱，计算地震记录频率域振幅对数谱的各频率的平均值曲线，对平均值曲线进行光滑，得到地震记录频率域振幅对数谱。4.根据权利要求1所述波阻抗谱反演方法，其特征在于，所述步骤S3包括：先对井波阻抗频率域振幅对数谱进行带通滤波，再结合地震记录频率域振幅对数谱计算算子的频率域振幅对数谱，将频率域振幅对数谱取幂得到算子频率域振幅谱后，并通过反傅立叶变换变换到时间域得到时间域算子，对所述时间域算子进行截短，得到算子。5.根据权利要求1所述波阻抗谱反演方法，其特征在于，所述步骤S4包括：将算子与地震记录进行褶积得到相对波阻抗，其后加载波阻抗低频分量，得到绝对波阻抗。6....

【专利技术属性】
技术研发人员：孟黎歌，朱凌燕，王婷婷，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：