[发明专利]一种地震数据重建方法有效
申请号: | 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:对所有分块数据进行数据重建,并将所有重建的数据拼接成完整的重建数据。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国海洋石油集团有限公司;中海油研究总院有限责任公司,未经中国海洋石油集团有限公司;中海油研究总院有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910212581.8/1.html,转载请声明来源钻瓜专利网。
- 数据显示系统、数据中继设备、数据中继方法、数据系统、接收设备和数据读取方法
- 数据记录方法、数据记录装置、数据记录媒体、数据重播方法和数据重播装置
- 数据发送方法、数据发送系统、数据发送装置以及数据结构
- 数据显示系统、数据中继设备、数据中继方法及数据系统
- 数据嵌入装置、数据嵌入方法、数据提取装置及数据提取方法
- 数据管理装置、数据编辑装置、数据阅览装置、数据管理方法、数据编辑方法以及数据阅览方法
- 数据发送和数据接收设备、数据发送和数据接收方法
- 数据发送装置、数据接收装置、数据收发系统、数据发送方法、数据接收方法和数据收发方法
- 数据发送方法、数据再现方法、数据发送装置及数据再现装置
- 数据发送方法、数据再现方法、数据发送装置及数据再现装置