基于高斯束反偏移的炮记录获取方法及系统.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110873894A(43)申请公布日2020.03.10(21)申请号201811024171.2(22)申请日2018.09.04(71)申请人中国石油化工股份有限公司地址100728北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司勘探分公司(72)发明人袁茂林杨鸿飞郑天发蒋福友王静波李苏光李彦奇肖亮(74)专利代理机构北京思创大成知识产权代理有限公司11614代理人董晓盈(51)Int.Cl.G01V1/28(2006.01)G01V1/36(2006.01)权利要求书3页说明书8页附图3页(54)发明名称基于高斯束反偏移的炮记录获取方法及系统(57)摘要公开了一种基于高斯束反偏移的炮记录获取方法及系统。该方法包括：建立地震观测系统，输入偏移结果和背景速度场，计算炮点处沿每一个方向出射的高斯束的走时、路径与地震波场；针对每一条出射的高斯束的路径，逐点判断是否存在扰动位置，若是，则计算扰动位置沿每一个方向出射的高斯束的走时和地震波场，若否，则针对下一条出射的高斯束的路径逐点判断是否存在扰动位置；根据炮点沿每一个方向出射的高斯束的走时和地震波场与扰动位置沿每一个方向出射的高斯束的走时和地震波场，计算接收点处的地震波场并进行加权叠加，获得地震炮记录。本发明通过输入偏移结果和背景速度场，计算走时、路径和复值动力学射线追踪参数，高效简便地获取地震炮记录。CN110873894ACN110873894A权利要求书1/3页1.一种基于高斯束反偏移的炮记录获取方法，其特征在于，包括：s建立地震观测系统，输入偏移结果和背景速度场，计算炮点xs处沿每一个方向p出射的高斯束的走时、路径与地震波场；针对每一条出射的高斯束的路径，逐点判断是否存在扰动位置，若是，则计算扰动位置x沿每一个方向pr出射的高斯束的走时和地震波场，若否，则针对下一条出射的高斯束的路径逐点判断是否存在扰动位置；根据炮点xs沿每一个方向出射的高斯束的走时和地震波场与扰动位置x沿每一个方向出射的高斯束的走时和地震波场，计算接收点xr处的地震波场并进行加权叠加，获得地震炮记录。2.根据权利要求1所述的基于高斯束反偏移的炮记录获取方法，其中，通过公式(1)计算出射的高斯束的走时与路径：其中，xi(s)为物理空间坐标(x,y,z)中的射线坐标，pi(s)为射线慢度矢量，τ为沿射线的走时，dτ为求积步长，V(s)为速度场沿二阶导数矩阵。3.根据权利要求1所述的基于高斯束反偏移的炮记录获取方法，其中，所述炮点xs处的地震波场为：其中，G(x,xs,ω)为炮点xs处的地震波场，ω为频率，为炮点xs处出射的高斯束的方向在x方向上的分量，为炮点xs处出射的高斯束的方向在y方向上的分量，为炮点xs处出射的高斯束的方向在z方向上的分量。4.根据权利要求3所述的基于高斯束反偏移的炮记录获取方法，其中，所述接收点xr处的地震波场为：其中，G(xr,x,ω)为接收点xr处的地震波场，为接收点xr处出射的高斯束的方向在x方向上的分量，为接收点xr处出射的高斯束的方向在y方向上的分量，为接收点xr处出射的高斯束的方向在z方向上的分量。5.根据权利要求4所述的基于高斯束反偏移的炮记录获取方法，其中，通过公式(4)进行加权叠加：u(xr,xs,ω)＝S(ω)∫voldxm(x)∫∫G(x,xs,ω)∫∫G(xr,x,ω)(4)其中，u(xr,xs,ω)为由炮点xs处激发、接收点xr处叠加形成的地震记录，m(x)为偏移结2CN110873894A权利要求书2/3页果；S(ω)为炮点函数。6.一种基于高斯束反偏移的炮记录获取系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：s建立地震观测系统，输入偏移结果和背景速度场，计算炮点xs处沿每一个方向p出射的高斯束的走时、路径与地震波场；针对每一条出射的高斯束的路径，逐点判断是否存在扰动位置，若是，则计算扰动位置x沿每一个方向pr出射的高斯束的走时和地震波场，若否，则针对下一条出射的高斯束的路径逐点判断是否存在扰动位置；根据炮点xs沿每一个方向出射的高斯束的走时和地震波场与扰动位置x沿每一个方向出射的高斯束的走时和地震波场，计算接收点位置xr处的地震波场并进行加权叠加，获得地震炮记录。7.根据权利要求6所述的基于高斯束反偏移的炮记录获取系统，其中，通过公式(1)计算出射的高斯束的走时与路径：其中，xi(s)为物理空间坐标(x,y,z)中的射线坐标，pi(s)为射线慢度矢量，τ为沿射线的走时，dτ为求积步长。8.根据权利要求6所述的基于高斯束反偏移的炮记录获取系统，其中，所述炮点xs处的地震波场为：其中，G