发明名称 |
一种适用于逆时偏移的吸收边界条件方法 |
摘要 |
本发明是一种适用于逆时偏移的吸收边界条件方法,叠前与深度域建模后炮点波场模拟,沿时间方向二阶差分,每一个时刻t对计算区域四个水平方向分别增加n层的计算层;采用Higdon法对每个计算层循环计算,更新波场;对计算区域顶面和底面方向用二阶混合差分近似吸收边界条件公式计算,采用相同的步骤进行检波点波场模拟;利用互相关成像对炮点波场与检波点波场互相关计算得到一炮偏移结果,并完成所有炮的地震数据偏移计算,叠加形成剖面。本发明能有效处理入射角度较大的入射波的反射,在没有增加过多计算的情况下有效消除虚假反射导致的成像噪音,计算效率高、吸收效果好、以及易于实现。 |
申请公布号 |
CN104280768B |
申请公布日期 |
2017.03.15 |
申请号 |
CN201310293637.X |
申请日期 |
2013.07.12 |
申请人 |
中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司 |
发明人 |
吴晓丰;戴南浔;武威;刘磊 |
分类号 |
G01V1/28(2006.01)I;G01V1/36(2006.01)I |
主分类号 |
G01V1/28(2006.01)I |
代理机构 |
北京康信知识产权代理有限责任公司 11240 |
代理人 |
吴贵明;张永明 |
主权项 |
一种适用于逆时偏移的吸收边界条件方法,特点是通过以下步骤实现:1)采集地震数据并完成叠前处理与深度域建模;2)进行炮点有限差分波场模拟,沿时间方向的二阶差分,在每一个时刻t,针对计算区域的四个水平方向,分别增加n层的计算层;同时对于每一个计算层,循环使用基于Higdon法的一阶混合差分近似吸收边界条件公式计算,更新波场U1、U2、U3,所述波场U1、U2、U3分别是延拓过程中前一时刻、当前时刻、以及下一时刻的波场;3)针对计算区域的顶面和底面两个方向,分别增加n层的计算层,同时对于每一个计算层,循环使用基于Higdon法的二阶混合差分近似吸收边界条件公式计算,更新所述波场U1、U2、U3;所述的n取值范围和更新所述波场U1、U2、U3与步骤2)一致;4)按照步骤2)和步骤3)进行检波点有限差分波场模拟;5)利用互相关成像条件,对步骤2)、步骤3)得到的炮点波场与步骤4)得到的检波点波场进行互相关计算得到一炮偏移结果,完成所有炮的地震数据偏移计算,将结果叠加形成逆时偏移成像剖面。 |
地址 |
100007 北京市东城区东直门北大街9号中国石油大厦 |