基于Hankel张量分解的地震信号重构方法技术

本发明专利技术公开了一种基于Hankel张量分解的地震信号重构方法，其包括构造Hankel化的目标函数，将Hankel化的目标函数转换为矩阵补全目标函数，采用交替最小化方法轮换求解。本发明专利技术通过引入低秩张量分解，将Hankel构造和低秩张量有效的结合起来，解决了仅仅依靠张量低秩分解进行正则化，对数据的低秩性要求过高的问题；同时，通过低秩张量分解的方式，避免了求张量奇异值分解过程，求解速度得到了很大提升。

【技术实现步骤摘要】
基于Hankel张量分解的地震信号重构方法
本专利技术地震数据处理
，具体涉及一种基于Hankel张量分解的地震信号重构方法。
技术介绍
随着步入21世纪，全球经济进入一个新的发展期，各行各业对能源的需求日益强烈，然而新能源并没有能够替代现在石油、天然气在能源体系的核心地位。伴随着全球经济高速发展，所有国家都对石油、天然气严重依赖，其需求日益增加，国内对石油天然气这种关系民生的资源的消耗更是达到上亿吨的数量。虽然我国的国土面积大，但是相对于中东而言，我国并不是处于石油带上，只有西藏的部分地区在全球的石泊带的尾巴上，这就造成我国可供开采的石油天然气、页岩气不足:随着不断的开采，我们引以为傲的大庆油田的石泊储量也只能使用几十年。在对能源的强烈需求催使下，我们必须加大对石油、天然气的勘探力度，对潜在的薄层中的油气资源需要倾注更多的力量。尤其随着我国石油勘探在国家的中西部地区进行，而这些地方往往拥有复杂的地表环境，石泊在地下岩层中埋藏较深，因为地壳活动的影响其分布并不可能是一大块区域，而是一片一片，造成勘测难度极大。经济成本也会相应的增加，同时勘探中存在的风险也是难以想象的:相对较深的本文档来自技高网...

【技术保护点】
1.一种基于Hankel张量分解的地震信号重构方法，其特征在于，包括以下步骤：A、对缺失叠后地震数据进行Hankel重构，构造Hankel化的目标函数；B、将Hankel化的目标函数分解为k个沿时间维度的独立子问题，再将每一个子问题降解为矩阵补全目标函数；C、采用交替最小化方法对步骤B中目标函数进行轮换求解，得到最优拟合张量。

【技术特征摘要】
1.一种基于Hankel张量分解的地震信号重构方法，其特征在于，包括以下步骤：A、对缺失叠后地震数据进行Hankel重构，构造Hankel化的目标函数；B、将Hankel化的目标函数分解为k个沿时间维度的独立子问题，再将每一个子问题降解为矩阵补全目标函数；C、采用交替最小化方法对步骤B中目标函数进行轮换求解，得到最优拟合张量。2.如权利要求1所述的基于Hankel张量分解的地震信号重构方法，其特征在于，所述步骤A中，对缺失叠后地震数据进行Hankel重构，具体为：通过低秩张量因子分解将缺失叠后地震数据分解为两个低秩张量的张量乘积，再将时域的张量乘积转换为频域的前切片矩阵乘积，并进行Hankel重构。3.如权利要求2所述的基于Hankel张量分解的地震信号重构方法，其特征在于，所述步骤A中，Hankel化的目标函数表示为其中，u和v为缺失叠后地震数据进行低秩分解后的低秩张量，和分别为u和v的傅里叶域形式，y为观测数据，为y的傅里叶域形式，和表示实数域，PΩ为在观测集Ω上的投影操作符，⊙为以元素为单位的乘积，为前切片乘积，||·||F为Frobenius范数，H(·)为张量Hankel化操作符。4.如权利要求3所述的基于Hankel张量分解的地震信号重构方法，其特征在于，所述步骤B中，矩阵补全目标函数表示为其中，U和V为低秩因子矩阵，Y为观测数据矩阵。5.如权利要求4所述的基...

【专利技术属性】
技术研发人员：钱峰，邓弢，陈全，胡光岷，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51