[发明专利]一种地震数据重建方法有效

专利信息
申请号: 201910212581.8 申请日: 2019-03-20
公开(公告)号: CN110161565B 公开(公告)日: 2021-05-11
发明(设计)人: 黄小刚;曹静杰;张金淼;杜向东;李维新;朱振宇;翁斌;刘志斌;张益明;丁继才;薛东川;王小六 申请(专利权)人: 中国海洋石油集团有限公司;中海油研究总院有限责任公司
主分类号: G01V1/36 分类号: G01V1/36;G01V1/28
代理公司: 北京纪凯知识产权代理有限公司 11245 代理人: 谢斌
地址: 100010 北*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 地震 数据 重建 方法
【权利要求书】:

1.一种地震数据重建方法,其特征在于,包括以下步骤:

步骤A:获取待重建的地震数据及其采集矩阵;

步骤B:对地震数据和采集矩阵进行分块处理;

步骤C:建立每个小块的地震数据及其对应的采集矩阵的双曲正切函数正则化的稀疏反演模型,其中,所述双曲正切函数正则化的稀疏反演模型为:

式中,min表示最小化符号,Jσ(x)为目标函数,F-1为傅里叶变换的逆,x为地震数据在傅里叶变换域的系数,:=表示将minJσ(x)定义为s.t.表示前面目标函数最小化的约束条件为SF-1x=d,其中σ为超参数,d为待重建的地震数据,S为采集矩阵,表示地震数据在傅里叶变换域系数的第i个分量的一种双曲正切函数;N为地震数据在傅里叶变换域的系数x的分量个数,xi表示地震数据在傅里叶变换域的系数x的第i个分量;

步骤D:通过每个小块数据求取该稀疏反演模型的最优解,得到重建后的地震数据;

步骤D21:设定循环次数IterN、初始超参数σ0和终止迭代超参数σIterN,并将观测得到的地震数据作为初始解d0

步骤D22:计算当前第k次迭代的迭代解在变换域的系数xk,其表达式为:

xk=Fdk (2)

式中,k表示迭代序号;F表示傅里叶变换;当k=0时,表示初始迭代,dk表示第k次迭代时重建所得的地震数据;

计算当前第k次迭代的修正牛顿方向gk,其表达式为:

式中,σk为第k次迭代的超参数;

求出更新解:

式中,λ为步长;

步骤D23:更新解在时间空间域的投影:

进行k的自增1运算,如果k<IterN,重复所述步骤D22到所述步骤D23的迭代过程;如果k≥IterN,则迭代结束;

步骤D24:输出小块数据的最后一次迭代时重建所得的地震数据dIterN

步骤E:对所有分块数据进行数据重建,并将所有重建的数据拼接成完整的重建数据。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国海洋石油集团有限公司;中海油研究总院有限责任公司,未经中国海洋石油集团有限公司;中海油研究总院有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201910212581.8/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top