基于时空变系数加权的优化叠加方法及系统技术方案

本发明专利技术提出了一种基于时空变系数加权的优化叠加方法及系统，该方法包括：由原始数据叠加求取参考序列；计算参考序列和叠前道集间的局部时空变相关系数；确定时空变的加权系数；将时空变的加权系数应用于对应的叠前道集；将加权滤波后的叠前道集进行叠加，得到最终叠加剖面。本发明专利技术介绍了一种利用叠前道集与参考序列集之间的时空变相关系数作为加权系数的优化叠加方法，加权系数滤波处理后的叠前道集中的噪声干扰得到有效压制，有效信号得到增强，提高了地震成像道集和叠加剖面的信噪比和整体质量。

Optimal Superposition Method and System Based on Weighting of Spatiotemporal Variable Coefficient

The present invention proposes an optimal stacking method and system based on weighted spatiotemporal variable coefficients, which includes: obtaining reference sequences by superimposing original data; calculating local spatiotemporal variable correlation coefficients between reference sequences and prestack gathers; determining weighted coefficients of spatiotemporal variation; applying weighted coefficients of spatiotemporal variation to corresponding prestack gathers; and stacking pre-stack gathers after weighted filtering. Finally, the superimposed section is obtained. The present invention introduces an optimal stacking method using the spatiotemporal correlation coefficient between prestack gathers and reference sequence sets as weighting coefficients. The noise interference in prestack gathers after weighted coefficient filtering is effectively suppressed, the effective signal is enhanced, and the signal-to-noise ratio and the overall quality of seismic imaging gathers and stacked sections are improved.

【技术实现步骤摘要】
基于时空变系数加权的优化叠加方法及系统
本专利技术属于地震勘探资料处理
，涉及高精度地震成像及地震解释处理技术，具体地涉及基于时空变系数加权的优化叠加方法及系统。
技术介绍
面对复杂地区的勘探开发，获得正确的构造认识和油气有利圈闭对地震解释工作的精度要求越来越高。由于得到共反射点道集的过程是在叠前实现的，因此其信噪比一般比较低，低信噪比成像道集的直接叠加会降低成像叠加剖面的质量，困扰着传统解释工作的进行。不能被现有方法有效利用的噪声干扰是造成信噪比低的原因之一，并且以随机噪声居多，可以采用叠前去除随机噪声方法来消除或压制，如f-k域倾向滤波、f-x域预测滤波、中值滤波等，但是这些滤波方法没有方向特性，易造成对倾斜同相轴振幅的压制，同时也会模糊小的断层和裂缝等不连续结构。结构滤波是一种具有地层倾向的滤波方法，能够用于叠前资料的噪声压制，并且具有保护边缘结构的良好特性，但如何有效地对地层倾向进行估计是该方法的关键和难点。
技术实现思路
本专利技术介绍了一种利用叠前道集与参考序列集之间的时空变相关系数作为加权系数的优化叠加方法，加权系数滤波处理后的叠前道集中的噪声干扰得到有效压制，有效信号得到增强，提高了地震成像道集和叠加剖面的信噪比和整体质量，对于推动更高精度的地震成像具有重要的现实意义。针对地震成像道集中存在的噪声干扰问题，本专利技术介绍了一种利用叠前道集与参考序列集之间的时空变相关系数作为加权系数的优化叠加方法，加权系数滤波处理后的叠前道集中的噪声干扰得到有效压制，有效信号得到增强，提高了地震成像道集和叠加剖面的信噪比和整体质量。根据本专利技术的一个方面，提供一种基于时空变系数加权的优化叠加方法，该方法包括：由原始数据叠加求取参考序列；计算参考序列和叠前道集间的局部时空变相关系数；确定时空变的加权系数；将时空变的加权系数应用于对应的叠前道集；将加权滤波后的叠前道集进行叠加，得到最终叠加剖面。进一步地，所述参考序列的求取公式为：式中，N为叠加道数，j为采样点数，αi,j(t)为叠前道集数据。进一步地，参考序列和叠前道集间的局部时空变相关系数表示为：c1＝[λ2I+S(ATA-λ2I)]-1SATb,c2＝[λ2I+S(BTB-λ2I)]-1SBTa.式中，A和B分别表示由a和b元素组成的对角矩阵算子，S为整形算子，a和b为两个离散信号的矢量表示，λ为正则化控制算子，c1和c2的积的各个分量即为两个信号的局部时空变相关系数。进一步地，时空变的加权系数的优化叠加公式表示为：式中，ωi,j(t)为加权系数，ε为局部相关系数的阈值，为第j个共中心点道集权重的和，Hj(t)为wi,j·ai,j(t)≠0的采样点的个数，ηi,j(t)为第j个共中心点道集第i道的叠前道集与参考序列的局部相关系数。进一步地，将时空变的加权系数应用于对应的叠前道集的公式为：式中，ωi,j(t)为加权系数，αi,j(t)为叠前道集数据。根据本专利技术的另一方面，提供一种基于时空变系数加权的优化叠加系统，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：由原始数据叠加求取参考序列；计算参考序列和叠前道集间的局部时空变相关系数；确定时空变的加权系数；将时空变的加权系数应用于对应的叠前道集；将加权滤波后的叠前道集进行叠加，得到最终叠加剖面。通过对理论模型和实际资料的处理，验证了本方法对叠前道集中非相关性干扰具有很好的压制效果，有效提高了叠前道集和叠加剖面的质量。理论模型及实际地震数据处理结果表明：对叠前道集进行基于本方法的时空变加权系数优化叠加处理，有效压制了叠前道集中的随机噪声和非相关性噪声的干扰，提高了叠前道集和叠加剖面的信噪比与整体质量。附图说明通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术实施例的处理流程图。图2示出了合成理论模型道集及对应的相关系数。图3示出了常规叠加与本专利技术得到的叠加道的对比。图4示出了实际资料原始成像道集。图5示出了本专利技术得到的加权滤波后的成像道集。图6示出了原始叠加剖面。图7示出了本专利技术处理后得到的叠加剖面。具体实施方式下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。本专利技术属于地震勘探资料处理
，涉及高精度地震成像及地震解释处理技术。成像道集中由于存在随机噪声和非相关性噪声的干扰使得道集信噪比较低，影响了叠加剖面的成像质量。本专利技术介绍了一种利用叠前道集与参考序列集之间的时空变相关系数作为加权系数的优化叠加方法，加权系数滤波处理后的叠前道集中的噪声干扰得到有效压制，有效信号得到增强，提高了地震成像道集和叠加剖面的信噪比和整体质量。常规地震数据处理中，叠加是将多个地震数据道集求和得到单个地震道，使相干信号得到加强，压制随机噪音。本专利技术提出了一种利用叠前道集和参考序列集计算出时空变相关系数，将此系数作为参与叠加的叠加道的权重，对道集进行叠加，对提高道集信噪比有良好的应用效果。两个长度为N的离散信号ai和bi的全局相关系数可定义为：上述公式计算了两列离散信号的一个相关系数，不具有时空变的特征，只能衡量全局的相似性。为了分析信号相似性的局部变化，需要将相似系数转化为时空变量的函数，因此，需要定义滑动时窗相关系数衡量两个信号的局部相似性，其相关系数定义为：其中w为滑动时窗的长度。相关系数的线性代数表达可以写为两个最小平方逆γ1与γ2的积：γ2＝γ1γ2其中a和b为ai和bi的矢量表示。令A和B分别表示由a和b的元素所组成的对角矩阵算子，对上述方程的反演过程加入整形正则化项(Fomel，2007b)，标量γ1和γ2转换为矢量c1和c2，定义为：c1＝[λ2I+S(ATA-λ2I)]-1SATb,c2＝[λ2I+S(BTB-λ2I)]-1SBTa.λ控制算子A和B的相对大小，S是一种整形算子，c1和c2的积的各个分量即为两个信号的局部相似系数。数学上的局部相关系数法优化叠加公式如下：ε为局部相关系数的阈值，为第j个共中心点道集权重的和，Hj(t)为wi,j·ai,j(t)≠0的采样点的个数，ηi,j(t)为第j个共中心点道集第i道的叠前道集与参考序列的局部相关系数，其为旅行时和偏移距的函数，参考序列由一个共中心点道集平均叠加得到。局部相关系数阈值ε的应用，使得所有低于阈值的局部相关系数都不参与叠加，只有与参考序列的相似性相对较大的一部分道集参与叠加，这个过程可以看做是提高同相轴相干性的一种非线性滤波。最后，再对所有共中心点道集或共成像点道集进行加权叠加得到最终的叠加剖面。根据本专利技术的一个专利技术，提供一种基于时空变系数加权的优化叠加方法，该方法包括：由原始数据叠加求取参考序列；计算参考序列和叠前道集间的局部时空变相关系数；确定时空变的加权系数；将时空变的加权系数应用于对应的叠前道集；将加权滤波后的叠前道集进行叠加，得到最本文档来自技高网...

