各向异性高斯束逆时偏移方法、装置、电子设备及介质.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN115980836A(43)申请公布日2023.04.18(21)申请号202111204023.0(22)申请日2021.10.15(71)申请人中国石油化工股份有限公司地址100027北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人肖建恩李博(74)专利代理机构北京思创毕升专利事务所11218专利代理师孙向民廉莉莉(51)Int.Cl.G01V1/28(2006.01)G01V1/30(2006.01)G01V1/36(2006.01)权利要求书3页说明书10页附图5页(54)发明名称各向异性高斯束逆时偏移方法、装置、电子设备及介质(57)摘要本申请公开了一种各向异性高斯束逆时偏移方法、装置、电子设备及介质，其中偏移方法包括：通过高斯束加权积分构建格林函数；通过所述格林函数实现正向波场的延拓和反向波场的延拓；在进行波场延拓的过程中，通过相速度的运动学射线追踪方程求解高斯束，以对地下介质的各向异性对成像的影响进行校正；基于复值旅行时对地下介质的粘滞性对成像的影响进行校正；将校正后的正向波场延拓和反向波场延拓沿时间方向进行互相关，求取单炮成像结果；将所有单炮成像结果叠加，得到最终的成像结果。本发明能够有效校正地下介质中粘滞性和各向异性的影响，进一步提升成像质量，为后续储层预测提供可靠数据。CN115980836ACN115980836A权利要求书1/3页1.一种各向异性高斯束逆时偏移方法，应用于地震波场成像，其特征在于，包括：通过高斯束加权积分构建格林函数；通过所述格林函数实现正向波场的延拓和反向波场的延拓；在进行波场延拓的过程中，通过相速度的运动学射线追踪方程求解高斯束，以对地下介质的各向异性对成像的影响进行校正；基于复值旅行时对地下介质的粘滞性对成像的影响进行校正；将校正后的正向波场延拓和反向波场延拓沿时间方向进行互相关，求取单炮成像结果；将所有单炮成像结果叠加，得到最终的成像结果。2.根据权利要求1所述的各向异性高斯束逆时偏移方法，其特征在于，所述格林函数为：G(x,x0；ω)＝Φ∫u(x,x0,p,ω)dθ(1)其中，θ为出射角，Φ为初始振幅系数，u(x,x0,p,ω)为高斯束，有如下形式：其中，G为格林函数，i为虚数单位，V(x0)为地震波传播速度，A为高斯束复值振幅，t为高斯束复值走时，X为地下任一点，X0为出射点，P为慢度，ω为圆频率。3.根据权利要求2所述的各向异性高斯束逆时偏移方法，其特征在于，采用如下关系式实现反向波场的延拓：(0)其中，UB为反向延拓波场，为闭合空间表面，nx为闭合空间外法线，Sx为积分面，U为地表接收地震记录，t0为任一时刻，t为传播时间；用式(1)高斯束表征的格林函数，满足Kirchhoff近似边界条件G|z＝0＝0，替代式(3)中的G(x,t；x0,t0)，可得：其中，UB为反向延拓波场。4.根据权利要求2所述的各向异性高斯束逆时偏移方法，其特征在于，采用如下关系式实现正向波场的延拓，其中xs为震源点位置，f(ω)为地震子波的傅里叶变换，G为格林函数，UF为正向传播波场，X为地下任一点，t为传播时间。5.根据权利要求1所述的各向异性高斯束逆时偏移方法，其特征在于，所述运动学射线2CN115980836A权利要求书2/3页追踪方程的形式如下：其中P和Q为动力学射线追踪参数，M和N为射线中心坐标系下的方向，A,B,C,D为相关系数，有如下形式：6.根据权利要求2所述的各向异性高斯束逆时偏移方法，其特征在于，所述复值旅行时的公式为：*其中，tns(x)为所述复值旅行时，t(x)为声波旅行时，t(x)为沿射线路径方向的走时，ω为为圆频率，ωr为参考频率。将式(10)中的复值旅行时替代式(2)中的t，实现对粘滞性的校正。7.根据权利要求1所述的各向异性高斯束逆时偏移方法，其特征在于，通过以下关系式求取所述单炮成像结果，I(x0,xs)＝∫UB(x0,t0)UF(x0,t；xs)dt0(11)其中，I(x0,xs)为单炮成像值，UB为反向延拓波场,UF为正向传播波场，t0为任一时刻，t为传播时间，Sx为积分面。8.一种电子设备，其特征在于，所述电子设备包括：至少一个处理器；以及，与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1‑7任一所述的各向异性高斯束逆时偏移方法。9.一种非暂态计算机可读存储介质，其特征在于，该非暂态计算机可读存储介质存储计算机指令，该计算机指令用于使计算机执行权利要求1‑7任一所述的各向异性高斯束逆时偏移方法。10.一种各向异性高斯束逆时偏移装置，应