提高地震资料分辨率效果评判优选方法技术

本发明专利技术公开了一种提高地震资料分辨率效果评判优选方法，步骤1，导入不同方法或参数的一组三维地震勘探单炮数据、道集数据和叠后剖面数据，确定目的层范围；步骤2，通过数学公式计算目的层数据的振幅、频率、相位、信噪比、时频特征、相干特征、振幅比特征、波形特征、断层特征等属性信息；步骤3，在属性计算公式基础上，利用分类混合交汇公式对地震资料提高分辨率效果进行评判；步骤4，根据计算得到的分类混合交汇效果评判参数，实现提高分辨率结果的评判优选分析，得出最佳的提高分辨率处理结果。该方法能够高效地实现提高分辨率处理结果的效果评判，优选出最理想的提高分辨率处理结果，从而为地震资料提高分辨率处理提供指导和全过程质量监控。

Evaluation and optimization method for improving resolution effect of seismic data

The invention discloses a method for evaluating and optimizing the effect of improving the resolution of seismic data. Step 1 introduces a group of three-dimensional seismic exploration single shot data, gather data and post-stack profile data with different methods or parameters to determine the target layer range; Step 2 calculates the amplitude, frequency, phase, signal-to-noise ratio, time-frequency characteristics, coherence characteristics and amplitude ratio of target layer data through mathematical formulas. Attribute information such as feature, waveform feature and fault feature; Step 3, based on attribute calculation formula, the effect of improving resolution of seismic data is evaluated by using classified mixed intersection formula; Step 4, according to the evaluation parameters of classified mixed intersection effect calculated, the evaluation and optimization analysis of improving resolution result is realized, and the best processing result of improving resolution is obtained. This method can efficiently evaluate the effect of improving resolution processing results, and optimize the most ideal result of improving resolution processing, so as to provide guidance and whole process quality monitoring for improving resolution processing of seismic data.

