一种地震勘探特征子波整形提高分辨率方法技术

本发明专利技术提供了地震勘探特征子波整形提高分辨率方法，该方法提出先进行子波处理提升高频成份的子波的振幅，用于提升子波能量；再用复数Paul子波(PW)替换原来的子波进行分频重构的处理思路，旨在最大程度提高分辨率和保持构造形态特征。在实际勘探实践工作中，应用本发明专利技术提供的方法后，不论地震勘探特征子波整形分辨率还是信噪比均有很大程度的提高。

【技术实现步骤摘要】

本专利技术涉及地震勘探
，尤其涉及一种地震勘探特征子波整形提高分辨率方法。
技术介绍
随着人类油气勘探活动的不断深入，物探、地质等多学科也在逐步交叉、融合发展，这些学科不仅是人们实践经验的总结，很多都已上升为勘探理论。地震资料处理是地震勘探工作中的重要一环，且地震资料处理是一个复杂的动态过程，地震资料处理需要提高叠前、叠后地震资料的分辨率，在此基础上进行储层预测、反演及含气性识别。在高分辨率资料处理中，地震资料的分辨率依赖于每个处理环节，因此，地震资料分辨率的提高与保持，是一个全方位的工作。有效地拓宽地震记录的频带和提高原有的分辨率是高分辨率资料处理中的两个重要环节。针对地震资料及数据处理过程中多个环节的仔细分析研究，发现通过对处理过程中会受到处理技术较大影响并导致分辨率无法提高的参数的研究，可以得出新的有价值的提高分辨率的方法。
技术实现思路
针对现有技术展宽频带、提升分辨率的需要和对地震数据处理过程的深入研究，本专利技术提供一种新的提高分辨率的方法。实际资料处理结果表明，应用该方法后，不论分辨率还是信噪比均有很大程度的提高。本专利技术采用的技术方案：一种地震勘探特征子波整形提高分辨率方法，该方法包含以下步骤：先进行子波处理提升高频成分的子波的振幅，再用复数Paul子波替换原来的子波。所述提升高频成分的子波的振幅包含以下步骤：(1)根据实际地震数据计算振幅谱；(2)计算高频振幅补偿曲线；(3)对数据进行小波变换得到小子波；(4)对所述小子波进行振幅补偿，以期望的高频振幅补偿曲线为基准，将每一小子波的最大振幅与所述期望的高频振幅补偿曲线对应值对应，求出该小子波的振幅补偿因子，然后计算每一小子波补偿后的振幅谱。所述步骤(1)中的数据是指在选择的分析时窗内，取多道数据：所述振幅谱是多道的傅氏振幅谱；所述步骤(2)的高频振幅补偿曲线计算方式是，首先根据所述傅氏振幅谱求得地层等效吸收衰减曲线，然后根据所述地层等效吸收衰减曲线计算高频振幅补偿曲线。高频振幅补偿曲线的数值可调。高频振幅补偿曲线的数值可调降为原来的2/3或1/2。所述步骤(4)中的计算每一小子波补偿后的振幅谱是通过求得振幅补偿因子来计算的。所述振幅补偿因子是通过如下步骤获得：1)进行多道数据的小波变换，得到不同频率的小波振幅谱；2)在所述小波振幅谱中选取最大主振幅值(Amax(f0))；3)在fi大于f0的范围内，求取其他小子波的最大主振幅Amax(fi)；4)每一个小子波的振幅补偿因子k(fi)为，其中，i＝1，2，3，……，0＜α＜1，fi＞f0；根据上式可计算高频振幅补偿曲线，α值即代表了调降该高频振幅补偿曲线，预设一条期望的高频振幅补偿曲线，将所述每一小波的最大振幅拟合到所述期望高频振幅补偿曲线上，从而求得上述公式中的α值，得到该条高频振幅补偿曲线针对频率fi的该小子波的振幅补偿因子k(fi)；5)在单道数据上，在选择的时窗内，进行小波变换，得到不同频率fi的小子波振幅谱A(fi)；针对频率fi的小子波，用振幅补偿因子k(fi)与该小子波中的所有振幅相乘，得到补偿后小子波的频率为fi的振幅，A'(fi)＝k×A(fi)，(k＝k(fi))对每一个频率fi的振幅作如上处理，得到所述小子波补偿后的振幅谱；6)对每一个小子波做与上述步骤相同的处理，得到每一个小子波的振幅谱。将这些新小子波重构，得到高频振幅补偿后的反射波。对所述重构的反射波进行随机噪声压缩处理。所述用复数Paul子波替换原来的子波，采用如下的公式进行复数Paul子波(即PW)分频重构：地震反射信号s表示为地震子波r和反射系数b序列的褶积，如下算式：本专利技术技术效果：本专利技术与传统提高地震勘探特征子波整形分辨率方法相比具有以下优点：该方法提出先进行子波处理提升高频成份的子波的振幅，再用复数Paul子波(PW)替换原来的子波进行分频重构的处理思路，既能够最大程度提高分辨率，又兼顾保持构造形态特征。在实际勘探实践工作中，应用本专利技术提供的方法后，不论地震勘探特征子波整形分辨率还是信噪比均有很大程度的提高。专利技术人发现，由于动静校正存在一定的误差，可使高频成份在叠加过程中相互抵消使其振幅减弱；更为严重的是，在叠前偏移处理时，为保证偏移处理的运行效率，而使用较大的网格，为此不得不在偏移处理过程中做高截滤波处理，使高频成份的振幅大幅度减小。对偏移处理数据的高通滤波扫描的测试，发现高频成份具有很高的信噪比，因而是有效成分，如果充分利用这些被过滤掉的有效的高频成份，提升高频成份的振幅，从而提升主频-展宽频带，就能达到提高分辨率的目的，根据前面的分析，在提升高频成分的振幅后，用复数Paul子波(PW)替换原来的子波进行分频重构，能有效地提高分辨率和信噪比。在先进行子波处理提升高频成份的子波的振幅的步骤中，本专利技术根据实际地震数据的地层等效吸收衰减曲线，计算高频振幅补偿曲线(参见图3中的点状虚线)，同时该曲线有很大的调整空间，可以根据实际数据的需要经过后面步骤的计算后再来进行调整以使子波整形振幅谱补偿达到期望；本专利技术根据实际地震数据的地层等效吸收衰减曲线，是本专利技术中提升高频成分的子波的振幅的一个重要依据，为实际工作提供了理想的评价标准。在地震信号的多子波分解与重构中，子波类型与子波分辨率的选择是十分重要的。针对这一问题，首先计算了复数Paul子波和复数Morlet子波的时频谱，并讨论了频率域内子波宽度控制参数对频率域分辨率的影响；然后应用匹配追踪算法对实际地震数据进行多子波分解与重构，从理论上确定了控制复数Paul子波和复数Morlet子波频率域分辨率的最佳子波宽度，对两种子波分解重构地震剖面的差异性进行了对比分析。理论和实际计算结果表明：复数Paul子波与复数Morlet子波的分辨率和频谱能量分布有所不同，在实际应用中需针对不同的研究目的选择不同类型的子波。基于复数Paul子波的重构剖面保留了更多的低频能量，适用于反射同相轴连续的地下构造研究；基于复数Morlet子波的重构剖面高、低频频谱差异大，更适用于地震波吸收衰减研究。因此本专利技术在提升高频成分的振幅后，用主频高-频带宽的复数Paul子波(PW)替换原来的子波进行分频重构，能有效地提高分辨率和信噪比。附图说明图1是实际地震数据的振幅谱示意图。图2是根据实际地震数据的振幅谱求地层等效吸收衰减曲线示意图。图3是根据实际地震数据的地层等效吸收衰减曲线，计算高频振幅补偿曲线。图4是调降高频振幅补偿曲线的数值示意图。图5是再次调降高频振幅补偿曲线的数值示意图。图6是实际地震数据的振幅谱高频补偿之后的期望振幅谱。图7是实际地震数据的傅氏振幅谱及其小波振幅谱。图8是实际地震数据的振幅谱高频补偿之后的期望振幅谱小波变换示意图。图9是实际地震数据原始剖面。图10是原始地震数据采用子波振幅谱调整后剖面(第一阶段)。图11是原始地震数据采用子波振幅谱调整后特征子波分频高频重建剖面示意图。(第二阶段)图12是原始数据频谱示意图。图13是实际地震数据采用子波振幅调整后(第一阶段)后频谱示意图。图14是实际地震数据采用Paul特征子波分频重构后(第二阶段)后频谱示意图。具体实施方式本专利技术涉及一种地震勘探特征子波整形提高分辨率方法。实践表明，有效地拓宽地震记录的频带和提高原有本文档来自技高网...

