[发明专利]一种反射声波测井的波场分离方法及装置有效
申请号: | 201610945365.0 | 申请日: | 2016-10-26 |
公开(公告)号: | CN106526678B | 公开(公告)日: | 2019-08-23 |
发明(设计)人: | 宫昊;陈浩;何晓;王秀明 | 申请(专利权)人: | 中国科学院声学研究所 |
主分类号: | G01V1/50 | 分类号: | G01V1/50 |
代理公司: | 北京亿腾知识产权代理事务所(普通合伙) 11309 | 代理人: | 陈霁 |
地址: | 100190 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 实测波形 预处理 直达波 矩阵 反射波波形 反射声波 测井 反演 慢度 分离方法及装置 波场 传播 构建模块 计算模块 偏移成像 偏移模块 反射波 相干法 有效地 构建 移除 | ||
1.一种反射声波测井的波场分离方法,其特征在于,包括以下步骤:
对实测波形进行预处理,并对预处理后的实测波形中的直达波的慢度范围进行计算;
根据计算出的慢度范围求出速度范围,根据所述速度范围构建传播矩阵;
根据所述传播矩阵与所述预处理后的实测波形,反演得到所述直达波的波形;
将反演得到的直达波的波形从所述预处理后的实测波形中移除,得到反射波波形,提高所述反射波波形的信噪比,并对提高信噪比后的反射波波形进行偏移成像,获得最终的反射波波形。
2.根据权利要求1所述的方法,其特征在于,所述慢度范围通过以下步骤获取:
运用慢度时间相干法,通过对阵列波形的相关函数的计算,得到所述阵列波形的到时与慢度。
3.根据权利要求2所述的方法,其特征在于,所述相关函数的计算公式为:
其中,ρ(S,T)为相关函数,S为慢度,T为计算时窗起始位置,m为接收器的总个数,Tω为对比窗口长度,ωi(t)为第i道波形,δ为各接收器之间的间距。
4.根据权利要求1所述的方法,其特征在于,所述传播矩阵为:
其中,N是接收器的总个数,n为不大于N的自然数,L为直达波阵型的总个数,l为不大于L的自然数,ω为角频率,d为相邻两个接收器之间的距离,Vl为直达波中第l个阵型的速度。
5.根据权利要求1所述的方法,其特征在于,所述直达波的波形计算公式为:
其中,~表示复共轭,T表示转置,E为所述传播矩阵,A为所述直达波的阵型组成的矩阵,W为所述预处理后的实测波形组成的矩阵。
6.根据权利要求1所述的方法,其特征在于,所述将反演得到的直达波的波形从所述预处理后的实测波形中移除,得到反射波波形的步骤包括:
将第n个接收器对应的直达波频谱从其对应全波形的频谱中移除,得到第n个接收器的剩余信号对应的频谱,
其中,Wn(ω)是第n个接收器接收到的实测波形,N为接收器的总个数,L为直达波阵型的总个数,Al(ω)为所述直达波的第l个阵型,l为不大于L的自然数,Rn(ω)为所述第n个接收器的剩余信号对应的频谱。
7.一种反射声波测井的波场分离装置,其特征在于,所述装置包括:
计算模块,用于对实测波形进行预处理,并对预处理后的实测波形中的直达波的慢度范围进行计算;
构建模块,用于根据计算出的慢度范围求出速度范围,根据所述速度范围构建传播矩阵;
反演模块,根据所述传播矩阵与所述预处理后的实测波形,反演得到所述直达波的波形;
偏移模块,用于将反演得到的直达波的波形从所述预处理后的实测波形中移除,得到反射波波形,提高所述反射波波形的信噪比,并对提高信噪比后的反射波波形进行偏移成像,获得最终的反射波波形。
8.根据权利要求7所述的装置,其特征在于,所述计算模块具体用于:对实测波形进行预处理,并利用慢度时间相干法对预处理后的实测波形中的直达波的慢度范围进行计算。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院声学研究所,未经中国科学院声学研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201610945365.0/1.html,转载请声明来源钻瓜专利网。