申请/专利权人：中国地质大学(武汉)

申请日：2019-09-17

公开（公告）日：2020-11-24

公开（公告）号：CN110716231B

主分类号：G01V1/36(20060101)

分类号：G01V1/36(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.11.24#授权;2020.02.21#实质审查的生效;2020.01.21#公开

摘要：本发明公开了一种基于共聚焦域稀疏反演的海上多震源波场分离方法及系统，该方法首先，将基于伪分离后的混合炮集，在炮集记录中选取若干条同相轴，提取对应的叠加速度和初至时间；然后，基于得到的时间‑速度对，计算进行聚焦变换所需的聚焦变换算子，并对时空域地震数据体进行聚焦正变换，将时空域地震记录变换到共聚焦域；然后，在共聚焦域建立稀疏约束反演目标函数，并采用谱投影梯度算法求解目标函数，得到共聚焦域的稀疏解；最后，在共聚焦域数据体的两端进行聚焦反变换，得到去除串扰噪声的时空域地震记录，从而达到混合源分离的目的。分离后的波场不仅有效地去除了串扰噪声，而且较好的保留了原始的振幅、相位等信息，有利于后续的处理。

主权项：1.一种基于共聚焦域稀疏反演的海上多震源波场分离方法，其特征在于，包括以下步骤：S1、获取混合炮集P以及所述混合炮集中包括的每个单震源的激发延迟时间；根据所述每个单震源的激发延迟时间，对所述混合炮集进行伪分离操作；S2、在经伪分离处理后的混合炮集中选取若干条同相轴，提取每条同相轴对应的叠加速度和初至时间，得到若干个时间-速度对；S3、先将混合炮集进行傅里叶正变换，再采用聚焦变换算子对混合炮集进行聚焦正变换操作后，得到若干个共聚焦子域；其中，聚焦正变换过程中需用到的聚焦算子均基于步骤S2处理得到的时间-速度对计算得到；所述聚焦算子定义为： 上述表达式中，为经过傅里叶反变换后得到的聚焦算子；和为相移算子，Δz＝|zn-zm|，其中，zi为炮点或检波点与焦点在垂直方向上的距离，c、t为基于同相轴得到的对应时间-速度对；其中，k是定义为k＝ωc的波数，ω和c分别代表角频率和时间-速度对中对应的纵波速度；kz是纵向波数，kx是横向波数，表示对“·”进行傅里叶反变换；S4、基于步骤S3得到的所述若干个共聚焦子域，建立多震源波场分离的稀疏约束反演目标函数，其为： 其中，Xk表示了地下深度为k及其附近的反射波场在共聚焦子域压缩的信息；σ表示混合炮集中的噪声大小，ω表示角频率；对于元素为Aij的矩阵A，||A||S和||A||F分别表示和范数以及弗罗贝尼乌斯范数，其中，||A||S＝∑i,j|Aij|，i,j分别表示元素A在矩阵中对应的行数、列数；K为聚焦算子的数量；是将地表震源的波场外推到地下深度为zk的正传算子，Wk-是将地下深度为zk的波场外推到地表检波点的反传算子；S5、采用谱投影梯度算法对所述稀疏约束反演目标函数进行求解，得到去除串扰噪声的共聚焦域数据；S6、针对所述去除串扰噪声的共聚焦域数据，采用聚焦变换算子对其进行聚焦反变换，再进行傅里叶反变换后，将其还原为时间-空间域的数据体，实现混合源波场分离。

全文数据：

权利要求：

百度查询： 中国地质大学(武汉) 基于共聚焦域稀疏反演的海上多震源波场分离方法及系统

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD