一种地震数据分形保幅方法
【专利摘要】本发明涉及一种地震数据分形保幅方法，其包括以下步骤：1）针对拓频前的地震数据体X0和拓频后的地震数据体Y，采用带通算子Bi分别相应提取它们的局部频段信息，记为BiX0和BiY，i=1,2,3…M，不同的i代表不同的带通频段；2）基于地震信号局部频段成分与整体的特征存在的相似性，在每一对数据体BiX0和BiY之间分别建立保幅关系，并相应提取局部保幅因子Φi；3）对所述步骤2）提取的多个局部保幅因子Φi进行平滑和组合处理，得到整体保幅因子Φ；4）采用整体保幅因子Φ作用于数据体Y，对全频段信号Y实现误差精度允许范围内的保幅。本发明可以广泛应用于对拓频后的地震数据进行振幅保真处理，进一步为后续地震资料岩性反演和高精度解释提供高保真数据。
【专利说明】一种地震数据分形保幅方法
【技术领域】
[0001]本发明涉及油气地震勘探领域中叠前及叠后地震数据高保真方法，特别是一种关于对拓频后的地震数据进行振幅保真处理，进而为后续地震资料岩性反演和高精度解释提供高保真数据的地震数据分形保幅方法。
【背景技术】
[0002]地震数据振幅保真是地震勘探领域处理的重点和难点，恢复地震数据的真实振幅对于岩性成像、储层预测等尤其重要。保幅通常分为两个概念:1、对数据采集和前期处理过程中因地层吸收衰减、去噪技术等因素造成的振幅变化进行振幅恢复；2、地震数据经过高分辨率处理后，展宽频带的保幅处理。
[0003]现有技术中针对第二种情况关于地震数据拓频后的振幅保真问题，传统的做法是将高分辨率处理后的宽频信号能量保幅到处理前的窄频能量水平。事实上，这种做法一定程度上破坏了地震数据的真实振幅，包括原有频段的振幅信息，相位信息也存在一定的损伤。
【发明内容】
[0004]针对上述问题，本发明的目的是提供一种地震数据分形保幅方法，能够更好地恢复全频段振幅和相位信息，为后续地震资料处理和解释提供高品质资料。
[0005]为实现上述目的，本发明采取以下技术方案:一种地震数据分形保幅方法，其包括以下步骤:1)针对拓频前的地震数据体Xtl和拓频后的地震数据体Y，采用带通算子Bi分别相应提取它们的局部频段信息，记为BiXc^P BiY, i=l, 2，3…M，不同的i代表不同的带通频段；2)基于地震信号局部频段成分与整体的特征存在的相似性，在每一对数据体BiXtl和BiY之间分别建立保幅关系，并相应提取局部保幅因子Oi ;3)对所述步骤2)提取的多个局部保幅因子Oi进行平滑和组合处理，得到整体保幅因子Φ:
[0006]Φ = f (Φ1； Φ2, —, ΦΜ)
[0007]式中，f为平滑组合算子；4)采用整体保幅因子Φ作用于数据体Y，对全频段信号Y实现误差精度允许范围内的保幅，即满足:I I ΦΥ-Χ| I2 < ε，ε为信号保幅的精度。
[0008]所述步骤2)中的局部保幅因子fI3i:
【权利要求】
1.一种地震数据分形保幅方法，其包括以下步骤: 1)针对拓频前的地震数据体Xtl和拓频后的地震数据体Y，采用带通算子Bi分别相应提取它们的局部频段信息，记为BiXtl和BiY, i=l, 2，3…M，不同的i代表不同的带通频段； 2)基于地震信号局部频段成分与整体的特征存在的相似性，在每一对数据体BiXtl和BiY之间分别建立保幅关系，并相应提取局部保幅因子Oi ； 3)对所述步骤2)提取的多个局部保幅因子Oi进行平滑和组合处理，得到整体保幅因子Φ: Φ =Φ2,…，ΦΜ) 式中，f为平滑组合算子； 4)采用整体保幅因子Φ作用于数据体Y，对全频段信号Y实现误差精度允许范围内的保幅，即满足:| |ΦΥ-Χ| |2< ε，ε为信号保幅的精度。
2.如权利要求1所述的一种地震数据分形保幅方法，其特征在于:所述步骤2)中的局部保幅因子O1:
3.如权利要求1或2所述的一种地震数据分形保幅方法，其特征在于:所述步骤3)采用两种中值滤波方法中的一种对局部保幅因子01进行平滑，滤除奇异点值噪音；两种中值滤波方法具体滤波方式分别为: ①假设存在一数值序列(X1,Xy…，xn}，若
4.如权利要求1或2所述的一种地震数据分形保幅方法，其特征在于:所述步骤3)采用两种组合方式中的一种对局部保幅因子Oi进行组合，得到整体保幅因子Φ ;两组组合方式分别为: ①保幅因子Φ是各频段保幅因子的算术加权平均:
【文档编号】G01V1/28GK103630932SQ201310434890
【公开日】2014年3月12日 申请日期:2013年9月23日 优先权日:2013年9月23日
【发明者】李绪宣, 曹思远, 胡光义, 范廷恩, 董建华, 王宗俊, 樊鹏军, 蔡文涛, 张会来, 袁殿 申请人:中国海洋石油总公司, 中海油研究总院, 中国石油大学(北京)