[发明专利]近地表Q偏移方法及装置有效
申请号: | 201811129861.4 | 申请日: | 2018-09-27 |
公开(公告)号: | CN108919354B | 公开(公告)日: | 2019-09-27 |
发明(设计)人: | 徐锦承;刘礼农;张剑锋 | 申请(专利权)人: | 中国科学院地质与地球物理研究所 |
主分类号: | G01V1/30 | 分类号: | G01V1/30 |
代理公司: | 北京超凡志成知识产权代理事务所(普通合伙) 11371 | 代理人: | 梁斌 |
地址: | 100000 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 近地表 偏移成像 偏移速度场 目标成像 扫描成像 速度场 偏移 检波 叠前地震数据 场计算 成像点 基准面 叠加 衰减 吸收 | ||
1.一种近地表Q偏移方法,其特征在于,包括:
获取目标工区的叠前地震数据;
根据所述叠前地震数据确定基准面,并建立所述基准面的上下两个偏移速度场,包括:近地表等效速度场和基准面以下偏移速度场;
根据所述偏移速度场确定出近地表Q值估计的目标成像线,所述基准面是连续光滑的时间深度面;
根据所述目标成像线处计算所述目标工区的目标检波点的常近地表Q补偿的常Q扫描成像道集;
根据所述常Q扫描成像道集计算对应的所述目标检波点的位置的部分近地表等效Q值,其中,一个检波点组对应一个近地表等效Q值,相邻检波点组成一个检波点组;
根据所述部分近地表等效Q值和所述近地表等效速度场建立所述目标工区的目标近地表等效Q值场;
根据所述近地表等效速度场和所述目标近地表等效Q值场计算所述目标工区各个成像点对应的补偿近地表吸收衰减后的偏移成像道集;
将所述偏移成像道集叠加形成所述目标工区的偏移成像剖面。
2.如权利要求1所述的近地表Q偏移方法,其特征在于,所述建立所述基准面的上下两个偏移速度场,包括:
将所述叠前地震数据利用初始叠加速度场计算各个成像点的成像道集;
通过反动校正将所述各个成像点的成像道集拉平以得到所述基准面以下的偏移速度,经过平滑处理得到基准面以下偏移速度场;
根据所述基准面以下偏移速度场按设定的百分比扫描近地表等效速度生成局部时窗范围的成像道集,根据反射同相轴聚焦程度选取百分比,得到近地表等效速度场。
3.如权利要求1所述的近地表Q偏移方法,其特征在于,所述根据所述目标成像线处计算所述目标工区的部分检波点的常近地表Q补偿的常Q扫描成像道集的步骤,包括:
根据所述近地表等效速度场的分布特征选定目标检波点;
将所述目标检波点利用共检波点组成像道集计算公式计算得到常Q扫描成像道集。
4.如权利要求3所述的近地表Q偏移方法,其特征在于,所述共检波点组成像道集计算公式表示为:
I(x,y,T)=(Ag/As)g(τs+τg,η),
其中,x和y表示成像点的坐标,T表示成像点的时间深度,n表示所述叠前地震数据包含的地震道总数目,以i标记该n个地震道中的第i个地震道,对于该第i个地震道,Ag表示与检波点位置相关的幅值,As表示与炮点位置相关的幅值,τg表示检波点到成像点的走时,τs表示炮点到成像点的走时,η表示与射线参数相关的参数,xg和yg表示检波点的坐标,Ig(x,y,T,xg,yg)表示共检波点成像道集,I(x,y,T)表示近地表等效Q补偿成像,g(τs+τg,η)表示地震道进行近地表Q补偿计算后的时间序列。
5.如权利要求1或4所述的近地表Q偏移方法,其特征在于,根据所述常Q扫描成像道集计算对应的所述目标检波点的位置的部分近地表等效Q值的步骤,包括:
根据所述常Q扫描成像道集选择有近地表Q值标定附近的共检波点成像道集;
确定一目标时窗,计算该目标时窗的频谱曲线的包络作为参考频谱;
根据所述参考频谱和其它近地表Q值扫描的共检波点成像道集的频谱,选择与参考频谱接近的近地表Q值作为所述目标检波点组对应的部分近地表等效Q值。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院地质与地球物理研究所,未经中国科学院地质与地球物理研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811129861.4/1.html,转载请声明来源钻瓜专利网。