【技术保护点】
一种地震勘探特征子波整形提高分辨率方法，该方法包含以下步骤：先进行子波处理提升高频成分的子波的振幅，再用复数Paul子波替换原来的子波。

【技术特征摘要】
1.一种地震勘探特征子波整形提高分辨率方法，该方法包含以下步骤：先进行子波处理提升高频成分的子波的振幅，再用复数Paul子波替换原来的子波。2.根据权利要求1所述的一种地震勘探特征子波整形提高分辨率方法，所述提升高频成分的子波的振幅包含以下步骤：(1)根据实际地震数据计算振幅谱；(2)计算高频振幅补偿曲线；(3)对数据进行小波变换得到小子波；(4)对所述小子波进行振幅补偿，以期望的高频振幅补偿曲线为基准，将每一小子波的最大振幅与到所述期望的高频振幅补偿曲线对应值对应，求出该小子波的振幅补偿因子，然后计算每一小子波补偿后的振幅谱。3.根据权利要求2所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于所述步骤(1)中的数据是指在选择的分析时窗内，取多道数据：所述振幅谱是多道的傅氏振幅谱；所述步骤(2)的高频振幅补偿曲线计算方式是，首先根据所述傅氏振幅谱求得地层等效吸收衰减曲线，然后根据所述地层等效吸收衰减曲线计算高频振幅补偿曲线。4.根据权利要求3所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于高频振幅补偿曲线的数值可调。5.根据权利要求4所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于高频振幅补偿曲线的数值可调降为原来的2/3或1/2。6.根据权利要求3所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于所述步骤(4)中的计算每一小子波补偿后的振幅谱是通过求得振幅补偿因子来计算的。7.根据权利要求6所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于所述振幅补偿因子是通过如下步骤获得：1)进行多道数据的小波变换，得到不同频率的小波振幅谱；2)在所述小波振幅谱中选取最大主振幅值(...

【专利技术属性】
技术研发人员：唐文榜，李治国，孔选林，李飞，樊佳芳，
申请(专利权)人：北京天工新源科技有限公司，
类型：发明
国别省市：北京;11