[发明专利]基于深度波数谱的瞬变波包提取方法、存储介质与终端有效
申请号: | 201910040925.1 | 申请日: | 2019-01-16 |
公开(公告)号: | CN109856672B | 公开(公告)日: | 2019-09-20 |
发明(设计)人: | 张繁昌;胥勋勇;朱磊;桑凯恒;唐杰;张佳佳 | 申请(专利权)人: | 中国石油大学(华东) |
主分类号: | G01V1/28 | 分类号: | G01V1/28 |
代理公司: | 济南圣达知识产权代理有限公司 37221 | 代理人: | 李琳 |
地址: | 266580 山*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 波数谱 积分变换 波包 瞬变 存储介质 地震记录 深度域 波数 终端 低通滤波 地震信号 反变换 复赛谱 观测 合成 返回 | ||
1.一种基于深度波数谱的瞬变波包提取方法,其特征是:包括:
对深度域地震信号进行积分变换,得到深度波数谱;
将每一深度点的波数谱变换到同态域,得到对应的波数复赛谱,经低通滤波后反变换为波数谱;
确定波数谱中每一深度点的波数谱最大值,根据该最大值对应的波数,计算对应深度的瞬变波包;
合成深度域地震记录,并与观测地震记录对比,若差别超过设定值则返回积分变换,改变积分变换的调节因子,重新得到深度波数谱,直到最终计算的差值满足设定值。
2.如权利要求1所述的一种基于深度波数谱的瞬变波包提取方法,其特征是:对深度域地震信号进行积分变换的具体过程为:
式中S(γ,ω)表示深度波数谱;γ表示目标深度点;ω表示圆波数;s(h)表示深度域地震信号;h表示对应的深度范围;W(h,ω)是深度域高斯窗函数:
其中β为调节因子。
3.如权利要求1所述的一种基于深度波数谱的瞬变波包提取方法,其特征是:波数谱变换到同态域,得到其对应的波数复赛谱的具体过程包括:对每一深度的波数谱进行傅里叶变换,对变换结果的对数值进行傅里叶反变换,进而得到每一深度的波数复赛谱。
4.如权利要求1所述的一种基于深度波数谱的瞬变波包提取方法,其特征是:计算瞬变波包的具体过程为:在经过低通滤波并反变换后的波数谱中,找到每一深度点的波数谱最大值,此最大值对应的波数即为该深度点波包的主波数,将主波数带入雷克子波公式得到该深度的瞬变波包。
5.如权利要求1所述的一种基于深度波数谱的瞬变波包提取方法,其特征是:合成深度域地震记录,并与观测地震记录对比的具体过程为:将得到的瞬变波包带入到瞬变褶积模型中,得到深度域合成地震记录,得到的深度域合成地震记录与观测地震记录进行对比。
6.如权利要求5所述的一种基于深度波数谱的瞬变波包提取方法,其特征是:瞬变褶积模型为:
式中Y表示长度为N的深度域合成地震记录;w1、w2,…,wN表示每一深度点的瞬变波包,这些波包构成波包矩阵;r1、r2,…,rN表示测井反射系数。
7.一种计算机可读存储介质,其特征是:其中存储有多条指令,所述指令适于由终端设备的处理器加载并执行权利要求1-6中任一项所述的基于深度波数谱的瞬变波包提取方法。
8.一种终端设备,其特征是:包括处理器和计算机可读存储介质,处理器用于实现各指令;计算机可读存储介质用于存储多条指令,所述指令适于由处理器加载并执行权利要求1-6中任一项所述的基于深度波数谱的瞬变波包提取方法。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(华东),未经中国石油大学(华东)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910040925.1/1.html,转载请声明来源钻瓜专利网。