【技术保护点】
1.一种基于时空变系数加权的优化叠加方法，其特征在于，该方法包括：由原始数据叠加求取参考序列；计算参考序列和叠前道集间的局部时空变相关系数；确定时空变的加权系数；将时空变的加权系数应用于对应的叠前道集；将加权滤波后的叠前道集进行叠加，得到最终叠加剖面。

【技术特征摘要】
1.一种基于时空变系数加权的优化叠加方法，其特征在于，该方法包括：由原始数据叠加求取参考序列；计算参考序列和叠前道集间的局部时空变相关系数；确定时空变的加权系数；将时空变的加权系数应用于对应的叠前道集；将加权滤波后的叠前道集进行叠加，得到最终叠加剖面。2.根据权利要求1所述的基于时空变系数加权的优化叠加方法，其特征在于，所述参考序列的求取公式为：式中，N为叠加道数，j为采样点数，αi,j(t)为叠前道集数据。3.根据权利要求1所述的基于时空变系数加权的优化叠加方法，其特征在于，参考序列和叠前道集间的局部时空变相关系数表示为：c1＝[λ2I+S(ATA-λ2I)]-1SATb,c2＝[λ2I+S(BTB-λ2I)]-1SBTa.式中，A和B分别表示由a和b元素组成的对角矩阵算子，S为整形算子，a和b为两个离散信号的矢量表示，λ为正则化控制算子，c1和c2的积的各个分量即为两个信号的局部时空变相关系数。4.根据权利要求1所述的基于时空变系数加权的优化叠加方法，其特征在于，时空变的加权系数的优化叠加公式表示为：式中，ωi,j(t)为加权系数，ε为局部相关系数的阈值，为第j个共中心点道集权重的和，Hj(t)为wi,j·ai,j(t)≠0的采样点的个数，ηi,j(t)为第j个共中心点道集第i道的叠前道集与参考序列的局部相关系数。5.根据权利要求1所述的基于时空变系数加权的优化叠加方法，其特征在于，将时空变的加权系数应用于对应的叠前道集的公式为：式中，ωi,j(t)为加权系数，αi,j(t)为叠前道集数据。6.一种基于时空变系数加权的优化叠加系统，其特征在...

【专利技术属性】
技术研发人员：许璐，李博，王小青，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11