[发明专利]一种频率域正演方法及装置

摘要

本发明涉及地质勘探技术领域，具体的讲是一种频率域正演方法及装置。其中方法包括，建立17点格式的差分公式，构造稀疏矩阵；读入子波参数和速度模型；频率循环得到单频波场；对所有频率波场求反傅里叶变换；得到正演结果。通过本发明实施例的方法及装置，四阶高精度正演克服了以往对dx必须等于dz的条件限制，不仅可应用于dx=dz,也可以适用于dx≠dz的情况，更适合在实际生产中运用，而且最小波长所需要的网格数在该方法中仅需要2.4个，也优于其他类似方法。

主权项

1.一种频率域正演方法，其特征在于包括，建立如下的17点格式的差分公式，其中P‾m+1,n=α2Pm+1,n+1-α22(Pm+1,n-1),]]>P‾m-1,n=α2Pm-1,n+1-α22(Pm-1,n+1+Pm-1,n-1)---(3)]]>P‾m+2,n=α1Pm+2,n+1-α12(Pm+2,n+2+Pm+2,n-2),]]>P‾m-2,n=α1Pm-2,n+1-α12(Pm-2,n+2+Pm-2,n-2),]]>P‾m,n=α3Pm,n+α4(Pm,n+1+Pm,n-1)+1-2α4-α32(Pm,n+2+Pm,n-2),]]>和P~m,n+1=β2Pm,n+1+1-β22(Pm+1,n+1+Pm-1,n+1),]]>P~m,n-1=β2Pm,n-1+1-β22(Pm+1,n-1+Pm-1,n-1),---(4)]]>P~m,n+2=β1Pm,n+2+1-β12(Pm+2,n+2+Pm-2,n+2),]]>P~m,n-2=β1Pm,n-2+1-β12(Pm+2,n-2+Pm-2,n-2),]]>P~m,n=β3Pm,n+β4(Pm+1,n+Pm-1,n)+1-2β4-β32(Pm+2,n+Pm-2,n),]]>其中P_m,n≈P(mΔz,nΔx)，Δx,Δz代表x方向采样间隔和z方向采样间隔，m和n代表z方向的网格坐标和x方向的网格坐标，α_i,β_i,b,c,d,e,f代表求导时候的加权系数；构造稀疏矩阵；读入子波参数和速度模型；频率循环得到单频波场；对所有频率波场求反傅里叶变换；得到正演结果。