[发明专利]振幅随偏移距变化道集获取方法和装置

权利要求书

1.一种振幅随偏移距变化道集获取方法，其特征在于，地球介质为各向异性介质，所述各向异性介质包括n+1层的层状介质，所述n为大于1的整数，所述方法包括：

使地震波从所述n+1层的层状介质中的第n+1层传播到第1层；

建立所述地震波在所述第n+1层中的波位函数，所述波位函数包括入射纵波的波位函数、入射横波的波位函数、反射纵波的波位函数和反射横波的波位函数；

根据所述入射纵波的波位函数、所述入射横波的波位函数、所述反射纵波的波位函数、所述反射横波的波位函数以及所述n+1层的层状介质的边界条件函数建立第一矩阵，所述第一矩阵表示所述地震波在每两层之间分界面上的位移分量与应力分量；

根据所述第一矩阵建立表示所述n+1层的层状介质中第n层地震参数与第n-1层的地震参数的关系的第二矩阵，所述地震参数包括位移分量和应力分量；

根据所述第二矩阵建立表示所述n+1层的层状介质中所述第n层的地震参数与第1层的地震参数的关系的第三矩阵；

根据所述第三矩阵获取纵波反射系数和横波反射系数；

根据测井曲线数据确定子波；

根据所述纵波反射系数、所述横波反射系数和所述子波获取振幅随偏移距变化AVO道集；

其中，所述第一矩阵为：

G⁽ⁿ⁾＝(S_cd)⁽ⁿ⁾[M₁ⁿ,N₁ⁿ,M₂ⁿ,N₂ⁿ]^Te^j(ωt-σx)，

其中，所述所述n表示第n层，所述M₁ⁿ为所述第n层的入射纵波的振幅，N₁ⁿ为所述第n层的入射横波的振幅，M₂ⁿ为所述第n层的反射纵波的振幅，N₂ⁿ为所述第n层的反射纵波的振幅，所述j为虚数符号，所述ω为角速度，所述t为时间，所述σ为波矢量，所述x为x轴上的位置坐标，所述u⁽ⁿ⁾为所述地震波在所述第n层中沿x轴方向上的位移分量，所述w⁽ⁿ⁾为所述地震波在所述第n层中沿z轴方向上的位移分量，所述σ_zz⁽ⁿ⁾为所述地震波在所述第n层中沿z轴方向上的应力分量，所述τ_xz⁽ⁿ⁾为所述地震波在所述第n层中沿x轴方向上的应力分量，

所述所述角标c表示矩阵的行序号，所述角标d表示矩阵的列序号，其中，

所述φ_iP1为所述入射纵波的偏振角，所述φ_is1为所述入射横波的偏振角，所述φ_p1为所述反射纵波的偏振角，所述φ_S1为所述反射横波的偏振角，所述F₁为入射纵波数的垂直分量，所述F₁ⁿ为第n层内入射纵波数的垂直分量，所述s₁为入射横波数的垂直分量，所述s₁ⁿ为第n层内入射横波数的垂直分量，所述F₂为反射纵波数的垂直分量，所述s₂为反射横波数的垂直分量，所述F₂ⁿ为第n层内反射纵波数的垂直分量，所述s₂ⁿ为第n层内反射横波数的垂直分量，所述z为z轴上的位置坐标，所述z轴垂直于所述层状介质且指向所述层状介质的第1层，所述c₁₃、所述c₁₅、所述c₃₃、所述c₃₅和所述c₅₅均为各向异性介质的弹性系数，所述σ为波矢量。