[发明专利]一种成像方法和装置有效
申请号: | 201610971350.1 | 申请日: | 2016-10-28 |
公开(公告)号: | CN107678058B | 公开(公告)日: | 2019-06-11 |
发明(设计)人: | 叶月明;李立胜;刘午牛;王启迪;邵萌珠 | 申请(专利权)人: | 中国石油天然气股份有限公司 |
主分类号: | G01V1/28 | 分类号: | G01V1/28 |
代理公司: | 北京三友知识产权代理有限公司 11127 | 代理人: | 李辉;徐焕 |
地址: | 100007 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 成像结果 地层 成像 补偿因子 倾角补偿 波场 方法和装置 起伏地表 采样点 累加 累加结果 预设参数 散射 屏蔽 地下 传播 分析 | ||
1.一种成像方法,其特征在于,包括:
获取待测地层多个采样点在倾角域多个倾角中各个倾角的成像结果;
获取所述待测地层多个采样点在倾角域多个倾角中各个倾角的补偿因子;
根据所述成像结果、所述补偿因子和预设参数,确定所述多个倾角中各个倾角补偿后的成像结果;
对所述多个倾角中各个倾角补偿后的成像结果进行累加,将累加结果作为所述待测地层的成像结果;
其中,获取所述待测地层多个采样点在倾角域多个倾角中各个倾角的成像结果,包括:依次获取所述多个采样点中各个采样点的第一震源方向波场和第一检波器方向波场,其中,所述第一检波器方向波场是检波器位置不存在震源时的波场;按照以下方式对所述多个采样点中各个采样点执行以下操作,以得到各个采样点的第一成像矩阵:将当前采样点的第一震源方向波场和第一检波器方向波场变换至频域,并将变换至频域的波场进行互相关成像,得到所述当前采样点的第一成像矩阵;将所述多个采样点中各个采样点的第一成像矩阵变换至倾角域,得到所述多个倾角中各个倾角的成像结果;
获取所述待测地层多个采样点在倾角域多个倾角中各个倾角的补偿因子,包括:依次获取所述多个采样点中各个采样点的第二震源方向波场和第二检波器方向波场,其中,所述第二检波器方向波场是检波器位置存在震源时的波场;按照以下方式对所述多个采样点中各个采样点执行以下操作,以得到各个采样点的第二成像矩阵:将当前采样点的第二震源方向波场和第二检波器方向波场变换至频域,并将变换至频域的波场进行互相关成像,得到所述当前采样点的第二成像矩阵;将所述多个采样点中各个采样点的第二成像矩阵变换至倾角域,得到所述多个倾角中各个倾角的补偿因子。
2.如权利要求1所述的方法,其特征在于,按照以下公式对所述多个倾角中各个倾角补偿后的成像结果进行累加:
其中,image(x,z+Δz)表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角域多个倾角中各个倾角补偿后的成像结果进行累加的结果,LiFinal(xi,zi+Δz)表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角为i时进行补偿后的成像结果,θmin≤i≤θmax,其中,θmin表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角域多个倾角的最小值,θmax表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角域多个倾角的最大值,Δz表示预设步长。
3.如权利要求1所述的方法,其特征在于,按照以下公式确定所述多个倾角中各个倾角补偿后的成像结果:
LiFinal(xi,zi+Δz)=Lia(xi,zi+Δz)/(Fia(xi,zi+Δz)+ε)
其中,LiFinal(xi,zi+Δz)表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角为i时进行补偿后的成像结果,Lia(xi,zi+Δz)表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角为i时的成像结果,Fia(xi,zi+Δz)表示所述待测地层中横向位置为x、深度为z+Δz的采样点在倾角为i时的补偿因子,ε表示所述预设参数,Δz表示预设步长。
4.如权利要求3所述的方法,其特征在于,所述预设参数在0.0001到0.1之间取值。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气股份有限公司,未经中国石油天然气股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201610971350.1/1.html,转载请声明来源钻瓜专利网。