买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：广东海洋大学

摘要：本发明属于海洋地震勘探技术领域，具体涉及一种针对低信噪比被动源虚炮集的鬼波压制方法，解决了虚炮集数据中的低信噪比与虚炮集与常规拖缆观测系统不同的问题。本发明鬼波压制方法分别从提升去鬼波算子的抗噪性和增加先验约束两方面入手达到提升方法稳定性的目的；针对噪声干扰问题，通过对数据体进行变观操作，构建数据矩阵，满足引入稀疏约束的条件，并通过构建稀疏约束鬼波方压制法的目标函数，以提升反演的准确性；另一方面，从反演类方法本身的稳定性入手，本发明将反演形式由原来的单炮循环变成了多炮同时反演，并将地层的连续性作为了反演过程中新的约束信息，使得反演结果的精度进一步提升，极大的扩充了鬼波压制算法的适用性。

主权项：1.一种针对低信噪比被动源虚炮集的鬼波压制方法，其特征在于，包括以下步骤：A、利用被动源数据Tobs干涉重构出系列单道记录Dvxi|xj,t，重新排序为以xi作为虚震源位置的共炮点虚炮记录Dsxi,t，并按炮点位置xi增大方向进行排列，获得炮点和检波点位置相同的3D数据体Dxi,xj,t；B、对步骤A所得数据体Dxi,xj,t按照共偏移距方向进行变观操作，以改变观测方位，获得数据D'xoff,xk,t；变观过程，通过如下公式表示为： 其中，D′offxoff,t＝Dxi,xj,t，4xoff＝xi-xj，5其中，D'xoff,xk,t表示变观后的3D数据体，xoff表示用偏移距进行编号的炮号，其取值范围[-xn,xn]，xk为检波点坐标，k∈1,n；C、以D't简化表示D'xoff,xk,t，对变观后的数据体D't沿时间方向进行傅里叶变换，得到其频率域形式数据体D'ω；D、由频率域数据体D'ω提取切片获取频率域切片，构建稀疏约束鬼波压制方法的目标函数FS；所述目标函数为公式7： 其中，S表示压制鬼波后的数据，λ为正则化算子，||||F表示求解矩阵的Frobenius范数，||S||1表示从第一列开始排列矩阵S的所有元素所获得的L1范数；E、通过引入低秩约束来引入地层连续性作为结果的先验约束，得到包含先验信息约束的目标函数，进而获得稀疏核范数目标函数FS；包含先验信息约束的目标函数为公式8： 其中，μ是低秩约束的正则化参数，rank表示求矩阵的秩，根据矩阵完备理论，使用以下目标函数作为最优凸逼近问题来代替上述公式，得到公式9 其中，是矩阵S的第i个奇异值；F、使用分裂Bergman迭代策略求解步骤E所得目标函数FS，引入辅助变量矩阵Θ，将目标函数FS调整为FS,Θ；G、设置参数初始值；H、迭代求解，并利用软阈值收缩算法求解目标函数，以获得辅助矩阵Θ；I、设置Θ＝Θl，并利用AMP算法求解目标函数，以获得第l次迭代的去鬼波结果Sl；J、判断鬼波压制结果Sl是否满足收敛条件，满足判断条件，输出经过l次迭代过程后的压制鬼波结果Sl；否则，重新设置参数，并返回步骤I。

全文数据：

权利要求：

百度查询： 广东海洋大学 一种针对低信噪比被动源虚炮集的鬼波压制方法