地震数据保幅重建方法和系统技术方案

本发明专利技术提供地震数据保幅重建方法和系统，所述方法包括：获取原始缺失地震数据，并对原始缺失地震数据进行预处理得到动校正数据；对动校正数据进行频域变换得到第一频域数据；利用快速高阶抛物Radon变换算法对第一频域数据进行计算得到第二频域数据；计算权系数，并对第二频域数据进行加权计算得到频域重建数据；判断频域重建数据是否满足预设条件；如果不满足预设条件，则将频域重建数据进行迭代计算得到时域重建数据。本发明专利技术可以改进远偏移距数据的重建效果，并提高方法的计算效率。

【技术实现步骤摘要】
地震数据保幅重建方法和系统
本专利技术涉及商用厨房灭火
，尤其是涉及地震数据保幅重建方法和系统。
技术介绍
在地震数据重建方面，一般有如下方法：基于预测滤波重建方法，基于波动方程重建方法和基于函数变换重建方法等。基于函数变换这类方法的特点是利用数学变换理论和信号分析原理，如Fourier变换、Wavelet变换、Curvelet变换、Radon变换和压缩感知等对数据进行重建。Radon变换又包括线性Radon、双曲Radon和抛物Radon这三种。在地震勘探中，利用水平层状介质或小倾角地质结构的反射同相轴走时特征在CMP道集上是双曲线的性质，双曲Radon变换在地震数据多次波压制、插值重建理论研究方面得到广泛关注，但是其计算效率的低下限制了双曲Radon变换在工业上的应用。在基于抛物Radon变换地震数据重建方面，传统的抛物Radon变换是同相轴振幅沿着不同的曲率进行叠加，没有考虑振幅随偏移距的变化特征，因此在处理实际地震数据时不具备保幅性；同时，以上方法中的迭代加权系数在不同位置的缺失道都相同，而实际上不同位置缺失的重建效果是不一样的，用相同的权系数会导致较大误差出现。专利技本文档来自技高网...

【技术保护点】
一种地震数据保幅重建方法，其特征在于，包括：获取原始缺失地震数据，并对所述原始缺失地震数据进行预处理得到动校正数据；对所述动校正数据进行频域变换得到第一频域数据；利用快速高阶抛物Radon变换算法对所述第一频域数据进行计算得到第二频域数据；计算权系数，并对所述第二频域数据和所述权系数进行加权计算得到频域重建数据；判断所述频域重建数据是否满足预设条件；如果不满足所述预设条件，则将所述频域重建数据进行迭代计算得到时域重建数据。

【技术特征摘要】
1.一种地震数据保幅重建方法，其特征在于，包括：获取原始缺失地震数据，并对所述原始缺失地震数据进行预处理得到动校正数据；对所述动校正数据进行频域变换得到第一频域数据；利用快速高阶抛物Radon变换算法对所述第一频域数据进行计算得到第二频域数据；计算权系数，并对所述第二频域数据和所述权系数进行加权计算得到频域重建数据；判断所述频域重建数据是否满足预设条件；如果不满足所述预设条件，则将所述频域重建数据进行迭代计算得到时域重建数据。2.根据权利要求1所述的地震数据保幅重建方法，其特征在于，所述频域变换包括傅里叶变换，所述利用快速高阶抛物Radon变换算法对所述第一频域数据进行计算得到第二频域数据包括：根据快速高阶抛物Radon反变换公式对所述第一频域数据进行反演求解得到变换域系数，其中，所述反演求解包括利用最小二乘框架下的奇异值分解SVD方法进行矩阵求逆计算；将所述变换域系数代入所述快速高阶抛物Radon反变换公式得到所述第二频域数据。3.根据权利要求1所述的地震数据保幅重建方法，其特征在于，所述计算权系数包括：根据下式计算权系数：其中，W(Tn)为所述权系数，Tn为第n道数据，xn，yn为缺失道的坐标；xmax，ymax为最大偏移距；γ1和γ2取值范围为[0,1]；an表示所述缺失道所处位置的缺失程度。4.根据权利要求1所述的地震数据保幅重建方法，其特征在于，所述预设条件为所述频域重建数据满足所述迭代计算的次数要求，所述将所述频域重建数据进行迭代计算包括：将所述频域重建数据填补入原始缺失道；重复进行所述快速高阶抛物Radon变换算法和所述加权计算，直至满足所述次数要求，其中，所述次数要求为三次；对满足所述次数的所述频域重建数据进行傅里叶反变换得到所述时域重建数据。5.根据权利要求2所述的地震数据保幅重建方法，其特征在于，根据快速高阶抛物Radon反变换公式对所述第一频域数据进行反演求解得到变换域系数包括：根据下式计算所述变换域系数：其中，M(f,λx,λy)为所述变换域系数，分别为X，Y方向的所述快速高阶抛物Radon变换算子，H为矩阵的共轭转置，...

【专利技术属性】
技术研发人员：毛伟建，唐欢欢，王亮亮，詹毅，张文波，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，中国科学院测量与地球物理研究所，
类型：发明
国别省市：河北,13