基于分数阶傅里叶域的子波提取方法技术

提供了一种基于分数阶傅里叶域的子波提取方法，包括：建立地层信息初始模型；利用输入子波和地层信息初始模型的反射系数合成地震记录；计算输入子波和合成的地震记录数据的互相关矩阵以及地震记录数据的自相关矩阵，并对所述互相关矩阵和自相关矩阵进行分数阶傅里叶变换；对所述互相关矩阵和自相关矩阵应用分数阶傅里叶域的最优滤波算法来计算均方误差最小的最优阶和最优滤波算子；基于所述最优滤波算子和最优阶提取子波。

【技术实现步骤摘要】

本专利技术涉及地震勘探
更具体地讲，涉及一种在地震勘探中使用的基于分数阶傅里叶(FRFT)域的子波提取方法。
技术介绍
地震勘探技术逐渐向高精度、高分辨率方向发展，精确地提取子波是反褶积、波阻抗反演以及正演模拟的基础。近二三十年，人们提出了很多子波提取的算法，并取得了一定的效果，但是现有的地震子波提取方法均有各自的不足，且提取的子波精度不高，使得子波提取技术成为影响地震反演技术进一步提高的关键因素之一。常规的算法都假设地层是个时不变系统，然而，由于地下地层信息随着深度的变化而改变，因此它是个时变系统，所以用线性时不变系统来描述它必然不够准确，从而造成子波提取精度较差，波阻抗反演结果 与钻井地质不吻合，储层预测困难。因此，要寻找一种方法，用时变系统来描述地下地层信息，弥补常规傅里叶域不能处理非时变、非平稳过程的不足，提高子波提取的精度，提高进行地震层位标定的准确性和地震反演的分辨率，更准确地进行初层预测和流体识别。目前地震子波的提取方法包括两大类，第一类是确定性提取子波，第二类是统计性提取子波。第一类是指假设反射系数可有测井计算得到，结合井旁地震道，由褶积理论求出地震子波，即本文档来自技高网...

【技术保护点】
一种基于分数阶傅里叶域的子波提取方法，包括：建立地层信息初始模型；利用输入子波和地层信息初始模型的反射系数合成地震记录；计算输入子波和合成的地震记录数据的互相关矩阵以及地震记录数据的自相关矩阵，并对所述互相关矩阵和自相关矩阵进行分数阶傅里叶变换；对所述互相关矩阵和自相关矩阵应用分数阶傅里叶域的最优滤波算法来计算均方误差最小的最优阶和最优滤波算子；基于所述最优滤波算子和最优阶提取子波，其中，最优算子f(m)由以下等式表示：其中，？和？分别表示所述互相关矩阵和自相关矩阵的分数阶傅里叶变换，并且通过以下等式来提取子波？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？...

【技术特征摘要】
1.一种基于分数阶傅里叶域的子波提取方法，包括 建立地层信息初始模型； 利用输入子波和地层信息初始模型的反射系数合成地震记录； 计算输入子波和合成的地震记录数据的互相关矩阵以及地震记录数据的自相关矩阵，并对所述互相关矩阵和自相关矩阵进行分数阶傅里叶变换； 对所述互相关矩阵和自相关矩阵应用分数阶傅里叶域的最优滤波算法来计算均方误差最小的最优阶和最优滤波算子；...

【专利技术属性】
技术研发人员：洪余刚，邹文，周晶晶，彭真明，袁圆，黄东山，张洞君，陈小二，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：