[发明专利]基于压缩感知的井约束地震低频恢复方法

权利要求书

1.基于压缩感知的井约束地震低频恢复方法，包括以下步骤：

a1、井资料处理：通过井资料估计反射系数

a2、提取全局静态地震子波，构建测量矩阵方程WF；

a3、初始化参数：任意给定一个初始反射系数r0，内部循环最大迭代次数Lin，外部循环最大迭代次数Lout，残差ε，初始迭代次数k＝1；

a4、r⁰＝r0，外循环：判断||G-WFr^k||₂＞ε and k≤L_out，r⁰＝r^k-1，内循环：j＝0 to Lin-1，终止时反射系数

a5、求取拓频后频谱

a6、求恢复后频谱L、H分别为提取低频和高频频谱值的算子；

a7、低频恢复后的地震数据：

所述T_λ为阈值函数，所述(WF)^T为WF的共轭转置，所述WF为测量矩阵由子波算子和傅里叶变换算子组成，所述为R_w的共轭转置，所述R_w为取样算子提取对应井的位置信息，所述α取值大于(WF)^TWF的最大特征值，所述λ为调节反射系数的稀疏度，所述内循环以及外循环执行：k＝k+1，所述反射系数所述G为含噪数据的傅里叶变换。