The invention discloses a method and system for separating effective wave from multiple wave in seismic data, which includes constructing a mixed objective function based on L1 norm and L2 norm, obtaining initial effective wave data and initial multiple data in seismic data based on the mixed objective function, and determining weighting coefficients according to the initial effective wave data and the initial multiple data. The optimal matching operator is obtained from the weighting coefficients and the mixed objective function, and the effective wave and multiple wave in the seismic data are separated according to the optimal matching operator. This method can sufficiently adapt to the diversity of multiple waves caused by complex and diverse geological conditions in actual seismic data, and can effectively improve the signal-to-noise ratio of pre-stack seismic data, and provide high-quality pre-stack data for subsequent seismic data processing.
【技术实现步骤摘要】
用于分离地震数据中有效波与多次波的方法及系统
本专利技术属于地震数据处理领域,尤其涉及一种用于分离地震数据中有效波与多次波的方法及系统。
技术介绍
压制与利用多次波是勘探地震学中研究的热点和难点问题之一。在油气勘探中,地震多次波普遍存在。由于多次波的存在,给速度分析、偏移成像带来极大的困难,并将导致假的成像构造,这将严重影响地震解释工作。针对地震有效波与多次波的分离问题,目前地球物理学者已研究出多种自适应分离方法,主要归纳为以下两类:一是基于L2范数的地震多次波自适应分离方法。该类方法基于剩余能量最小准则,在最小二乘意义下求取自适应匹配算子,然后从原始数据中减去匹配后的多次波。一是基于L1范数的地震多次波自适应分离方法。该方法假定地震有效波具有最小L1范数,即地震有效波具有很好的稀疏性。当地震有效波不满足剩余能量最小准则时,该方法可在实现有效波与多次波分离的同时,也能相对更好的保持有效波能量。对于实际地震数据,由于产生多次波的地质条件复杂多样,不仅有低阶次多次波,还有高阶次多次波,甚至有层间多次波。关于地震有效波具有最小L2范数或者L1范数假设条件常常并不成立,而且有效波与多次波的能量对比随时间和空间变化很大,现有技术中的分离方法经常会导致压制多次波后的地震数据中剩余较强的多次波能量。当地震有效波与多次波同相轴重叠时,压制多次波的自适应减过程甚至会伤害地震有效波,影响地震数据处理的保真度。针对上述复杂情况下的多次波分离问题,目前还没有一种方法能够较好解决。
技术实现思路
本专利技术所要解决的技术问题之一是需要提供一种用于实现复杂情况下的多次波分离有效波与多次波 ...
【技术保护点】
1.一种用于分离地震数据中有效波与多次波的方法,包括:基于L1范数与L2范数构建混合目标函数;基于所述混合目标函数获取地震数据中的初始有效波数据与初始多次波数据;根据所述初始有效波数据与所述初始多次波数据确定加权系数;基于所述加权系数以及所述混合目标函数获取最优匹配算子,并根据所述最优匹配算子分离得到所述地震数据中的有效波与多次波。
【技术特征摘要】
1.一种用于分离地震数据中有效波与多次波的方法,包括:基于L1范数与L2范数构建混合目标函数;基于所述混合目标函数获取地震数据中的初始有效波数据与初始多次波数据;根据所述初始有效波数据与所述初始多次波数据确定加权系数;基于所述加权系数以及所述混合目标函数获取最优匹配算子,并根据所述最优匹配算子分离得到所述地震数据中的有效波与多次波。2.根据权利要求1所述的方法,其特征在于,根据如下表达式构建所述混合目标函数Q(f):其中,d为原始的地震数据,M为预测的多次波数据,f为匹配滤波算子,λ为加权系数,||·||1表示L1范数,表示L2范数。3.根据权利要求2所述的方法,其特征在于,所述基于所述混合目标函数获取地震数据中的初始有效波数据与初始多次波数据,包括:在[0,1]的取值范围内任取一个值作为所述加权系数的初始值;采用迭代加权最小二乘算法将所述混合目标函数Q(f)重新定义为:其中,W为加权矩阵,且其初始值取为单位矩阵I;基于所述加权系数以及所述加权矩阵的初始值,采用最小二乘算法求取所述混合目标函数Q(f)′的最小值,进而得到匹配算子的初始值;根据所述匹配算子的初始值获取所述地震数据中的初始有效波数据与初始多次波数据。4.根据权利要求3所述的方法,其特征在于,所述加权系数的初始值取为1。5.根据权利要求2至4中任一项所述的方法,其特征在于,所述根据所述初始有效波数据与所述初始多次波数据确定加权系数,包括:根据所述初始有效波数据与所述初始多次波数据获取所述有效波与所述多次波的能量比,并基于所述能量比确定所述加权系数。6.根据权利要求5所述的方法,其特征...
【专利技术属性】
技术研发人员:薛永安,黄新武,王勇,陈习峰,张浩东,汤国松,庞全康,刘立民,管文华,潘成磊,
申请(专利权)人:中国石油化工股份有限公司,中国石油化工股份有限公司江苏油田分公司,
类型:发明
国别省市:北京,11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。