【技术实现步骤摘要】
提高地震资料分辨率效果评判优选方法
本专利技术涉及油气田地震勘探与开发
，特别是涉及到一种提高地震资料分辨率效果评判优选方法。
技术介绍
随着地震勘探与开发程度的不断深入，特别是薄互层油气藏、岩性油气藏、隐蔽油气藏等，对地震资料的分辨率要求越来越高，地震资料提高分辨率处理是地震资料处理的重要环节之一，提高分辨率的处理方法有很多，主要包括：单道预测反褶积、地表一致性预测反褶积、道统计反褶积、盲反褶积、零相位反褶积、混合相位反褶积、子波整形反褶积、井控反褶积、反Q滤波、谱模拟反褶积、时频域反褶积、谱白化等方法。各种反褶积方法各有优缺点，并对地震资料特征具有不同的适应条件，在进行地震资料提高分辨率处理时，需要进行大量的提高分辨率方法和参数试验，得到了大量的提高分辨率处理结果，而且所得到的提高分辨率处理效果的好坏难以准确衡量与评判。特别是随着高精度、高密度三维地震勘探的不断发展，高精度、高密度三维地震资料数据量巨大，对地震数据进行提高分辨率处理后，可以得到大量的处理结果，但是无法确定哪一种提高分辨率处理结果是最理想的，在进行方法和参数试验时，对处理结果进行分析占用了大量的处理时间，制约着地震资料的高效快速处理，拖长了整个处理周期。
技术实现思路
本专利技术的目的在于针对现有技术的不足，提供一种高效实现提高分辨率处理结果的效果评判及最理想结果优选的提高地震资料分辨率效果评判优选方法。本专利技术则通过数学公式计算目的层数据的振幅、频率、相位、信噪比、时频特征、相干特征、振幅比特征、波形特征、断层特征等属性信息，并提出了时频谱互相关、频率变化、时间变化、断层特征等属性信息，并赋予每一类信息一个加权系数，综合考虑所有属性信息和相对应的加权系数，根据计算得到的分类混合交汇参数对提高分辨率处理结果进行效果评判和优选，从大量结果中优选最理想效果，最终可以快速高效的得到最优提高分辨率处理结果，缩短了地震资料处理周期，从而达到高效勘探与开发的目的。本专利技术采用的技术方案如下。提高地震资料分辨率效果评判优选方法，包括如下步骤：步骤1：导入不同提高分辨率方法或参数的一组提高分辨率处理后三维地震勘探单炮数据、道集数据和叠后剖面数据，并确定目的层范围，或者导入拾取的层位信息；步骤2：通过数学公式计算目的层数据的振幅、频率、相位、信噪比、时频特征、相干特征、振幅比特征、波形特征、断层特征等属性信息，并可根据地震资料的特征，对属性信息进行加减乘除的混合组合，得到反映提高分辨率效果的属性信息集合体；步骤3：在属性计算公式的基础上，对所有属性信息进行分类，并赋予每一类属性信息一个加权系数，频率类加权系数为α，信噪比类加权系数为β，频率变化类加权系数为φ，时间变化类加权系数为时频互相关类加权系数为ε，利用分类混合交汇公式对地震资料提高分辨率效果进行评判，综合考虑所有属性信息和相对应的加权系数，分类混合交汇公式为：其中：αi为频率类加权系数，Mi为频率类属性个数,fi(x)为频率类相关的属性信息，βj为信噪比类加权系数，Nj为信噪比类属性个数,sj(x)为信噪比类相关的属性信息，φk为频率变化类加权系数，Kk为频率变化类属性个数,bk(x)为频率变化类相关的属性信息，为时间变化类加权系数，Pp为时间变化类属性个数,tp(x)为时间变化类相关的属性信息，εq为时频互相关类加权系数，Qq为时频互相关类属性个数,cq(x)为互相关类相关的属性信息；步骤4：根据计算得到的分类混合交汇效果评判参数Dω，从而高效准确的实现提高分辨率处理结果的效果评判，优选得到最佳提高分辨率处理结果，为地震资料提高分辨率处理提供指导和全过程质量监控。作为优选技术方案，在步骤1中，单炮数据、道集数据和叠后剖面数据中的任意一种数据即可，或者三种数据同步输入，导入数据的种类越多，分析结果越准确可靠。并输入相对应数据的层位信息或目的层范围，从而确定效果评价范围，所有输入数据必须具有相同的采样点数、采样间隔和地震道数，并具有相同的数据格式、读写方式、增益水平和预处理流程。作为优选技术方案，在步骤2中，通过地震资料波形互相关公式和波形自相关公式计算稳定信噪比参数SNR＝Cor/Cor和波形相似参数其中：h0(t)为提高分辨率处理前地震波形数据，h1(t)和h2(t)为不同方法或参数的提高分辨率处理后地震波形数据，t为时间，单位ms，Cor为波形互相关，Corself为波形自相关。作为优选技术方案，在步骤2中，通过傅里叶变换公式计算地震数据的主频、低频段、高频段、有效频带、优势频带等属性信息，傅里叶变换公式为根据傅里叶变换振幅谱可以计算得到相应的属性信息，主要包括：主频参数为fmain＝maxf[H1(ω)]，低频段参数为flow＝maxf[H1(ω)]×low，高频段参数为fhigh＝maxf[H1(ω)]×high，有效频带宽度参数为feff＝maxf[H1(ω)]×0.8high-maxf[H1(ω)]×0.8low，优势频带宽度参数为fdom＝maxf[H1(ω)]×high-maxf[H1(ω)]×low。其中：h1(t)为提高分辨率处理后地震道数据，t为时间，单位ms；H1(ω)为傅里叶变换振幅谱，ω为频率，单位Hz；maxf[]为振幅谱最大值所对应的频率值，low为主频位置处振幅值的低频百分比，high为主频位置处振幅值的高频百分比，有效频带相对于优势频带具有一个占比系数0.8。作为优选技术方案，在步骤2中，通过修正S变换计算信号占比、主频变化、频宽变化、频谱相似系数等属性信息，修正S变换计算公式为：信号占比参数为主频变化参数为fmain(t)＝maxf[GNST1(t,f)]，频宽变化参数为fband(t)＝maxhigh[GNST1(t,f)]-maxlow[GNST1(t，f)]，二维时频谱相似参数为其中：h1(t)为提高分辨率处理后地震道数据，t为时间，单位ms；σ和ρ为修正S变换的分辨率参数，f为频率值；GNST1(τ,f)为提高分辨率处理后地震数据的二维时频谱，f1为信号低频下限，f2为信号高频上限，fN为地震数据采样间隔对应的最大频率值，单位Hz，maxf[]为时频谱中每个时刻最大值所对应的频率值，maxhigh[]为时频谱中每个时刻所对应的高频值，maxlow[]为时频谱中每个时刻所对应的低频值。作为优选技术方案，在步骤3中，通过数学公式计算得到属性信息后，对属性信息进行归类，并赋予加权系数，频率加权系数为α，信噪比加权系数为β，频率变化加权系数为φ，时间变化加权系数为时频互相关类加权系数为ε。针对不同的地震资料数据，加权系数需要进行调整和优化，加权系数的取值范围为[1,100]，并且其中：αi为频率类加权系数，Mi为频率类属性的个数，βj为信噪比类加权系数，Nj为信噪比类属性的个数,φk为频率变化类加权系数，Kk为频率变化类属性的个数,为时间变化类加权系数，Pp为时间变化类属性的个数,εq为时频互相关类加权系数，Qq为时频互相关类属性的个数。作为优选技术方案，在步骤4中，将所有输入的地震数据和提高分辨率处理前的数据相互进行分类混合交汇效果评判公式的计算，得出每个数据的分类混合交汇效果评判参数Dω。该评判参数Dω的大小表示地震数据提高分辨率处理效果的好坏，评本文档来自技高网...

