[发明专利]基于快速傅立叶变换的虚源重构方法及系统有效
| 申请号： | 201811139507.X | 申请日： | 2018-09-28 |
| 公开（公告）号： | CN110967734B | 公开（公告）日： | 2022-03-08 |
| 发明（设计）人： | 雷朝阳;陈国金;吴永栓 | 申请（专利权）人： | 中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院 |
| 主分类号： | G01V1/28 | 分类号： | G01V1/28 |
| 代理公司： | 北京思创毕升专利事务所 11218 | 代理人： | 孙向民;廉莉莉 |
| 地址： | 100027 北*** | 国省代码： | 北京;11 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 基于 快速 傅立叶 变换 虚源重构 方法 系统 | ||
1.一种基于快速傅立叶变换的虚源重构方法，其特征在于，包括：
步骤1：输入地震数据，确定虚源A与虚源检波器B；
步骤2：利用快速傅立叶变换，分别计算在每一个实际震源激发下对应的A与B的波场记录的互相关值；
步骤3：将所有的互相关值进行叠加，得到A点激发B点接收的虚源波场；
步骤4：重复进行步骤1-3，完成所有虚源波场的计算；
其中，所述步骤2包括：
根据每一个实际震源激发下对应的A与B的波场记录，获得时间序列x与h；
确定x与h的采样点数均为M，构造时间序列其中，采样点数均为N＝2*M；
对所述时间序列进行快速傅立叶变换得到频谱
将的振幅谱相乘，相位谱相减，得到互相关的频谱
对所述频谱进行反快速傅立叶变换得到循环相关结果
根据所述循环相关结果计算A与B的波场记录的互相关值。
2.根据权利要求1所述的基于快速傅立叶变换的虚源重构方法，其中，通过公式(1)计算所述时间序列
其中，xi为时间序列x中的第i个参量。
3.根据权利要求1所述的基于快速傅立叶变换的虚源重构方法，其中，通过公式(2)计算所述时间序列
其中，hi为时间序列h中的第i个参量。
4.根据权利要求1所述的基于快速傅立叶变换的虚源重构方法，其中，通过公式(3)计算所述A与B的波场记录的互相关值：
其中，rxh(i)为互相关值，为循环相关结果，N为的采样点数。
5.一种基于快速傅立叶变换的虚源重构系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
步骤1：输入地震数据，确定虚源A与虚源检波器B；
步骤2：利用快速傅立叶变换，分别计算在每一个实际震源激发下对应的A与B的波场记录的互相关值；
步骤3：将所有的互相关值进行叠加，得到A点激发B点接收的虚源波场；
步骤4：重复进行步骤1-3，完成所有虚源波场的计算；
其中，所述步骤2包括：
根据每一个实际震源激发下对应的A与B的波场记录，获得时间序列x与h；
确定x与h的采样点数均为M，构造时间序列其中，采样点数均为N＝2*M；
对所述时间序列进行快速傅立叶变换得到频谱
将的振幅谱相乘，相位谱相减，得到互相关的频谱
对所述频谱进行反快速傅立叶变换得到循环相关结果
根据所述循环相关结果计算A与B的波场记录的互相关值。
6.根据权利要求5所述的基于快速傅立叶变换的虚源重构系统，其中，通过公式(1)计算所述时间序列
其中，xi为时间序列x中的第i个参量。
7.根据权利要求5所述的基于快速傅立叶变换的虚源重构系统，其中，通过公式(2)计算所述时间序列
其中，hi为时间序列h中的第i个参量。
8.根据权利要求5所述的基于快速傅立叶变换的虚源重构系统，其中，通过公式(3)计算所述A与B的波场记录的互相关值：
其中，rxh(i)为互相关值，为循环相关结果，N为的采样点数。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院，未经中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201811139507.X/1.html，转载请声明来源钻瓜专利网。
- 上一篇：加热装置
- 下一篇：一种基于质谱的金团簇离子反应活性分析方法
