[发明专利]基于傅里叶变换的地震数据重建方法及系统在审
| 申请号： | 201710892670.2 | 申请日： | 2017-09-27 |
| 公开（公告）号： | CN109557581A | 公开（公告）日： | 2019-04-02 |
| 发明（设计）人： | 霍志周;刘喜武;张远银 | 申请（专利权）人： | 中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院 |
| 主分类号： | G01V1/28 | 分类号： | G01V1/28;G01V1/36 |
| 代理公司： | 北京思创毕升专利事务所 11218 | 代理人： | 孙向民;廉莉莉 |
| 地址： | 100728 北*** | 国省代码： | 北京;11 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 地震数据 重建 傅里叶 傅里叶变换 目标函数 线性方程 最小二乘 自回归 叠加效果 叠前偏移 速度分析 应用效果 反演 稀疏 | ||
1.一种基于傅里叶变换的地震数据重建方法，包括：
基于初始地震数据，通过多步自回归方法，建立傅里叶重建线性方程；
基于所述傅里叶重建线性方程，通过最小二乘反演，建立目标函数；
基于所述目标函数，获得稀疏解，进而重建地震数据。
2.根据权利要求1所述的基于傅里叶变换的地震数据重建方法，其中，所述多步自回归方法包括：
将所述初始地震数据转换至频率域，获取近似线性的同相轴地震数据；
基于所述同相轴地震数据，通过自回归模型，获取谐波函数叠加的地震数据；
基于所述同相轴地震数据与所述谐波函数叠加的地震数据，获得地震数据概率，进而建立所述傅里叶重建线性方程。
3.根据权利要求2所述的基于傅里叶变换的地震数据重建方法，其中，所述同相轴地震数据为：
其中，S(mΔx,nΔf)表示同相轴地震数据，Δx表示空间域采样间隔，Δf表示频率域采样间隔，pj表示第j个线性同相轴的斜率，Aj表示振幅，m表示地震数据的道数，n表示每道地震数据的采样点数。
4.根据权利要求3所述的基于傅里叶变换的地震数据重建方法，其中，所述谐波函数叠加的地震数据为：
其中，L表示谐波函数的数量，P(j,nΔf)表示预测滤波因子。
5.根据权利要求4所述的基于傅里叶变换的地震数据重建方法，其中，所述地震数据概率为：
P′(j,nΔf/α)＝P(j,nΔf)j＝1,2,…,L (3)
其中，P′(j,nΔf/α)表示地震数据概率。
6.根据权利要求1所述的基于傅里叶变换的地震数据重建方法，其中，所述目标函数为：
其中，J表示目标函数，d表示数据空间的向量，m表示模型空间的向量，矩阵A表示傅里叶反变换运算，表示噪音的协方差矩阵，表示先验模型的协方差矩阵，m0表示先验模型的初始解。
7.根据权利要求1所述的基于傅里叶变换的地震数据重建方法，其中，所述稀疏解为：
其中，W为权系数组成的对角阵。
8.一种基于傅里叶变换的地震数据重建系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
基于初始地震数据，通过多步自回归方法，建立傅里叶重建线性方程；
基于所述傅里叶重建线性方程，通过最小二乘反演，建立目标函数；
基于所述目标函数，获得稀疏解，进而重建地震数据。
9.根据权利要求8所述的基于傅里叶变换的地震数据重建系统，其中，所述多步自回归方法包括：
将所述初始地震数据转换至频率域，获取近似线性的同相轴地震数据；
基于所述同相轴地震数据，通过自回归模型，获取谐波函数叠加的地震数据；
基于所述同相轴地震数据与所述谐波函数叠加的地震数据，获得地震数据概率，进而建立所述傅里叶重建线性方程。
10.根据权利要求8所述的基于傅里叶变换的地震数据重建系统，其中，所述目标函数为：
其中，J表示目标函数，d表示数据空间的向量，m表示模型空间的向量，矩阵A表示傅里叶反变换运算，表示噪音的协方差矩阵，表示先验模型的协方差矩阵，m0表示先验模型的初始解。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院，未经中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201710892670.2/1.html，转载请声明来源钻瓜专利网。
- 上一篇：一种断层的地震解释方法
- 下一篇：一种二维多分量地震资料偏移成像方法及系统