【技术保护点】
1.提高地震资料分辨率效果评判优选方法，其特征是包括：步骤1，导入提高分辨率处理后三维地震勘探的单炮数据、道集数据和叠后剖面数据的一种或几种，输入相对应数据的层位信息或目的层范围；步骤2，通过数学公式计算目的层数据的属性信息，包括振幅、频率、相位、信噪比、时频特征、相干特征、振幅比特征、波形特征和断层特征，并根据地震资料的特征，对属性信息进行加减乘除的混合组合，得到反映提高分辨率效果的属性信息集合体；步骤3，在属性计算公式的基础上，对所有属性信息进行分类，并赋予每一类属性信息一个加权系数，频率类加权系数为α，信噪比类加权系数为β，频率变化类加权系数为φ，时间变化类加权系数为

【技术特征摘要】
1.提高地震资料分辨率效果评判优选方法，其特征是包括：步骤1，导入提高分辨率处理后三维地震勘探的单炮数据、道集数据和叠后剖面数据的一种或几种，输入相对应数据的层位信息或目的层范围；步骤2，通过数学公式计算目的层数据的属性信息，包括振幅、频率、相位、信噪比、时频特征、相干特征、振幅比特征、波形特征和断层特征，并根据地震资料的特征，对属性信息进行加减乘除的混合组合，得到反映提高分辨率效果的属性信息集合体；步骤3，在属性计算公式的基础上，对所有属性信息进行分类，并赋予每一类属性信息一个加权系数，频率类加权系数为α，信噪比类加权系数为β，频率变化类加权系数为φ，时间变化类加权系数为时频互相关类加权系数为ε；利用分类混合交汇公式对地震资料提高分辨率效果进行评判，综合考虑所有属性信息和相对应的加权系数，分类混合交汇公式为：其中：αi为频率类加权系数，Mi为频率类属性个数,fi(x)为频率类相关的属性信息，βj为信噪比类加权系数，Nj为信噪比类属性个数,sj(x)为信噪比类相关的属性信息，φk为频率变化类加权系数，Kk为频率变化类属性个数,bk(x)为频率变化类相关的属性信息，为时间变化类加权系数，Pp为时间变化类属性个数,tp(x)为时间变化类相关的属性信息，εq为时频互相关类加权系数，Qq为时频互相关类属性个数,cq(x)为互相关类相关的属性信息；步骤4，根据计算得到的分类混合交汇效果评判参数Dω，实现提高分辨率处理结果的效果评判，优选得到最佳提高分辨率处理结果。2.根据权利要求1所述的提高地震资料分辨率效果评判优选方法，其特征是：在步骤1中，同步导入所要分析的单炮数据、道集数据和叠后剖面数据，导入数据须具有相同的采样点数、采样间隔和地震道数，并具有相同的数据格式、读写方式、增益水平和预处理流程。3.根据权利要求2所述的提高地震资料分辨率效果评判优选方法，其特征是：在步骤2中，通过地震资料波形互相关公式和波形自相关公式计算稳定信噪比参数SNR＝Cor/Corself和波形相似参数其中：h0(t)为提高分辨率处理前地震波形数据，h1(t)和h2(t)为不同方法或参数的提高分辨率处理后地震波形数据，t为时间，单位ms，Cor为波形互相关，Corself为波形自相关。4.根据权利要求3所述的提高地震资料分辨率效果评判优选方法，其特征是：在步骤2中，通过傅里叶变换公式计算地震数据的属性信息包括主频、低频段、高频段、有效频带和优势频带，傅里叶变换公式为根据傅里叶变换振幅谱计算得到相应的属性信息，包括：主频参数为...

【专利技术属性】
技术研发人员：刁瑞，刘浩杰，韩站一，芮拥军，冯玉苹，关键，王蓬，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：北京,11